
Calhoun 

imU'uiioiu! A Kim of die Nav.il Po<ic^rdduacc School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 



Theses and Dissertations 


1. Thesis and Dissertation Collection, all items 


1994-12 

Three-dimensional computer graphics 
visualization of target detection 

Gorgulu, Mehmet; Yilmaz, Mustafa 

Monterey, California. Naval Postgraduate School 


http://hdl.handle.net/10945/42804 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 

Downloaded from NPS Archive: Calhoun 



DUDLEY 

KNOX 

LIBRARY 


http ://w w w. nps.edu/lEbrary 


Calhoun is the Naval Postgraduate Sc hoof's public access digital repository for 
research mate hate and institutional publications created by tire NPS community, 
Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
appointed — and published — scholar^ author. 

Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Cirde 
Monterey, California USA 93943 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



THREE-DIMENSIONAL COMPUTER GRAPHICS 
VISUALIZATION OF TARGET DETECTION 

by 

Mehmet Gorgulu 
Mustafa Yilmaz 
December 1994 

Thesis Advisor: Richard Chris Olsen 

Co-Advisor : David R. Pratt 


Approved for public release; distribution is unlimited. 







REPORT DOCUMENTATION PAGE 


Form Approved OMB No. 0704 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing 
instruction, searching existing data sources, gathering and maintaining the data needed, and completing and reviewing the collection of 
information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions 
for reducing this burden, to Washington headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson 
Davis Highway, Suite 1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 

1 . AGENCY USE ONLY (Leave blank) 

2. REPORT DATE 

December 1994 

3. REPORT TYPE AND DATES COVERED 
Master’s Thesis 

4. TITLE AND SUBTITLE Three-Dimensional Computer Graphics Visualization Of 
Target Detection 

UNCLASSIFIED 

5. FUNDING NUMBERS 

6. AUTHOR(S) 

Mehmet Gorgtilii - Mustafa Yilmaz 



7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Naval Postgraduate School 

Monterey CA 93943-5000 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

10. SPONSORING/MONITORING 
AGENCY REPORT NUMBER 

11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official policy or 
position of the Department of Defense or the U.S. Government. 

12a. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; 
distribution unlimited 

12b. DISTRIBUTION CODE 

13. ABSTRACT The purpose of this thesis is to visualize the sensor performance for a generic missile. We simulate the proceses 
performed by a missile using IR or TV sensors. Two generic scenes (background) were created, one for each generic sensor. The 
program simulates the scene from the point of view of a missile sensor. A graphical user interface was included for user input. These 
inputs provide the initial environmental conditions and the structural specifications of the sensor and the targets. Depending on these 
inputs, the sensor will show a detection and a lock-on range to the user. The detection range for the IR sensor was based on the intensity 
of the signal, above a specific threshold. For the TV system, target contrast was used. Atmospheric extinction was included. Several 
aspects of the SGI hardware and software capability were used to mimic physical problems and processes at considerable savings in 
computational effort. One was the use of the SGI Gouraud shading capability to establish the temperature distribution for IR targets; a 
second was use of the hardware (screen) projection to map from 3-D to 2-D. For further work, this program can be integrated to the 
EOTDA (Electro-optical Tactical Decision Aid) software. The graphics part of the program was written by using OpenGL graphics 
library and the user interface was implemented by using OSF/Motif The main program was implemented in C++ on Silicon Graphics 
Reality Engines. 

14. SUBJECT TERMS Infrared, TV, detection, lock-on, target 




15. NUMBER OF PAGES 

178 







16. PRICE CODE 

17. SECURITY 

CLASSIFICATION OF 

REPORT 

18. SECURITY 

CLASSIFICATION OF THIS 
PAGE 

19. SECURITY 

CLASSIFICATION OF 

ABSTRACT 

20. LIMITATION OF 

ABSTRACT 

UL 

Unclassified 

Unclassified 

Unclassified 




NSN 7540-01-280-5500 


Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 






























Approved for public release; distribution is unlimited. 


Three-Dimensional Computer Graphics Visualization of Target Detection 


by 

Mehmet Gorgulu 

B.S., Turkish Naval Academy , 1988 
Mustafa Yilmaz 

B.S., Turkish Naval Academy, 1988 

Submitted in partial fulfillment 
of the requirements for the degree of 

MASTER OF SCIENCE IN APPLIED PHYSICS 

from the 

NAVAL POSTGRADUATE SCHOOL 
December 1994 


Authors: 


Approved by: 



Richard Chris Olsen, Thesis Advisor 


David R. Pratt, Co-Advisor 


#,8C988i©n i'o? 
"”IiH QEA&I 

Unannounced 
1 Justification. 


William R. Colson, GHairman, Department of Physics 

V 


| Distribution/_ 

f Avaiiabllit's. Coflss_ 
I L&vail. siuct/w? 

. I Special. 


»: mm g •mr- 














ABSTRACT 


The purpose of this thesis is to visualize the sensor performance for a 
generic missile. We simulate the processes performed by a missile using IR or 
TV sensors. Two generic scenes (background) were created, one for each generic 
sensor. The program simulates the scene from the point of view of a missile 
sensor. A graphical user interface was included for user input. These inputs 
provide the initial environmental conditions and the structural specifications of the 
sensor and the targets. Depending on these inputs, the sensor will show a 
detection and a lock-on range to the user. The detection range for the IR sensor 
was based on the intensity of the signal, above a specific threshold. For the TV 
system, target contrast was used. Atmospheric extinction was included. Several 
aspects of the SGI hardware and software capability were used to mimic physical 
problems and processes at considerable savings in computational effort. One was 
the use of the SGI Gouraud shading capability to establish the temperature 
distribution for IR targets; a second was use of the hardware(screen) projection to 
map from 3-D to 2-D. For further work, this program can be integrated to the 
EOTDA (Electro-optical Tactical Decision Aid) software. The graphics part of the 
program was written by using OpenGL graphics library and the user interface was 
implemented by using OSF/Motif. The main program was implemented in C++ on 
Silicon Graphics Reality Engines. 
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I. INTRODUCTION 


A. A BRIEF DESCRIPTION OF OUR WORK 

Visual effects have increased in importance in the last decade in both program¬ 
ming and interacting with computers. With the help of rapidly improving technology, 
we can put highly detailed graphics on our computer screens instead of only text char¬ 
acters. It has been shown that human perception is very high to visually presented 
data. As a result, computer graphics and visualization of data have been gaining an in¬ 
creasing importance in representing data visually. [Ref. 3] 

Military platforms need sensors to perform their functions. These sensors are 
expensive to build and their performance is highly dependent on the environmental 
conditions and sensor characteristics. Any change in these characteristics changes the 
performance of sensors, and the ability of a sensor to perform its targeting function. 
Our work simulates this process and allows the study of the sensor performance at low 
cost, while varying environmental conditions. 

For our work, we chose generic targets and backgrounds. We began with the 
physics problem of calculating the atmospheric absorption depending on the range, an¬ 
gle and wavelength and target projection area. We calculated the radiation for each 
node of the target and then found the power radiated by the target. By applying at¬ 
mospheric extinction to this signal, we computed the signal level at the sensor. The 
temperature level at the sensor determined the temperature distribution on the target. 
The contrast value between the target and the background established the detection 
criteria for the target. There are different detection algorithms for the detection of the 
targets. These detection algorithms were examined and cross-box technique was ap¬ 
plied in the program. Different computer graphics techniques were used in the pro¬ 
gram to be able to visualize the IR and TV scene and to be as close to reality as 
possible. 
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B. CHAPTER LAYOUT 

Chapter II discusses the background and previous work in this area. It talks 
about the theory behind the modeling, different computer graphics and target visualiza¬ 
tion methods and the atmospheric extinction coefficient. The calculation of the atmos¬ 
pheric extinction coefficient and the use of Lowtran 6 including the reasons are 
described in this part. 

Chapter III defines the problem and presents our assumptions. Chapter IV is 
an examination of the design decisions and the tools. It also covers the data structures 
for IR and TV sensors. Chapter V talks about the logical flow of the program and the 
components of the program, some of the problems and results. Chapter VI outlines the 
achievements in the program and makes suggestions on the future work of the pro¬ 
gram. Supplementary and detailed information about the program is presented in the 
Appendices. 
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II. PREVIOUS AND RELATED WORK 


A. INTRODUCTION 

A basic military problem in warfare is determining engagement ranges for 
weapon systems. Engagement is dependent on detection and recognition. The factors 
which affect the engagement range are environmental conditions, sensor characteris¬ 
tics, and the platforms which the sensor is on. 

One tool which is used to prepare for engagement is the EOTDA (Electro- 
Optical Tactical Decision Aid) Software. This software has three sensor elements: 
Infra-red, TV and laser. There is a user interface to the program and data is provided 
to the program by this interface. This data includes the sensor type, environmental 
data and target type and target information. The program calculates the sensor lock- 
on range. The calculated data is displayed on the screen. The advantage of this pro¬ 
gram is it gives a good estimate of the detection and the lock-on range. The disadvan¬ 
tage is it was implemented on a PC. Hughes Aircraft Corporation implemented this 
program in Unix environment, but our brief Experience showed the program was not 
transportable. 

We planned to work with the EOTDA software, and design a 3-D visualization 
of the EOTDA output. Unfortunately, the software is not allowed to be released to 
non-U. S. officers (allies). In line with our original intention, we decided to build some 
generic targets and backgrounds, and then work on detection ranges with generic IR 
and TV sensors. We modeled a tank, an aircraft, a ship and a building target with four 
possible backgrounds: soil, grass, concrete and sea. 

We designed the program so that the user could choose the initial conditions 
for the program and then set the environment. The purpose of the user is to determine 
the detection and lock-on range for the IR and the TV system. In the program, the 
user was given the full control for changing the distance from the sensor to the target. 
As the distance changes, the calculations are done again depending on the changing 
distance. If the signal coming from either the background or the target exceeds the 
threshold level, the program gives detection. 
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The above implementation resulted in a computer program that can measure 
and visualize the performance of a generic IR or TV missile system depending on the 
input data. (The measure of performance for a missile system is detection distance.) 
The reliability of the result is completely dependent on the input data. The missile sys¬ 
tem gets the input and then tries to distinguish the target from the background and 
whenever there is detection, this is presented on the screen. The IR sensor part of the 
program we have, gives an idea of the temperature distribution on the target and the 
background. The user can see the temperature distribution on the target from differ¬ 
ent angles. The calculation method and the detection algorithm of the program rely on 
physical algorithms. 

B. THEORY 

I. Heat And Temperature 

Heat can be defined as thermal energy in transition. It flows from one place to 
another as a result of temperature difference. The flow of heat changes the energy lev¬ 
els in the objects. Temperature is a property of the object and not a complete measure 
of the internal energy of the matter. Heat always flows from the object which is at the 
higher temperature to the object which is at the lower temperature. This is the princi¬ 
ple governing the heat transfer between a target (heat source) and a sensor. [Ref. 14] 
There are three models of heat transfer: conduction, convection and radiation. 
All heat transfer processes occur by one or more of these three modes. Infrared meas¬ 
urement is based on the radiative heat flow and is most closely related to the radiation 
mode of heat transfer. [Ref. 14] 

Radiative heat transfer can take place across a vacuum and it occurs by electro¬ 
magnetic emission and absorption. It occurs at the speed of light, and the energy 
transferred is proportional to the difference of the fourth power of the temperature. 

The bulk of the transport takes place in the infrared portion of the spectrum, from 0.75 
to 100 pm. [Ref. 14] 
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2. Infrared Radiation 

The measurement of thermal infrared radiation is the basis for non-contact tem¬ 
perature measurement and thermography. The thermal infrared radiation leaving a sur¬ 
face is called exitance (radiosity). 

Every object in nature emits wavelengths corresponding to its temperature. It 
is easy to get a temperature signature of an emitting object. The IR sensors are less 
susceptible to atmospheric conditions, and the frequency which the terrestial objects 
emit generally falls into IR frequency range. These are some of the fundamental rea¬ 
sons why the IR window is used. 

All targets radiate energy in the infrared spectrum. The hotter the target, the 
more energy radiated. Very hot targets radiate in the visible spectrum. Infrared detec¬ 
tors can sense the infrared radiant energy and produce useful electrical signals propor¬ 
tional to the temperature of target surfaces. Every target surface above absolute zero 
(0 K or -273 °C) radiates energy in the infrared. When the targets are hot enough, they 
radiate or glow in the visible part of the spectrum. 

As surfaces cool, not only do they emit less energy, but the wavelength distri¬ 
bution shifts to longer infrared wavelengths. Even though the human eye is incapable 
of sensing this energy, infrared sensors can sense these invisible longer wavelengths. 
They enable us to measure the self-emitted radiant energy from even very cold targets. 
[Ref. 14] 

A blackbody radiator is an idealized source of radiation. Blackbody radiation 
has two important characteristic: Total energy, and the shape of the curve of intensity 
vs. wavelength (Location of peak.) Blackbody radiation is perfectly diffuse and radi¬ 
ates at all wavelengths. A blackbody source has an emissivity of 1.0 but the real 
sources in the nature have varying emissivities. In our program , we accepted our tar¬ 
gets as blackbodies and made our calculations accordingly. 

Every blackbody has a spectral emittance curve. Figure 1 shows a family of 
blackbody spectra, varying with temperature. These spectral curves were explained by 

Plank, at the turn of this century. [Ref 13] Equation 1 gives the power radiated per 
unit surface area per unit wavelength, as a function of X. 
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M x(T) =Watts-cm- 3 ; (1) [Ref.7] 
6 1 


In this equation h is Plank's constant, c is the speed of the light, k is the Boltz- 
man constant, and X is the wavelength at which the blackbody radiates energy. Equa¬ 
tion 1 can be integrated to give the total radiation in a given wavelength range when 
integrated over all X we obtain 

W = J 2 W^ = 5eT4 ( 2 ) [ Ref7 J 

This is Stefan-Boltzman Law, obtained empirically by Stefan in 1879. The peak 
wavelength of the radiated energy is dependent of the temperature and is given by 
Wien's displacement law: 

X m = b/T (3) [Ref. 5] 

In the above Equations : 

W= Radiant Power emitted per unit area 
e = Emissivity 

8 = Stephan - Boltzman cons tan t 
T= Absolute temperature of the target (K) 

X m = Wavelength of maximum radiation 
b = Wien's displacement constant = 2897 K 

The sensor design we implemented responds to 8-12 Jim bandpass. The proper 
approach would be to integrate Equation 1 over the X range. Because the integration 

process slows the program down substantially, we assumed that the majority of the en¬ 
ergy was 8-12 Jim band, and simply used the Stefan-Boltzman Low to calculate the to¬ 
tal radiance of the targets. Note that for terrestial temperature ranges, the radiation in 

any IR band varies linearly with temperature [Figure 1, and also Ref 13, page 121.] 

The result of Equation 2 gives the radiation at zero range. The radiance we 

can use is the radiance normal to the surface of the target and this is found by dividing 
total radiant energy by it. So, N=W/7t, N being the irradiance of the surface. Then 
the irradiance of the surface at a given range R from the object is: 

Pr = ^ “j 2 - W/m 2 (4) [Ref. 5] 

In this equation, R is the distance between the target and the sensor, and \ is 
the projectional area of the target on the sensor. The energy at the sensor is 
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calculated by multiplying the energy at the target by the atmospheric transmission fac¬ 
tor. In our program, we did not assume that the effect of atmospheric extinction was 
simply inversely quadratic. We found the atmospheric extinction by using the Lowtran 
6 code described below. This software gave us the atmospheric extinction coefficient. 

The energy which is sensed by the sensor is also dependent on the projection 
area of the target on the sensor. The reason for target area dependency and how to 
calculate the projection area can be understood from Figure 2. The equation to calcu¬ 
late the projection area of the target is as in Equation 5. 

A T =l*h*cos0*cos<j)+w*h*cos0*sin<|)+l*w*sin0 (5) [Ref. 5] 

In this equation Aj. is the projected area of the target, 0 is the elevation angle, 
<|> is the azimuth angle, 1 is the length of the target, w is the width of the target, h is the 

height of the target. 

3. Visual Wavelength Signals 

In daylight reflected solar energy in the visible wavelength, can be used by TV 
sensors to detect the targets. Our implementation of TV detector depends on a con¬ 
trast model. In this sense, contrast means the difference between the gray level of a 
target and the background. There is a sudden change of gray levels in transition from 
the environment to the target. 

The energy coming from the target and the environment produces different illu¬ 
mination on the TV sensor screen. The energy striking the sensor plate induces differ¬ 
ent voltage levels for target and background. An electron beam scans this plate and 
creates different illumination values on the screen for each pixel as shown in Figure 3. 
Atmospheric factors affect this type of detection more than it does to other types of 
detections. We simulated this by applying fog in our program using an extinction 
model described below. 

The software we used (OpenGL) provides us three types of fog functions. A 

fog function can be linear, quadratic or exponential. Since the atmospheric extinction 
occures exponentially, we chose the exponential fog function e _dR . In this equation d 
is the fog density coefficient which corresponds to the atmospheric extinction coeffi¬ 
cient in the IR part and R is the distance in eye coordinates from origin to the object. 
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TV tubes are a kind of electron device which converts an optical image into an 
electrical signal. These tubes are used to generate a train of electrical pulses which 
represent light intensities present in an optical image focused on the tube. These inten¬ 
sities correspond to the luminance or gray levels on the screen. Each picture element 
(pixel) on the screen has a luminance value. 

These tubes use an electron beam to scan a photoconductive target which is the 
light sensor. A transparent conductive layer applied to the front side of the photocon¬ 
ductor serves as the signal (target) electrode. When a light pattern is focused on the 
photoconductor, its conductivity increases in the illuminated areas and the back side of 
the target charges to more positive values. The electron beam then reads the signal by 
depositing electrons on the positively charged areas thereby providing a capacitively 
coupled signal at the signal electrode. [Ref. 7] 

The scene radiance is represented on the screen with different gray levels. The 
reflection of the scene on the screen will be represented by the pixels which have dif¬ 
ferent brightness (luminosity) values. So there is going to be a pattern of contrast on 
the screen as shown in Figure 4 . The contrast among the pixels is related with their 
gray levels, i.e., with their luminosity. As a result of this there will be a distribution of 
contrast on the screen among the pixels. This contrast can be expressed in several 
ways. We define it as the difference of the maximum brightness and the minimum 
brightness. C=GL max -GL min where GL max is the maximum gray level, GL min is the mini¬ 
mum gray level. [Ref. 6], More information about the detection algorithm of TV sen¬ 
sors can be found in Chapter V. 

4. Extinction Coefficient 

The performance of military systems for imaging, target detection, tracking, 
target designation, and warning is strongly dependent on the transmission of the me¬ 
dium. The important atmospheric effects to be dealt with are refraction, absorption 
and scattering by the molecular constituents of the atmosphere. [Ref. 5] 

The transmission characteristics of the medium in the measurement path be¬ 
tween the target and the detector need to be considered in making noncontact meas¬ 
urements. No loss of energy is encountered when measuring through a vacuum. The 
atmosphere, however does not allow transmission in every band. As seen in the Figure 
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5, there are two spectral intervals which have high transmission. These are located in 
the 3-5 (im and 8-12 pm wavelength range. The IR sensors are built to scan in one of 
these atmospheric windows. 

The effect of the atmosphere on the signal is represented by the extinction coef¬ 
ficient, which is the sum of the coefficients for total absorption and total scattering. 

^. =k m +k a and Hs = °m + °a (6) [Ref. 7] 

where molecular absorption coefficient 

k aerosol coefficient 

a 

a m molecular scattering coefficient 
c a aerosol scattering coefficient 

In real world conditions, there are many factors that affect the transmission of 
atmosphere. For our program, we ignored scattering and the molecular absorption 

both for simplicity and lack of real data. To get the atmospheric coefficient we ran 
Lowtran in 8-12 pm. frequency band and got extinction coefficient values which are 
changing from 0.0703 to 1.166. 

The TV sensor is affected by the same underlying rules. TV systems work in 
the visible part of the electro magnetic spectrum, as seen in Figure 5. 

5. Lowtran 

Lowtran is a FORTRAN Computer code designed to calculate atmospheric 
transmittance and radiance, averaged over a 20 cm'' intervals in steps of 5 cm' 1 in the 
spectral range of 350 to 40,000 cm' 1 . The code uses single parameter band models for 

molecular absorbtion. [Ref. 8] Figure 6 illustrate the resolution of the code near 2200 
cm' 1 wave number. This is the 4,4 - 4,6 pm wavelength range with absorbtion due to 
N 2 0 and C0 2 [Ref. 12]. 

In our program, we used Lowtran to calculate the atmospheric extinction coef¬ 
ficient for different paths. The change of path in the atmosphere changes the extinc¬ 
tion coefficient. For our program, the coefficients have been calculated for each 5° 
increments from 0° to 90°. The calculated data is kept in a different file and is read 
into the program interactively depending on the angle value coming from the top view 
angle slider. In Figure 7, the plot of the calculated atmospheric extinction coefficients 
are seen. These are the calculated values at 4 km range for each top view angle theta 
where the measured angle is the angle from zenith to horizon as seen in Figure 8. 
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6. The Detectors And The Detection Algorithms 

The detector is a sensor which measures the energy that reaches to it. The en¬ 
ergy which reaches the detector differs from the zero range value because of atmos¬ 
pheric extinction. The IR detector scans within a certain frequency range. The energy 
is collected by a lens and directed on the detector. An infrared interference filter is 
placed in front of the detector to limit the spectral region or band of the energy reach¬ 
ing the detector. The detector generates voltage which is proportional to the energy 
arriving from the target. This voltage is processed through some electronic circuits 
and then changed to target surface temperature. [Ref. 5] 

Sensors scan the field of view. The scan techniques of sensor changes the 
structures of the sensors. There are three scanning techniques: Parallel, serial, serial- 
parallel. [Ref. 7] 

There are a variety of detection algorithms for IR sensors. The most known 
ones of these algorithms are :[Ref. 9] 

♦ Contrast Box Algorithm, as shown in Figure 9 which depends on the con¬ 
trast difference of the background and the target within the gate. 

♦ Double gated filter algorithm is like a contrast box but uses a non-linear 
double gated contrast filter to localize the target. 

♦ Spoke Filter depends on scanning the image of the target and it examines 
the locally complex gradient phase angles and gradient magnitudes. 

♦ The Superslice algorithm employs multiple gray-shade thresholding and 
edge-matching to generate and segment possible target regions. In this al¬ 
gorithm, scanning is directional as shown in Figure 10. 

The detection algorithm we used is similar to the contrast box model. The pro¬ 
gram reads the red value in the screen pixels, averages them and decides on the detec¬ 
tion of the target. 

While IR detectors work with temperature, TV detectors work with luminos¬ 
ity. Luminance is often called brightness, L=dI/dA cos0. This is the luminous inten¬ 
sity per projected area normal to the line of observation. For TV systems, the 

radiation is at the visible part of the electromagnetic spectrum. Sources of light 
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enlightens the targets and targets reflect the incoming light which is picked up by the 
sensors. [Ref. 7] 

C. COMPUTER GRAPHICS 

Computer Graphics started with the display of data on hardcopy plotters and 
cathode ray tube screens. Today, it has reached an amazing point with the help of so¬ 
phisticated hardware and improved software and rendering techniques. It has grown 
to include the creation, storage and manipulation of models and images of objects. 
These models may come from diverse and expanding set of fields, and include physical, 
mathematical, engineering, architectural and natural phenomena. Computer graphics 
today is becoming largely interactive with the advent of virtual reality. [Ref. 2] 

In the early 1980s, computer graphics was a small, specialized field because of 
the expensive hardware and software. Then personal computers with built-in raster 
graphics displays popularized computer bitmap graphics. In time, different techniques 
were developed and real objects were started to be drawn by polygons on computer 
graphics screens. The greater the number of the polygons, the more realistic was the 
picture of the object. Today, 80 million polygons per picture is accepted as the thresh¬ 
old for mimicing reality. [Ref. 3] 

1. Visualizing a 3-D Target 

As we mentioned above for the visualization of targets or objects, polygons are 
used as the simplest element (component) of the drawings. These polygons have verti¬ 
ces and the information about the polygons is stored in them. Generally, these poly¬ 
gons are in triangular shape. A user who wants to create an object must calculate the 
vertices of each polygon and then store these as data or must use special software to 
create the object. All of the polygons come together and they form the model. We 
can do lighting, shading and texturing on these models. [Ref. 3] 

Every object drawn on the screen is approximated by polygons, and then these 
polygons can be rendered and with the help of lightning, texturing, and special graphics 
techniques, objects look real and 3-D. Each polygon in the object is called a facet. 
Facets have vertices. Our targets in the program are formed of polygons and the data 
about the target facets and vertex counts are presented in Table 1. 
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TARGET TYPE 

VERTEX COUNT 

FACET COUNT 

T-62 TANK 

1,384 

960 

F-16 AIRCRAFT 

533 

1,068 

DESTROYER 

628 

590 

HOUSE 

84 

28 


Table 1. Target Geometric Information 


The target models used in our work were developed by other students for pre¬ 
vious projects either by Multigen modeling software or by hand-sketching. Multigen 
has a database which contains pre-drawn target files and different models. This data¬ 
base can be used to visualize different targets, but these models can not be used di¬ 
rectly with OpenGL library functions. Therefore, we converted them to a standard 
form which can be understood by OpenGL library routines. 

2. Texturing, Lighting, Shading and Fog 
a. Texturing 

Texturing is mapping a scanned image onto the surface of an object. 
This makes the image look more realistic. The techniques we apply texturing to a poly¬ 
gon are: [Ref. 3] 

♦ The underlying color can be modulated with the texture image to get shaded 
textured models. 

♦ The textured images can be mapped onto polygons, as a decal, obscuring 
the underlying polygon color. We used this method in our program. So, 
when light and texturing are on we do not see the underlying color. When 
the lights are off we do not see the color but the texture. 

♦ A textured image can be blended with a single color. 

By specifying the s&t texture mapping coordinates, [Ref. 10] we can cover the 
polygon with the desired image. In our program, we used texturing in the TV part to 
be able to get a real reflection of the world. In this part of the program, we used dif¬ 
ferent textures to cover the environment and the objects. [Ref. 3] 
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b. Lighting 

As in nature, in computer graphics everything reflects light and color. 
For the reflection of light, there must be a light source. In nature, the light source is 
the sun. In computer graphics the light source is created by the programmer. Lighting 
is required to give the objects a three dimensional and real effect. To use a shading 
model for 3-D objects, there must be a lighting model. Ambient, diffuse and specular 
illumination models are used in developed systems. In the real world, the perception of 
the eye depends on the distribution of photon energies that arrive and trigger the cone 
cells in the eye. The photons come from a light source, and these are either reflected 
or transmitted. This forms the basic idea of lighting in OpenGL. [Ref. 10] 

c. Shading 

Shading is a result of lighting. There are two types of shading: constant 
and Gouraud. In constant shading, a single color is computed for an entire polygon, 
based on the position of the light source and the normal vector of the polygon. 
Gouraud shading implies evaluating the illumination at each pixel level. The applica¬ 
tion of the shading to the objects gives them a realistic appearance. The shading model 
is used to calculate the intensity of the light at each pixel. [Ref. 3] In our work, we 
applied diffused and specular reflection rules to the targets in the TV part. 

d Applying Fog 

Fog is a natural phenomenon and it affects the propagation of waves 
and the energy in nature. It works as an absorbent in the atmosphere. Even though, 
we do not have the absorption models in our program, we simulated this natural effect 
in our program. OpenGL provides a fog function. In our program, we determined the 
fog coefficient as a variable and, the user can control the density of the fog by using a 
slider in the TV (visible) mode. [Ref. 10] 
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III. PROBLEM 


A. DEFINITION OF THE PROBLEM 

The problem being approached here is two fold. We have a sensor detection 
problem, and a visualization task. There are two sensors with different detection 
algorithms. 

As illustrated in Figure 11, for both the IR and TV sensor codes there is a tar¬ 
get, a transmission medium, and a sensor. For our program, we implemented four tar¬ 
gets: a tank, an aircraft, a ship and a building. This implementation can be broadened 
by creating more target files and making the necessary changes both in the interface 
and in the related programming routines. 

The detection algorithms utilize emission and reflections for the ER and TV 
sensors, respectively. In each case, we model the effect of sensor field-of-view, sensi¬ 
tivity, and spectral range. Depending on their design properties and the ambient ef¬ 
fects these sensors have different detection and lock-on ranges. 

The energy radiated by the objects travels through the atmosphere toward the 
sensor. During this travel, due to the gaseous and particle material which are available 
in the atmosphere, the energy is attenuated and only a portion of it reaches to the sen¬ 
sor. Some portion of this energy is transmitted to the detector depending on the sen¬ 
sor parameters. This is the basic story of what is happening in our computer program. 

Our work is a visualization of a real physical phenomenon. In summary, our 
problem can be defined as: the integration of a user interface design, the visualization 
of IR and TV targets under real physical conditions, and creating a tool which can be 
easily modified to measure the performance of different sensors. 

B. ASSUMPTIONS 

We made a number of assumptions to make the implementation easier in the 
short time period available. These assumptions do not divert the real output of the 
data from real calculations so much, but they do affect the precision. These assump¬ 
tions can be taken as different aspects of the project and can be implemented to expand 
the project. The assumptions which are used in our work are as follows: 
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♦ The background of the target has a constant temperature. So, a fixed color 
value corresponding to this was applied to the background. 

♦ The targets were assumed to be formed of nodes connecting facets to each 
other. Each of these nodes radiates energy. 

♦ Each of the targets was assumed to be a blackbody. 

♦ The environmental conditions for the targets were assumed to be mid¬ 
latitude summer conditions. Extinction coefficient calculations were made 
according to this assumption. 

♦ The temperature range for the targets was assumed to be between 10-90 °C. 

♦ The angle of approach of the missile was assumed to increase or decrease as 
multiples of 5 degrees within the range of 0-90 degrees. 

♦ It is assumed that each target is a Lambertian surface, i.e. a surface which 
radiates in all directions. 

♦ It is assumed that the energy which the sensor sees is proportional to the 
projection area of the target on the sensor. 

♦ It is assumed that in the atmospheric conditions in which the energy travels, 
there are no aerosol. 

♦ In the program, the delta t (time) interval is assumed to be very short for the 
events happening. There is no change in the temperature of the target dur¬ 
ing the flight of the sensor, and all of the extinction calculations are done de¬ 
pending on the initial temperature of the target. Making this kind of an 
assumption is not illogical, because the flight time of a missile is very short 
and the temperature change which may occur at the target during this period 
is negligible. 

♦ It is assumed that the detection of an IR target occurs when the intensity of 
the signal (apparent temperature) exceeds the threshold level. 

♦ It is assumed that detection of a TV target occurs when the contrast range 
in a small subset of pixels exceeds the threshold level. 

♦ Lock on is accepted to occur if there are two separate detections on the 
same area on the screen. 
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IV. DESIGN 


A. SELECTION OF WORKING ENVIRONMENT 

1. Software 

The program consists of three main components. The first part is the imple¬ 
mentation of the user interface design, the second part is the sensor/detection physics 
involved, and the third part is the implementation of the (computer graphics) 
visualization. 

For the user interface part of the program, we adopted OSF/Motif since this is 
as close to a graphics standard as is available. We used OSF/Motif for the user inter¬ 
face and to control the program input and output. 

For the visualization aspect of the program, there were two tools we could 
choose IRIS Motif-GL or OpenGL. We used OpenGL because it is the most devel¬ 
oped computer graphics library at hand. This software provides an interface to the 
graphics hardware and its implementation is easy. There are various commands and 
functions in the OpenGL library. It uses primitives and different algorithms for pro¬ 
gramming. It is more developed than the IRIS Graphics Library for developing a vis¬ 
ual simulation program. 

For the sensor physics element of the computer program, we used the C++ 
programming language. The main flow of the program was designed by using the C++ 
programming techniques. The rest of the code which uses special programming func¬ 
tions, was embedded into the main C++ code. The structure we used in the program 
was a simple multi-dimensional array structure to store and process the data. 

2. Hardware 

Computer graphics utilizes special hardware. In addition to rapid graphics 
processing, we utilized texturing for the TV sensor. Texturing requires fast computing 
and runs slowly on IRIS Indigo Workstations. Texturing is implemented in hardware 
on the IRIS Reality engines and runs much more quickly. So, our hardware working 
environment is the IRIS Reality engines. These workstation are very powerful and the 
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data process rate of these machines are very high. The explanatory data about these 
machines are seen in Table 2. 


IRIS POWER SERIES 

4D/340 VGX 

Number of the Microprocessors 

4 33 Mhz IP7 

FPU 

MIPS R2010/R3010 Floating 
Point Chip 

CPU 

MIPS R2000A/R3000 Processor 
Chip 

On-Board Serial Ports 

2 per CPU board 

Data Cache Size 

64 KBytes 

Instruction Cache Size 

64 KBytes 

Secondary Cache Size 

256 KBytes 

Main Memory 

40 MBytes 


Table 2. Specifications of IRIS Reality Engine 


B. SENSOR DESIGN AND STRUCTURE 

1. Infrared Sensor Design 

The IR sensor is designed to mimic a threshold measurement technique in the 
8-12 pm. atmospheric window. Objects are modeled using temperature distribution as 
described below. Only emitted radiation is modeled so, the sensor might best be con¬ 
sidered a night - time model. The threshold level for sensors were chosen arbitrarily 
because of lack of real sensor characteristics. The unit for the threshold level is pixel 
intensity value within the range of 0 - 255. 

The main idea for the IR part of the program depends on the representation of 
the temperature with colors. For this, there is a temperature scale on the screen, and 
the colors correspond to different temperature values of the pixels. 

To implement the IR part of the program, we intended to write a program that 
simulates a generic IR sensor. We created some 3-D targets and applied color values 
to each vertex of these targets. This color on the target was interpolated from vertex 
to vertex and for each pixel on the screen between two points the temperature was ex¬ 
trapolated by the OpenGL software according to Gouraud shading as seen in Figure 
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12. This automatically provides an interpolated temperature distribution on the sur¬ 
faces of the targets. 

Since the radiated target temperature changes depending on the distance and 
angle between the target and the sensor, the color on each surface of the target 
changes as we increase or decrease the distance. The reason for this change is the at¬ 
mospheric extinction. The program reads in the initial data file. Prior to drawing the 
object, the temperature value of each vertex is converted to a power value and extinc¬ 
tion is applied to this power to calculate the power at sensor. The power at the sensor 
is converted back to temperature values. These temperature values are used as input 
parameters to figure out the color values from a look-up table. These color values are 
applied to the related vertices and the color of the whole target is changed. 

The IR model depends on the relation between the temperature and the color 
values of the computer. The detection in the IR part of the program depends on this 
reading the color value of each pixel on the screen. The detection is determined by 
the hottest point or area on the screen. The hottest point on the screen is the point 
which has the highest averaged R value, i.e. the point which has the highest 
temperature. 

The relation between the temperature and the color holds because in the color 
table presented in Chapter V, the R values of all colors were sorted in order. A color 
which has a small R value corresponds to a low temperature value. 

2. Infrared Data Structure 

In the program, we have four main arrays. One of them is used for storing the 
vertex data, one for storing the temperature of each vertex, one for storing the initial 
temperature of the vertex, and the last is used for storing the current (instantaneous) 
temperature of the vertices. Data is read into these arrays at the beginning of the pro¬ 
gram and then they are processed in these arrays. Vertex and initial temperature arrays 
are used once by the program at the beginning of the program. The current tempera¬ 
ture and color arrays are used repeatedly during the program. Also the color values of 
the ground are stored in a global, one-dimensional array and updated with the target 
calculation procedures. The data transfer among these arrays is as seen in Figure 16. 
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3. TV Sensor Design 

As we mentioned before, TV sensors generate electrical pulses which represent 
light intensities present in the optical image focused on the tube. The light intensities 
(luminance) are called gray level of the targets. The difference between the gray levels 
of a target is called contrast. 

In nature, how we recognize objects is by detecting the contrast between the 
object and its background. So, we used the same idea in our TV detection algorithm. 
This idea is similar to IR detection, too. We used fog feature to emulate the atmos¬ 
pheric extinction because we were not using false color graphics implementation used 
for the IR sensor. 

In the program, the whole screen of the computer was scanned and the contrast 
difference between pixels was checked. The area which gives the highest contrast 
difference was accepted to be the possible target. 

C. VISUALIZATION DESIGN 

The user interface in our program has two purposes. 

♦ To control the input for the program. 

♦ To test performance for different sensors whose parameters are different 
from each other. 

There are a lot of parameters that control the infrared signature of the targets. 
These parameters are related with the temperature of the background, temperature of 
the target, atmospheric propagation conditions and the properties of the sensor. These 
parameters are in the user interface. The user can enter the data and measure the per¬ 
formance depending on these parameters. Inclusion of all of the programming rou¬ 
tines that perform these calculations have not been completed in the program. 

On the interface, the basic options for the user at the beginning are to choose 
the type of the target, the type of the sensor and background. The default values are a 
tank for the target and the IR for the sensor with the background of soil. After select¬ 
ing these initial data, the display button displays the initial appearance of the target at 
the range of 4 km. We designed this user interface to give different options to the user 
and make the control of the program easier. 
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Another part of the user interface is on the display screen. There are three slid¬ 
ers on this screen. These sliders have been designed to control the input data of the 
program and the position of the program interactively and rapidly. One of these slid¬ 
ers controls the distance of the target to the sensor, the other slider controls the ap¬ 
proach angle of the missile to the target. The last slider controls the viewing angle of 
the target temperature distribution. These sliders are not only used for controlling the 
appearance of target, but also as an input source to the program. We selected sliders 
as an input and control mechanism because, they are easy to use, fast and interactive. 
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V. THE PROGRAM 


A. LOGICAL STRUCTURE 

The problem of visualizing the output of EOTDA software consists of two 
parts. These sections are the IR part and the TV part. Therefore, we separated the 
problem into two general subsections at the beginning. 

For the IR part of the program, the main hardship was first to be able to calcu¬ 
late the energy for the target model, second to represent of the temperature with the 
corresponding color value. To be able to solve these, we accepted a temperature 
range which can be expressed by the color Table 3. 


COLOR NAME 

RVALUE 

G VALUE 

B VALUE 

TEMP. RANGE 

Dark Blue 

0 

0 

1 

T<=10 

Dirty Green 

0.06 

0.31 

0 

10<T<=15 

Light Green 

0.13 

0.88 

0 

15<T<=20 

Light Blue 

0.19 

1 

1 

20<T<=25 

Brown 

0.25 

0.25 

0 

25<T<=30 

Purplish 

0.31 

0 

1 

30<T<=35 

Gray 

0.38 

0.38 

0.38 

35<t<=40 

Dark Brown 

0.44 

0.58 

0 

40<T<=45 

Cherry 

0.5 

0 

0.5 

45<T<=50 

Fading Red 

0.56 

0.16 

0.56 

50<T<=55 

Purple 

0.63 

0 

1 

55<T<=60 

Dark Gold 

0.69 

0.82 

0 

60<T<=65 

Dirty Red 

0.75 

0 

0.4 

65<T<=70 

Dirty Pink 

0.81 

0.37 

1 

70<T<=75 

Pinkish 

0.88 

0 

1 

75<T<=80 

Gold 

0.94 

0.95 

0 

80<T<=85 

Dark Red 

1 

0 

0 

85<T 


Table 3. Temperature Color Values for the IR Part of the Program 


The colors in this table were loaded into an array according to the temperature 
distribution on the target. Then the problem was to decide detection and lock-on 
precisely. 
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Since there is a one-to-one correspondence between the pixel colors on the 
screen and the temperature of the target, we can determine detection by using the R 
values of the pixels. In Table 3, it can be seen that the color values are sorted accord¬ 
ing to their increasing R values. 

After working on the general theory of the problem, we made a design decision 
depending on the flow of the events in the main program. The main flow of the pro¬ 
gram for the TV and the IR part is as seen in Figure 13 and 14. 

For the TV part of the program, what we see is close to the world as perceived 
by human eye. The detection is determined depending on the contrast difference be¬ 
tween the target and the background. So, what we did was, we used a multi¬ 
dimensional array to store the contrast values of the screen, and then by finding the 
highest contrast box, we decided that there was a detection or a lock-on. The criteria 
for detection is, if the signal value of the screen pixel is above the chosen threshold 
level, then the program gives detection. If there is detection at the same screen location 
for a second time, then it is accepted as lock-on. 
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B. COMPONENTS OF THE PROGRAM 


Table 3 contains a brief overview of the files in our program. The source code 
listings are contained in Appendix B. 


FILE NAME 

DESCRIPTION / EXPLANATION 

ExtinctionCoef. dat 

The calculated Lowtran output for atmospheric extinction 
coefficient. 

NPSimage.[C,.h] 

Image handling routines. 

*.rgb files 

Scanned images of real objects which are used for textur¬ 
ing to cover the models. 

displayjist.h 

Special OpenGL file which includes the drawing routines 
for the models. 

draw support. [C, h] 

The functions and their declarations which handle the 
drawing duty of the objects. 

eotdamain.[C,.h] 

Main flow of the program including the interface and the 
physical calculations of the program. 

eotda_j 5 lobals.il 

This header file includes the global variables and the con¬ 
stants for the whole program. 

houselist.C 

Data for drawing the building target. 

shiplist.C 

Data for drawing the ship target. 

planelist.C 

Data for drawing the aircraft target. 

materialsupport. [C, .h] 

Properties related with the lightning and the material prop¬ 
erties of the targets. 

texturesupport.C 

Functions related with texturing. 


Table 4. The main files in the program 
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C. PROGRAM INPUT AND OUTPUT 


The program takes files as input and outputs the visual color values on the 
screen and the numerical distance and angle values on the sliders. But as future work, 
a cursor connected to mouse movement can read pixel temperature values on the 
screen. So, there is no output file which is created by the computer. The input files 
which are read by the program load the arrays and these files are as listed in Table 5. 


FILE NAME 

DESCRIPTION / EXPLANATION 

AircraftlnitialFile. dat 

Vertex number, 3-D coordinates and the initial tempera¬ 
tures of the facets of the aircraft model. 

BuildingInitialFile.dat 

Vertex number, 3-D coordinates and the initial tempera¬ 
tures of the facets of the building model. 

ShipInitialFile.dat 

Vertex number, 3-D coordinates and the initial tempera¬ 
tures of the facets of the ship model. 

TankInitialFile.dat 

Vertex number, 3-D coordinates and the initial tempera¬ 
tures of the facets of the tank model. 


Table 5. Input Files and their Contents. 


The files in the table above contains 3-D information about the coordinates of 
the targets and also information about the temperature of each node on the facets. The 
format of the input files is presented in Figure 15. 

D. OPERATION 

Our modeling and implementation of the program is summarized as follows: 

The targets in the figures consist of facets. The more complex the target is, 
the more facets it has. The number of the facets and vertices belonging to each target 
has been presented in Table 1. Each facet on the target has vertices. These vertices 
are used to store information. The temperature value and the corresponding color val¬ 
ues are stored in these vertices. 
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1. Infra-Red Sensor 

♦ At the beginning of the program, we stored the temperature information of 
the targets in the vertices of the targets. 

♦ We took each vertice on the target separately and calculated the energy ra¬ 
diated by this vertex at the stored temperature. This gave the energy of the 
vertex at zero range from the target. 

♦ The energy value we found from this vertex travels in the atmosphere. Dur¬ 
ing its travel, it attenuates. This attenuation is caused by atmospheric ex¬ 
tinction. So, we multiplied this energy value with e -uR . This multiplication 
gives varying values depending on the extinction coefficient and the distance 
between the sensor and the target. The result of this multiplication gives the 
energy value at the sensor. 

♦ Sensors work by converting the energy value to temperature value. Hence, 
we converted this energy value to the corresponding temperature value. 
This corresponding temperature value is loaded into the current temperature 
array and the corresponding color is loaded into the color matrix. Then this 
color value is applied to this point for the next drawing cycle. 

♦ The same cycle is applied to all of the vertices on the target successively. 
This procedure is the same for all the IR targets. After this procedure is 
done, the next step is to scan the screen for detection and lock-on. 

At the beginning of the program, when we run the program, the program loads 
the arrays by reading data from the initial target files. The flow of data during the ini¬ 
tialization phase of the program as shown in Figure 16. 

In Figure 16, we can see what happens in the program when it is run. The 
first thing which is done in the program is to read the related target initialization data 
file. The data stored in these files are required for the initial visualization of the target. 
If the program is connected to another program like EOTDA software, and the output 
of the EOTDA software is provided in the same format, then would work as the initial 
target data. 

The target initialization files contain data according to the sorted vertex num¬ 
bers. Each of the targets in the program is formed of a lot of vertices, and data is 
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stored in these vertices. In this file, data is in the order of vertex number, x-coordinate 
of the vertex, y-coordinate of the vertex, z-coordinate of the vertex, and the tempera¬ 
ture value of the vertex. This data is in one line and the temperature in this line is the 
initial and the basic temperature of the target during the flight or cruise of the sensor. 
This temperature is the main temperature value for the rest of the extinction 
calculations. 

After the vertex coordinates are read into the vertex array (the index of the ar¬ 
ray is accepted as the corresponding vertex number) the last value which is read in is 
the temperature value of that vertex. This is the temperature value of the vertex at 
zero range. These values which are in the initial temperature array are subjected to ex¬ 
tinction calculations for each vertex, and then the newly found values are put into the 
current temperature array. 

After the calculation of the new temperatures, the corresponding color values 
are stored into the color matrix and then the draw routine in the program is invoked. 
This provides an active color temperature calculation depending on the input parame¬ 
ters. The temperature related color values are applied to the vertices and appears on 
the screen representing the temperature of each point on the target. 

The color values which are applied to the screen are the real values corre¬ 
sponding to the real estimated color values of the vertices i.e. the target. Because the 
color values have been sorted increasingly according to their R values, at low tempera¬ 
tures the corresponding color value on the screen is going to have a smaller R value, at 
high temperatures, the corresponding color value is going to have a high R value. 
This idea helps us when scanning the screen colors for detection. 

The screen scan pattern was inspired from the contrast box method which was 
invented by Texas Instruments and the scanning method which is used in the program 
is as seen in Figure 17. [Ref. 9] 

The scanning box on the screen covers an area of 200x200 pixels on the screen. 
This box moves on the screen from left to right like a scanner. At the end of each 
move, the red color values within the pixel area which is covered by the box are read, 
and the average value of the red values for the pixels in the box is found. These values 
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are stored in an array and they are compared with each other in order to be able to find 
the hottest (the reddest) area on the screen. 

The criteria for lock-on is the detection of the same point twice. Detection oc¬ 
curs when the biggest averaged red color value of a box exceeds the pre-accepted R 
value. 

2. Television Sensor 

TV sensors are designed to respond to a band of electromagnetic energy. The 
incoming radiation creates different level of illumination on the screen. This illumina¬ 
tion is called the gray level, and the relation of the gray levels of the pixels gives the 
contrast. 

Gray level for a TV system means the summation of the R, G and B values di¬ 
vided by three, i.e[(R+G+B)/3], After finding the luminance values for each pixel in 
the scanning box on the screen, we took the value of the lowest and highest luminance 
in the box. By using the relation C=GL max -GL min , we found the C (contrast) value for 
each scanner box. Then this value is stored in the multi-dimensional array for compari¬ 
son with the other box areas on the screen. 

The values in each cell of the array, as shown in Figure 18, are the contrast val¬ 
ues, x-coordinate and y-coordinate of the scanner box. We don't need the z- 
coordinate because we read the pixels of the screen. After storing all of the C (con¬ 
trast) values for each box, these values are compared with each other and the largest 
contrast value is found. The x and y screen coordinates of the largest contrast box are 
accessible to us at any time, so we can detect at which point there is possible detection 
and we can plot this area on the screen. 

Another condition for the detection to occur is that the C value must be above 
a threshold level. Having the same area as the highest contrast area for the second 
time results in lock-on. 
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E. SOME OBSERVATIONS 


In Table 6 and Table 7, some observations of the program output are pre¬ 
sented. In Table 6, the detection range for the tank target is shown. The distance val¬ 
ues in the tables depend on the sensor characteristic and the environmental conditions. 
In this table the only chancing value is the approach angle of the sensor to the target. 
Since the approach angle changes the extinction coefficient, this angle affects the de¬ 
tection range of the sensor. 


APPROACH ANGLE (Degrees) | 

SENSOR TYPE 

0 

10 

45 

90 

IR 

0.6 km 

1.2 km 

2.4 km 

2.8 km 

TV (No Fog) 

4 km 

4 km 

4 km 

4 km 

TV (50% Fog) 

2 km 

1.8 km 

2 km 

2 km 


Table 6. The Detection Ranges of IR and TV sensors for Tank target 


In Table 6 only the approach angle of the sensor changes. In Table 7, detection 
ranges for different targets for 15° approach angle and soil background are presented. 


SENSOR/TARGET 

TANK 

AIRCRAFT 

SHIP 

BUILDING 

IR 

1.6 km 

3.2 km 

1.2 km 

2.8 km 

TV (No Fog) 

4 km 

4 km 

4 km 

4 km 

TV (50% Fog) 

2.2 km 

2 km 

2 km 

2 km 


Table 7. The Detection Ranges of IR and TV Sensors for Different Targets 


F. HOW DOES THE INTERFACE WORK 

The user interface, as shown in Figure 19, provides input to the program and it 
consists of three main parts. The first part lets the user choose the target type and enter 
data about the target area, target speed and target course. From the target menu, we 
can choose from six targets. But we have implemented four of these targets in our 
program. The targets which are ready to be displayed are tank, ship, building and air¬ 
craft. The target area, direction and speed of the target areas are put on the interface 
for future use. In our program, we do not need these values. 
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The second part of the interface is the sensor part. In this part, we provide the 
sensor information to the program. The sensor type can have three different values. 
IR, TV and laser. We have implemented the first two of these sensors. For the IR 
sensor we have three different coloring scales. These scales are bi-level scale, gray- 
level scale and the cyclic scale. We designed the cyclic scale to present more accurate 
data on the screen. The user can decide the temperature of a point on the target by 
comparing the color values on the scale and on the target. Also in the sensor section 
of the interface, we have parameters which can change the detection possibility of the 
sensor. These parameters were put in the interface for the future development of the 
program. 

The third part of the interface is the environment section. In this section, we 
provide information about the environmental conditions. The first information we pro¬ 
vide is the background. We have a four background options in the program. These 
are soil, grass, cement and asphalt. The targets can be matched to one of these 
backgrounds. 

The other parts of the environment adjusts the day of the year, hour of the day 
dependent on the sun's position, and the percentage of fog in the environment. 

The display option displays the chosen target with the chosen sensor. Depend¬ 
ing on the option which is chosen from the target and the sensor parts of the interface, 
the program switches to that output screen. 

G. HOW DOES THE MAIN PART OF THE PROGRAM WORK 

1. Infrared Scene 

All of the implemented targets can be displayed in the IR screen. The IR 
screen shows the temperature distribution on the target. In the display screen, we have 
three sliders which let the user control the program and which also provide information 
to the program. Since the purpose of the program is to find the detection and the lock- 
on range depending on the sensor parameters, the user can play with the distance and 
the viewing angle of the target. Each distance and angle provides different input to 


31 



the program because these parameters affect the atmospheric extinction coefficient for 
the program. 

In Figure 20, the tank target is seen on the IR screen. On this screen, the exit 
button lets the user to get out of the program, the input button lets the user go back to 
the input interface and provide new data or chance information. The viewing distance 
slider lets the user control the distance between the sensor and the target. 

We can visualize the target in different color scales. The cyclic scale in the dis¬ 
play gives detailed information about the temperature distribution on the temperature. 
In Figures 21, 22 and 23, you can see the IR screen outputs of other targets. 

2. Television Scene 

The TV scene works as does in the IR screen. It is a reflection of the real 
world. The sliders available on the screen, function the same as it is in the IR screen. 
The TV screen outputs for the implemented targets are as shown in Figure 24, 25, 26, 
and Figure 27. 

H. PROBLEMS IN THE PROGRAMMING 

During our work, we experienced some compile errors which are machine de¬ 
pendent. In some machines we had z-buffer problem. By adding the line : 

"XtSetArg(wargs[n],GLwNdepthSize,l);n++" in the main program, we got rid 
of this problem. This line initializes the z-buffer to a depth of 1 pixel. 

We also observed that some of the scanned *.rgb files could not be read in 
some machines because of the difference in the graphics file. We saw that using dis¬ 
play lists in OpenGL programming boosts the program. So, all texturing and hypertext 
fonts are encapsulated of display lists. We observed that as the size of the program ex¬ 
pands, the program becomes more increasingly error prone. 

To write error-free programs, we applied the following useful software engi¬ 
neering techniques: 

♦ Modulation 

♦ Encapsulation 

♦ Top to bottom design 

♦ Testing by control loops 
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♦ Well documentation 

♦ Testing of each unit separately. 
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VI. SUMMARY AND CONCLUSIONS 


The work which has been presented in this thesis is a visualization of IR and 
TV detection systems. This work was started with the intention of writing a visualiza¬ 
tion part of the EOTDA software. Due to the lack of real data, the project turned out 
to be a generic program for IR and TV sensors. 

The program was written in C++ with OpenGL programming techniques and 
OpenGL graphics library functions. The interface part of the program was written in 
C++ and with OSF/Motif library function. The whole program was implemented un¬ 
der Unix operating system and the program can be transferred to other Unix systems if 
that system has the appropriate operating system and C++ and correct version of 
OpenGL libraries. 

Computer programs similar to our project have been written by private compa¬ 
nies, and their work took seven years to develop. This shows that there is a lot more 
to do on this computer program. 

The program we wrote can be used for two purposes. The first usage might be 
the visualization of EOTDA software output, which would require that the output files 
of the EOTDA software be converted to the used format by the visualization program. 
The second usage might be a generic performance tester for different sensors and en¬ 
vironment conditions. 

The user interface provided at the beginning of the program provides input to 
the program. But some of the routines related with this input have not been developed 
in the program because of time constraints. So, these routines can be developed or im¬ 
proved and the whole program can be used as a performance tester of sensors which 
have different construction parameters. 

The output file of the EOTDA file contains the facet number of each facet and 
the corresponding temperature to this facet. But the input file of our visualization pro¬ 
gram includes the vertex number, x, y, and z coordinates, and temperature of this facet. 
Since we were unable to obtain EOTDA data at the beginning of our project, we cre¬ 
ated our input file format and data for our own program. 
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The first thing that can be done on the program is that the program can be in¬ 
tegrated to EOTDA software. This requires a small change in the output files of the 
EOTDA software. 
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APPENDIX A. FIGURES AND SCREEN OUTPUTS 



Figure 1. The change of the blackbody curves with temperature "From Ref[7]" 
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CALCULATION OF PROJECTED AREA OF A TARGET 



Figure 2. The calculation of the projection area of a target 
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Figure 3. Schematical appearance of a TV camera tube "From Ref. 7" 
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Transmission (%) 


Figure 4. The Representation of Contrast Values on Computer Screen 
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Figure 5. Atmospheric Transmission Windows "From Ref 16' 
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Figure 6 . Transmittance of the Athmosphere "From Ref. 12" 



Figure 7. The plot of the atmospheric extinction at the range of 4 km. for 8-12 IR window 
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Figure 8. The measurement of theta in Lowtran Code 



Figure 9. The Schematical Representation of the Contrast Box Algorithm 









Figure 10. The Super-Slice Algorithm 


THE RADIATION MEDIUM 



Atmospheric Conditions Target 


Figure 11. The Radiation Medium 



Figure 12. Automatic application of the Gouraud Shading by hardware 



















Figure 13. Infrared Sensor Logic Flow Chart 
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Figure 14. TV sensor Logic Flow Diagram 
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Figure 19. The Screen Output Of The User Interface 
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Figure 22. The JR Screen Output Of The Aircraft Target. 
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Figure 24. The TV Screen Output Of The Tank Target. 


53 


































m 













































Figure 26. The TV Screen Output Of The Aircraft Target. 






































































Figure 27. The TV Screen Output Of The Building Target. 
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APPENDIX B. 


/************************************************************************ 

* This is eotda_main.C * 

* * 

* This program is written to form a core program for EOTDA software and * 

* utilizes the 3-D graphical target models.The program basically written * 

* by using openGL programming language .The program consist of basically * 

* three major parts.First part is related with the physical inputs and * 

* calculations, second part is User Interface Design,and finally the third * 

* part is related with displaying the graphical model based on the inputs * 

* entered by the user. Since texturing costs a high computing time of CPU * 

* this program is developed under Reality Engines. * 


Authors : Ltjg. Mustafa YELMAZ 

Ltjg. Mehmet GORGULU 




#include <iostream.h> // C++ I/O subsystem.. 

#include <stdio.h> 

#include <stdlib.h> // Get the exit() functiion definition 

#include <math.h> 

#include <fstream.h> 

#include <stream.h> 

#include <strstream.h> 

#include <iomanip.h> 

#include <Xm/MainW.h> 

#include <Xm/Xm.h> // Get the Motif library. 

#include <Xm/Form.h> 7/ We are going to use a Form container widget. 

#include <Xm/Frame.h> // Get the Frame widget. 

#include <Xm/RowColumn.h> //Get the RowColumn widget. 

#include <Xm/CascadeB.h> // Get the CascadeButton widget. 

#include <Xm/PushB.h> // Get the PushButton widget. 

#include <Xm/PushBG.h> 

#include <Xm/DialogS.h> // DialogShell widgets 

#include <Xm/Label.h> // Label widgets 

#inchrde <Xm/LabelG.h> // Label Gadets. 

#include <Xm/Separator.h> // Separator widgets 

#include <Xm/Scale.h> //Get the Scale widget. 

#include <Xm/MessageB.h> // Get the MessageBox. 

#include <Xm/Text.h> // Get the Text Widget. 

#include <Xm/TextF.h> //Get the Text Widget. 

#include <Xm/ToggleBG.h> // Get the ToggleBG Widget. 

#include <Xm/ArrowBG.h> // Get the ArrowBG Widget. 

#include <GL/GLwMDrawA.h> // We are going to use an OpenGL Motif Draw widget 
#include <Xll/StringDefs.h> 

#include <Xll/keysym.h> 

#include <GL/gl.h> // Get the OpenGL required includes. 

#include <GL/glu.h> 

#include <GL/glx.h> 
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#include "materialsupport.h" // Get material names. 

^include "materialsupportfuncs.h" // Get the material support functions. 

#include "drawsupportjuncs.h" // Get the drawing functions. 

^include "texturesupportfuncs.h" // Get the texture support functions, 

^include "TextureDisplayList.h" 

^include "eotdafuncs.h" // Get this program's function declarations. 

#include "eotda_main.h" // Get the variables and others for "main" program. 

#include "global targets.h" // Get the globals for this program, 

^include "makeRasterFont.h" // Get the font support routines. 

//#include "NPSimage.h" // Get the NPS image class definition. 

#define _COMMON_ 

^include "eotda_globals.h" 


winddata wind; // Create an object for winddata class 

timedataTime of day; // Create an object for timedata class 

static GLuint fontHandlel,fontHandle2; // Base of font display lists. 

char bufferl[10000], // Buffer for pixel red color values ..(100x100) 

bu£fer2[10000]; // Buffer for pixel luminance values...(100x100) 
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void main(int argc, char **argv) 

{ 


XmString message_string; // Compound strings. 

Arg wargs[ 15]; // Args used with XtSetArg below. 

message_string = XmStringCreateLtoR ("Please Wait Initializing the Display Lists.", charset); 

// Fallback resources for the application. 

static GLbyte * fallback_resources[] = { 

"*form*background: SGISlateBlue", // /usr/lib/rgb.txt 

"*frame*shadowType: SHADOW_IN", 

"♦glwidget* width: 1400", 

"*glwidget*height: 900", 

"*glwidget*rgba: TRUE", 

"*glwidget*doublebuffer: TRUE", 

"*glwidget*allocateBackground: TRUE", 

"eotda_main*geometiy: 1400x1000+0+0", //width, height, xoff, yoff 
"InputWindow*geometry: 1400x1000+0+0", 

"info*geometry: 10x10+800+800", 

NULL 

}; 


// Create the top level widget, 
toplevel = XtAppInitialize( 

&app_context, // Application context. 

"eotda main", // Application class. 

NULL, 0, // Command line option list. 

&argc, argv, // Command line args. 

fallbackresources, 

NULL, // Argument list. 

0); // Number of arguments. 

// Create a pop_up display window to inform the user during the initialization proccess of the 
// program. 


n = 0; 

XtSetArg (wargs[n], XmNmessageString, message_string); n++; 
info = XmCreatelnformationDialog (toplevel, 

"Please Wait Initializing the Display Lists.", wargs, n); 


XtManageChild (info); 


// Position the pop_up display in the center of the screen.. 
n = 0; 

XtSetArg (wargs[n], XmNx, 200); n++; 

XtSetArg (wargs[nj, XmNy, 200); n++; 

XtSetValues (XtParent(info), wargs, n); 


// Get rid of compound strings. 
XmStringFree(message_string); 


// Create Form widget. This is the top level container for below widgets. 
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n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 30); n++; 

MainWindow = XmCreateForm(toplevel, "MainWindow", wargs, n); 
XtManageChild(MainWindow); 

// Call the initializing function.This the first window we see when the program starts execution.. 
BuildlnputScreenDialogQ; 


// Create widget for display. But do NOT manage till "Display" button is pushed 
n = 0; 

XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNbottomPosition, 27); n++; 

DisplayWindow = XmCreateForm(MainWindow, "Display Window", wargs,n); 


// Create a frame that holds the display window 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2 ); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_IN); n++; 

MainDisplayFrame = XmCreateFrame (DisplayWindow, "MainDisplayFrame", wargs, n); 
XtManageChild (MainDisplayFrame); 


// Create Form widget. This is the container widget for rowcolumn widgets.Do NOT manage 
// this widget either.Manage it when "Display" button is pressed... 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNtopPosition, 27); n++; 

XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH FORM); n++; 
XtSetArg(wargsfn],XmNleftAttachment, XmATTACH FORM); n++; 

ControlWindow = XmCreateForm(MainWindow, "ControlWindow", wargs, n); 

// Create rowcolumn widget to display information. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargstn],XmNbottomAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg(wargs[n],XmNpacking, XmPACK_TIGHT); n++; 
re = XmCreateRowColumn(Control Window, "rowcol", wargs,n); 

XtManageChild (rc); 

// Call the function that creates the Scales under the bulletin board and manages them.. 

Create The Scales(rc); 

// Create a GL widget & initialize the z-buffer."XtSetArg(wargs[n], GLwNdepthSize,l); n++;" 

// is essential to initialize the z-buffer. Without this line z-buffering doesn't work and there will be // no hidden 
surface elemination... 
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n = 0; 

XtSetArg(wargs[n], GLwNdepthSize, 1); n++; 

glw = GLwCreateMDrawingArea(MainDisplayFrame, "glwidget", wargs, n); 
XtManageChild (glw); 

// Add callbacks to the glw widget. 

XtAddCallback(glw, GLwNginitCallback, initCB, (XtPointer) NULL); 
XtAddCallback(glw, GLwNexposeCallback, exposeCB, (XtPointer) NULL); 
XtAddCallback(glw, GLwNresizeCallback, resizeCB, (XtPointer) NULL); 
XtAddCallback(glw, GLwNinputCallback, inputCB, (XtPointer) NULL); 


// Add in the work procedure. A work procedure is called repeatedly whenever there are 
// events to process. 

workprocid = XtAppAddWorkProc(app_context, (XtWorkProc)drawWP, 
(XtPointer)NULL); 

// Instantiate it now the toplevel container widget. 

XtRealizeWidget(toplevel); 

// Loop for events. 

XtAppMainLoop(app_context); 


} 
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/************************************************************************* 


* Here we create all the widgets for input screen and manage the parent * 

* InputWindow widget so that we will see this screen first when the * 

* program is executed. When the "Display" button is pressed on this screen * 

* the InputWindow widget is unmanaged and DisplayWindow widget is managed * 

* so we will see the graphical display. * 

**************************************************************************/ 


void BuildlnputScreenDialogO 

{ 


Widge rc , //Locally defined container widget, 

rcl, // All other child widgets 

rc2, 
rc3 , 
rc4 ; 

Arg wargs[ 15]; // Args used with XtSetArg below. 

// Here is the global form that holds four different forms under it. 
n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 40); n++; 

InputWindow = XmCreateForm(MainWindow, "InputWindow", wargs, n); 
XtManageChild(InputWindow); 

// This the form that holds target related widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNfiractionBase, 20); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNrightPosition, 20); n++; 

XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNbottomPosition, 20); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 

Widget TargetForm = XmCreateForm(InputWindow, "TargetForm", wargs, n); 
XtManageChild(TargetForm); 

// Create the target related widgets here.. 
create_target_screen_widgets (TargetForm); 

// This the form that holds sensor related widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 20); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNleftPosition, 22); n++; 

XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNbottomPosition, 20); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 

Widget SensorForm = XmCreateForm(InputWindow, "SensorForm", wargs, n); 
XtManageChild(SensorForm); 

// Create the sensor related widgets here.. 
create_sensor_screen_widgets (SensorForm); 
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// This the form that holds Background related widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 40); n++; 

XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNtopPosition, 20); n++; 

Widget BackgroundForm = XmCreateForm(InputWindow, "BackgroundForm", wargs, n); 
XtManageChild(BackgroundForm); 

// Create the background related widgets here.. 
create_background_screen_widgets (BackgroundForm); 


}// End of BuildlnputScreenDialogO- 




y************** ******* ** ************************************** *********** 

* This function creates the necessary widgets for the target part * 

* of the input_screen. * 

*************************************************************************/ 
void create_target_screen_widgets(Widget parent) 

{ 

Widget pulldown, // Pulldown widget. 

option; // Cascade for the pulldown. 


XmString label string; // String. 

Arg wargs[10]; //Same old args stuff. 

// This the rowcolumn that holds Target puldown widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNleflAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget TargetRowColumnl = XmCreateRowColumn(parent, 

"TargetRowColumnl", wargs, n); 
XtManageChild(TargetRowColumn 1); 


// Add a pulldown to the RowColumn. 
n = 0; 

pulldown = XmCreatePulldownMenu (TargetRowColumnl, "pulldown", wargs, n); 


// Add the entries to the pulldown. 

make_pulldown_entry(pulldown, "TANK ", targetCB, (XtPointer) & TANK); 
make_pulldown_entry(pulldown, "TRUCK ", targetCB, (XtPointer) & TRUCK); 
make_pulldown_entry(pulldown, "SHIP ", targetCB, (XtPointer) & SHIP); 
make_pulldown_entry(pulldown, "HELO ", targetCB, (XtPointer) &HELO); 
make_pulldown_entry(pulldown, "BUILDING", targetCB, (XtPointer) & HOUSE); 
make_pulldown_entry(pulldown, "AIRCRAFT", targetCB, (XtPointer) & PLANE); 

// Create the option menu that has the pulldown under it. 
label string = XmStringCreateSimple("TARGET : "); 

n = 0; 

XtSetArg(wargs[n], XmNlabelString, label string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldown); n++; 

option = XmCreateOptionMenu (TargetRowColumnl, "optionMenu", wargs, n); 

// Manage the label widget. 

XtManageChild (option); 

// Free the string.. 

XmStringFree(label_string); 
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// This the rowcolumn that holds Material puldown widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetRowColumnl); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget TargetRowColumnll = XmCreateRowColumn(parent, 

"TargetRowColumnll", wargs, n); 
XtManageChild(TargetRowColumn 11); 


// Add a pulldown to the RowColumn. 
n = 0; 

XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 

Widget puIldown2 = XmCreatePulldownMenu (TargetRowColumnll, ''pulldown2", wargs, n); 

// Add the entries to the pulldown. When material is included into the computations 
// material names are to defined in global targets.h file and used here as user data 
// pased into the callback. 

make_pulldown_entry(pulldown2, "STEEL ", materialCB, (XtPointer) &TANK); 
make_pulldown_entiy(pulldown2, "ALIMINIUM", materialCB, (XtPointer) & TANK); 
makejulldown_entry(pulldown2, "CONCERETE", materialCB, (XtPointer) & TANK); 
make_pulldown_entry(pulldown2, "COMPOSITE", materialCB, (XtPointer) & TANK); 

// Create the option menu that has the pulldown under it. 
label_string = XmStringCreateSimple("MATERIAL: "); 
n = 0; 

XtSetArg(wargs[n], XmNlabelString, label_string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldown2); n++; 

Widget option2 = XmCreateOptionMenu (TargetRowColumnll, "optionMenu2", wargs, n); 
XtManageChild (option2); 

// Free the string.. 

XmStringFree(label_string); 

// This the rowcolumn that holds TargetArea widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetRowColumnll); n++; 
XtSetArg(wargs[n],XmNleflAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK_TIGHT); n++; 

Widget TargetRowColumn2 = XmCreateRowColumn(parent, 

"TargetRowColumn2", wargs, n); 
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// This is the Target area gadget. 

XtVaCreateManagedWidget("Target Area (A): ",xmLabelGadgetClass, 

TargetRowColumn2,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Target Area. 

Widget TextWidgetl = Xt VaCreateManagedWidget("Text Widget 1", 
xmTextWidgetClass ,TargetRowColumn2 ,NULL); 

// Add Callbacks to the TextWidgetl widget. 
XtAddCallback(TextWidgetl,XmNmodilyVerrfyCallback,TargetAreaCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(TargetRowColumn2); 


// This the rowcolumn that holds TargetHeading widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetRowColumn2); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACKTIGHT); n++; 
TargetRowColumn3 = XmCreateRowColumn(parent, 

"TargetRowColumn3", wargs, n); 


// This is the Target heading gadget. 

XtVaCreateManagedWidget("Target Heading : ",xmLabelGadgetClass, 

TargetRowColumn3,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the TargetHeading.... 

Widget TextWidget2 = XtVaCreateManagedWidget("TextWidget2”, 
xmTextWidgetClass ,TargetRowColumn3 ,NULL); 

// Add Callbacks to the TextWidget2 widget. 
XtAddCallback(TextWidget2,XmNmodifyVerifyCallback,TargetHeadingCB,NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(TargetRowColiunn3); 

// This the rowcolumn that holds Target Speed widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetRowColumn3); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACK_TIGHT); n++; 

TargetRowColumn4 = XmCreateRowColumn(parent, 

"TargetRowColumn4", wargs, n); 
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// This is the Target speed gadget. 

XtVaCreateManagedWidget("Target Speed (km/hr) :",xmLabelGadgetClass, 

TargetRowColumn4,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Target Speed. 

Widget TextWidget3 = XtVaCreateManagedWidget("TextWidget3", 
xmTextWidgetClass ,TargetRowColumn4 ,NULL); 

// Add Callbacks to the TextWidget3 widget. 
XtAddCallback(TextWidget3,XmNmodifyVerifyCallback,TargetSpeedCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(TargetRowColumn4); 


// Create the bulletin board widget 
n=0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetRowColumn4 ); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
TargetBB5= XmCreateBulletinBoard(parent, "TargetBB5", wargs, n); 
XtManageChild(TargetBB5); 

// Create a Frame widget to make it so we can resize the window, 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2 ); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_OUT); n++; 
Widget TargetFrame5 = XmCreateFrame (TargetBB5, "TargetFrame5", 

wargs, n); 

XtManageChild (TargetFrame5); 


// This the rowcolumn that holds Target Status widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 
XtSetArg (wargs[n], XmNradioBehavior, TRUE); n++; 

Widget TargetRowColumn5 = XmCreateRowColumn(TargetFrame5, 

"TargetRowColumn5", wargs, n); 


// Create the label for Engine Status, 
n = 0; 

Widget Label = XtVaCreateManagedWidget("Engine Status : ",xmLabelGadgetClass, 
TargetRowColumn5,NULL); 


// Here we create a Radio Button Gadget for Engine Status. 

Widget EngineOn = XtVaCreateManagedWidget("EngineOn", 

xmToggleButtonGadgetClass ,TargetRowColumn5 ,NULL); 

// Add Callbacks to the ON gadget. 

XtAddCallback(EngineOn,XmNvalueChangedCallback,EngineStatusCB,(XtPointer) 1); 
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// Here we create a Radio Button Gadget fot Targut Status. 

Widget EngineOff = XtVaCreateManagedWidget("Engine01F’, 

xmToggleButtonGadgetClass ,TargetRowColumn5 ,NULL); 


// Add Callbacks to the OFF gadget. 

XtAddCallback(EngineOff,XmNvalueChangedCallback,EngineStatusCB,(XtPointer) 0); 

// Now manage the parent RowColumn Widget. 

XtManageChild(TargetRowColumn5); 


// Create the bulletin board widget 

n=0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetBB5); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
TargetBB6= XmCreateBulletinBoard(parent, "TargetBB6", wargs, n); 
XtManageChild(TargetBB6); 

// Create a Frame widget 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2 ); n++; 

XtSetArg (wargsfn], XmNshadowType,XmSHADOW_OUT); n++; 
Widget TargetFrame6 = XmCreateFrame (TargetBB6, "TargetFrame6", 

wargs, n); 

XtManageChild (TargetFrame6); 

// This the rowcolumn that holds Target Status widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargsfn], XmNradioBehavior, TRUE); n++; 

Widget TargetRowColumn6 = XmCreateRowColumn(TargetFrame6, 

"TargetRowColumn6", wargs, n); 


// Create the label for Fire Status, 
n = 0; 

Widget Label2 = XtVaCreateManagedWidget("Fire Status : ",xmLabelGadgetClass, 
TargetRowColumn6,NULL); 


// Here we create a Radio Button Gadget fot Fire Status. 

Widget Fired = XtVaCreateManagedWidget("Fired", 

xmToggleButtonGadgetClass ,TargetRowColumn6 ,NULL); 

// Add Callbacks to the ON gadget. 

XtAddCallback(Fired,XmNvalueChangedCallback,FireStatusCB,(XtPointer) 1); 

// Here we create a Radio Button Gadget fot Targut Status. 
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Widget NotFired = XtVaCreateManagedWidget("NotFired", 
xmToggleButtonGadgetClass ,TargetRowColumn6 ,NULL); 


// Add Callbacks to the OFF gadget. 

XtAddCallback(NotFired,XmNvalueChangedCallback,FireStatusCB ; (XtPointer)0); 


// Now manage the parent RowColumn Widget. 
XtManageChild(TargetRowColumn6); 


// Create the bulletin board widget 
n=0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,TargetBB6); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
TargetBB7= XmCreateBulletinBoard(parent, "TargetBB7", wargs, n); 

// Don't manage it here. Manage with respect to Target selected 
// XtManageChild(TargetBB7); 


// Create a Frame widget 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2 ); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_OUT); n++; 
Widget TargetFrame7 = XmCreateFrame (TargetBB7, "TargetFrame7", 

wargs, n); 

XtManageChild (TargetFrame7); 


// This the rowcolumn that holds Target Status widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNradioBehavior, TRUE); n++; 

Widget TargetRowColumn7 = XmCreateRowColumn(TargetFrame7, 

"TargetRowColumn7", wargs, n); 


// Create the label for Fire Status, 
n = 0; 

Widget Label3 = XtVaCreateManagedWidget("Heat Isolated: ", 
xmLabelGadgetClass,TargetRowColumn7,NULL); 

// Here we create a Radio Button Gadget for Isolation Status. 

Widget Isolated = XtVaCreateManagedWidgetfYes", 

xmToggleButtonGadgetClass ,TargetRowColumn7 ,NULL); 

// Add Callbacks to the Isolated gadget. 
XtAddCallback(Isolated,XmNvalueChangedCallback,IsolationStatusCB, 
(XtPointer)l); 
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// Here we create a Radio Button Gadget fot Targut Status. 

Widget Notlsolated = XtVaCreateManagedWidget("No ", 

xmToggleButtonGadgetClass ,TargetRowColumn7 ,NULL); 


// Add Callbacks to the Notlsolated gadget. 
XtAddCallback(NotIsolated,XmNvalueChangedCallback,IsolationStatusCB, 
(XtPointer)0); 


// Now manage the parent RowColumn Widget. 
XtManageChild(TargetRowColumn7); 


} // End of create_target_screen_widgets 
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/************************************************************************ 

* This function creates the necessary widgets for the sensor part * 

* of the input_screen. * 

***************************************************************** ********/ 


void createsensorscreenwidgets (Widget parent) 

{ 


Arg wargs[ 10]; // Same old args stuff. 


// This the rowcolumn that holds Sensor puldown widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 

// XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget SensorRowColumnl = XmCreateRowColumn(parent, 

"SensorRowColumnl", wargs, n); 
XtManageChild(SensorRowColumn 1); 

// Add a pulldown to the RowColumn. 
n = 0; 

Widget pulldown = XmCreatePulldownMenu (SensorRowColumnl, "pulldown", wargs, n); 
// Add the entries to the pulldown. 

make_pulldown_entry(pulldown, "IR ", SensorCB, (XtPointer) &_IR); 
make_pulldown_entry(pulldown, "TV ", SensorCB, (XtPointer) &_TV); 
make_pulldown_entry(pulldown, "LASER ", SensorCB, (XtPointer) & LASER); 


// Create the option menu that has the pulldown under it. 
n = 0; 

XmString label_string = XmStringCreateSimple("SENSOR:"); 

XtSetArg(wargs[n], XmNlabelString, label_string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldown); n++; 

Widget option = XmCreateOptionMenu (SensorRowColumnl, "optionMenu3", wargs, n); 
XtManageChild (option); 

XmStringFree(label_string); 

// This the rowcolumn that holds Sensor puldown widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH WIDGET); n++; 

// XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg(wargs[n],XmNleftWidget,SensorRowColumnl); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK_TIGHT); n++; 

SensorRCll = XmCreateRowColumn(parent, 

"SensorRCll", wargs, n); 

XtManageChild(SensorRCl 1); 
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// Add a pulldown to the RowColumn. 
n = 0; 

Widget pulldownl = XmCreatePulldownMenu (SensorRCl 1, "pulldownl", wargs, n); 

// Add the entries to the pulldownl. 

make_pulldown_entry(pulldownl, "TEMP ", IRScaleCB, (XtPointer) 

&JIEMPERATURESCALE); 

make_pulldown_entry(pulldownl, "GRAY ", IR ScaleCB, (XtPointer) &GRAYSCALE); 
make_pulldown_entry(pulldownl, "SYCL.", IR ScaleCB, (XtPointer) & SYCLIC SCALE); 


// Create the option menu that has the pulldown under it. 
n = 0; 

label_string = XmStringCreateSimple("IR_SCALE :"); 

XtSetArg(wargs[n], XmNlabelString, label_string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldownl); n++; 

Widget optionMenul = XmCreateOptionMenu (SensorRCl 1, "optionMenul", wargs, n); 
XtManageChild (optionMenul); 

XmStringFree(label_string); 


// This the rowcolumn that holds Relative Aperture widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHJWIDGET); n++; 
XtSetArg(wargs[n] ,XmNtopWidget,SensorRowColumn 1); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget SensorRowColumn2 = XmCreateRowColumn(parent, 

"SensorRowColumn2", wargs, n); 


// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Relative Aperture :",xmLabelGadgetClass, 

SensorRowColumn2,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Relative Aperture. 

Widget TextWidgetl = XtVaCreateManagedWidget("TextWidgetl", 
xmTextWidgetClass ,SensorRowColumn2 ,NULL); 

// Add Callbacks to the TextWidgetl widget. 
XtAddCallback(TextWidgetl,XmNmodifyVerifyCallback,RelativeApertureCB,NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(SensorRowColumn2); 


// This the rowcolumn that holds Aperture Dram widget under it. 
n = 0; 

XtSetArg(wargs[n].XmNtopAttachment, XmATTACH_WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn2); n++; 
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XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACKTIGHT); n++; 

Widget SensorRowColumn3 = XmCreateRowColumn(parent, 

"SensorRowColumn3", wargs, n); 

// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget('ApertureDram :",xmLabelGadgetClass, 

SensorRowColumn3,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Relative Aperture. 

Widget TextWidget2 = XtVaCreateManagedWidget("TextWidget2", 
xmTextWidgetClass ,SensorRowColumn3 ,NULL); 

// Add Callbacks to the TextWidget2 widget. 
XtAddCallback(TextWidget2,XmNmodifyVerifyCallback,ApertureDramCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn3); 


// This the rowcolumn that holds Magnification widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn3); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACK TIGHT); n++; 

Widget SensorRowColumn4 = XmCreateRowColumn(parent, 

"SensorRowColumn4", wargs, n); 

// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Magnification :",xmLabelGadgetClass, 

SensorRowColumn4,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Magnification. 

Widget TextWidget3 = XtVaCreateManagedWidget("TextWidget3", 
xmTextWidgetClass ,SensorRowColumn4 ,NULL); 

// Add Callbacks to the TextWidget3 widget. 
XtAddCallback(TextWidget3,XmNmodifyVerifyCallback,MagnificationCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn4); 


// This the rowcolumn that holds Optical Transmittance widget under it. 
n = 0; 

XtSetArg(wargsfn],XmNtopAttachment, XmATTACHJWIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn4); n++; 
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XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACKTIGHT); n++; 

Widget SensorRowColumn5 = XmCreateRowColumn(parent, 

"SensorRowColumn5", wargs, n); 


// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Optical Transmittance :",xmLabelGadgetClass, 

SensorRowColumn5,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Optical Transmittance. 

Widget TextWidget4 = XtVaCreateManagedWidget("TextWidget4", 
xmTextWidgetClass ,SensorRowColumn5 ,NULL); 

// Add Callbacks to the TextWidget4 widget. 
XtAddCallback(TextWidget4,XmNmodifyVerifyCallback, 

OpticalT ransmittanceCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn5); 


// This the rowcolumn that holds Eln. Band Width widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn5); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget SensorRowColumn6 = XmCreateRowColumn(parent, 

"SensorRowColumn6", wargs, n); 

// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Eln. Band Width :",xmLabelGadgetClass, 

SensorRowColumn6,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Eln. Band Width. 

Widget TextWidget5 = XtVaCreateManagedWidget("TextWidget5", 
xmTextWidgetClass ,SensorRowColumn6 ,NULL); 

// Add Callbacks to the TextWidget5 widget. 
XtAddCallback(TextWidget5,XmNmodifyVerifyCallback, 
ElectronicBandwidthCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn6); 
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// This the rowcolumn that holds Sensor Height widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn6); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACKTIGHT); n++; 

Widget SensorRowColumn7 = XmCreateRowColumn(parent, 

"SensorRowColumn7", wargs, n); 

// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Sensor Height :",xmLabelGadgetClass, 

SensorRowColumn7,XmNorientation, XmHORIZONTAL,NULL); 
// Create the text widget for the Relative Aperture. 

Widget TextWidget6 = XtVaCreateManagedWidget("TextWidget6", 
xmTextWidgetClass ,SensorRowColumn7 ,NULL); 

// Add Callbacks to the TextWidget6 widget. 

XtAddCallback(TextWidget6,XmNmodifyVerify Callback,SensorHeightCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn7); 


// This the rowcolumn that holds Detectivity widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn7); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget SensorRowColumn8 = XmCreateRowColumn(parent, 

"SensorRowColumn8", wargs, n); 


// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Detectivity :",xmLabelGadgetClass, 

SensorRowColumn8,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Detectivity. 

Widget TextWidget7 = XtVaCreateManagedWidget("TextWidget7", 
xmTextWidgetClass ,SensorRowColumn8 ,NULL); 

// Add Callbacks to the TextWidget7 widget. 
XtAddCallback(TextWidget7,XmNmodifyVerifyCallback,DetectivityCB,NULL); 

// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumn8); 
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// This the rowcolumn that holds Detector Elements widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 

XtSet Arg(wargs[n] ,XmNtopWidget, SensorRowColumn8); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACKTIGHT); n++; 

Widget SensorRowColumn9 = XmCreateRowColumn(parent, 

"SensorRowColumn9", wargs, n); 

// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Detector Elements :",xmLabelGadgetClass, 

SensorRowColumn9,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Detector Elements. 

Widget TextWidget8 = XtVaCreateManagedWidget("TextWidget8", 
xmTextWidgetClass ,SensorRowColumn9 ,NULL); 

// Add Callbacks to the TextWidget8 widget. 
XtAddCallback(TextWidget8,XmNmodiiyVerifyCallback,DetectorElementsCB,NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(SensorRowColumn9); 


// This the rowcolumn that holds Firing Angle widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 
XtSetArg(wargs[n],XmNtopWidget,SensorRowColumn9); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK_TIGHT); n++; 

Widget SensorRowColumnlO = XmCreateRowColumn(parent, 

"SensorRowColumnlO", wargs, n); 


// This is the Relative Aperture gadget. 

XtVaCreateManagedWidget("Firing Angle :",xmLabelGadgetClass, 

SensorRowColumnlO,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Firing Angle. 

FiringAngleTextWidget = XtVaCreateManagedWidget("FiringAngleTextWidget", 

xmTextWidgetClass, SensorRowColumnlO ,XmNvalue,"10.0",NULL); 

// Add Callbacks to the FiringAngleTextWidget widget. 
XtAddCallback(FiringAngleTextWidget,XmNmodifyVerifyCallback,FiringAngleCB,NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(SensorRowColumnlO); 

// This the rowcolumn that holds Relative Aperture widget under it. 
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n = 0; 

XtSetArg(wargs[n],XmNtopAttachment XmATTACHWEDGET); n++; 
XtSetArg(wargs[n] ,XmNtopWidget, SensorRowColumn 10); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACKTIGHT); n++; 

Widget SensorRowColumn 11 = XmCreateRowColumn(parent, 

"SensorRowColumnll", wargs, n); 


// This is the Firing Distance gadget. 

XtVaCreateManagedWidget("Firing Distance :",xmLabelGadgetClass, 

SensorRowColumn 1 l,XmNorientation, XmHORIZONTAL,NULL); 

// Create the text widget for the Firing Distance. 

FiringDistanceTextWidget = XtVaCreateManagedWidget("FiringDistanceTextWidget", 
xmTextWidgetClass ,SensorRowColumn 1 l,XmNvalue,"4000.0",NULL); 

// Add Callbacks to the FiringDistanceTextWidget widget. 
XtAddCallback(FiringDistanceTextWidget,XmNmodifyVerifyCallback,FiringDistanceCB, 

NULL) 


// Now manage the parent RowColumn Widget. 
XtManageChild(SensorRowColumnl 1); 


} 
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/************************************************************************ 

* This function creates the necessary widgets for the background * 

* part of the input_screen. * 


*********1^**************************************************************/ 


void createbackgroundscreenwidgets (Widget parent) 

{ 


Widget rowcolumn, 

button ; 

XmString labelstring; 

Arg wargs [10]; // Same old args stuff. 


// This the form that holds Background puldown widgets under it. 
n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 40); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH_POSITION); n++; 
XtSetArg(wargs[n],XmMeftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNrightPosition, 15); n++; 
XtSetArg(wargs[n],XmNbottomPosition, 30); n++; 

Widget Backgroundforml = XmCreateForm(parent, 

"Backgroundforml", wargs, n); 
XtManageChild(Backgroundform 1); 


// This the rowcolumn that holds Pull down Menu widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_POSrriON); n++; 
XtSetArg(wargs[n],XmNleflAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNtopPosition, 20); n++; 

XtSetArg(wargs[n] ,XmNleftPosition, 5); n++; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACK_TIGHT); n++; 

Widget BackgroundRowColumnl = XmCreateRowColumn(Backgroundforml, 

"BackgroimdRowColumnl", wargs, n); 

// Add another pulldown to the RowColumn (for object selection), 
n = 0; 

Widget pulldown = XmCreatePulldownMenu (BackgroundRowColumnl, 

"pulldown2", wargs, n); 

// Add the entries to the pulldown. 

make_pulldown_entry(pulldown, "SOIL ", Background 1CB, (XtPointer) & SOIL); 
make_pulldown_entiy(pulldown, "GRASS ", Background2CB, (XtPointer) & GRASS); 
make_pulldown_entry(pulldown, "SEA ", Background3CB, (XtPointer) & SEA); 
make_pulldown_entry(pulldown, "ASPHALT", Background4CB, (XtPointer) & ASPHALT); 

// Create the option menu that has the pulldown under it. 
n = 0; 
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String string 1 = "BACKGROUND:"; 
label_string = XmStringCreate(stringl); 

XtSetArg(wargs[n], XmNlabelString, label string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldown); n++; 

Widget option = XmCreateOptionMenu (BackgroundRowColumnl, "optionMenu2", 

wargs, n); 


XtManageChild (option); 


// Free the label string.... 

XmStringFree(label_string); 

// Here manage the widget. 

XtManageChild (BackgroundRowColumnl); 

// This the form that holds Display widget under it. 
n = 0; 

XtSetArg(wargs[n],XmNfractionBase, 16); n++; 

XtSetArg(wargs[n],XmNbottomAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNtopAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNtopPosition, 30); n++; 

Widget Backgroundform3 = XmCreateForm(parent, 

"Backgroundform3", wargs, n); 
XtManageChild(Backgroundform3); 

// Here define a rowcolumn widget for display, 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_POSITION); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNtopPosition, 4); n++; 
XtSetArg(wargs[n],XmNleftPosition, 8); n++; 
XtSetArg(wargs[n],XmNrightPosition, 10); n++; 
XtSetArg(wargs[n],XmNbottomPosition, 12); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK_TIGHT); n++; 

Widget BackgroundRC31 = XmCreateRowColumn(Backgroundform3, 

"BackgroundRC31", wargs, n); 


// Add an Display PushButton. 
label_string = XmStringCreateLtoR("DISPLAY", charset); 
n = 0; 

XtSetArg(wargs[n], XmNalignment, XmALIGNMENTCENTER); n++; 
XtSetArg(wargs[n], XmNlabelString, label_string); n++; 

DisplayButton = XmCreatePushButton(BackgroundRC31, "Display", wargs, n); 

// Initial callback is finished so we can manage Display button so the user can use. 
XtManageChild(DisplayButton); 


XmStringFree(label_string); 
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// Now add an Display callback for the PushButton. 

XtAddCallback(DisplayButton, XmNarmCallback, DisplayCB, (XtPointer) NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(BackgroundRC31); 


// This the form for Background widgets., 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACHPOSITION); n++; 
XtSetArg(wargs[n],XmNleftAttachment, XmATTACH POSITION); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNleftPosition, 15); n++; 
XtSetArglwargsfnJ.XmNbottomPosition, 30); n++; 

Widget Backgroundform2 = XmCreateForm(parent, 

"Backgroundform2", wargs, n); 
XtManageChild(Backgroundform2); 


// This the rowcolumn that holds Background widgets under it. 

n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_FORM); n++; 

XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 

XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargsfn], XmNpacking, XmPACKTIGHT); n++; 

Widget BackgroundRC3 = XmCreateRowColumn(Backgroundform2, 

"BackgroundRC3", wargs, n); 

// This is the Wind Speed label gadget. 

XtVaCreateManagedWidget("Wind Speed :",xmLabelGadgetClass, 

BackgroundRC3,XmNorientation, XmHORIZONTAL,NULL); 

// This is the label widget which holds the wind speed and it is modified when the wind speed is 

// changed. 

WindSpeedLabel =XtVaCreateManagedWidget("WindSpeedLabel", 
xmLabelGadgetClass,BackgroundRC3, 
XtVaTypedArg,XmNlabelString,XmRString,"0 ",3, 

XmNuserData, &wind, 

NULL); 

// Here we create an arrow to increment the wind velocity.. 

Widget UpArrow = XtVaCreateManagedWidgetf'UpArrow", 
xmArrowButtonGadgetClass,BackgroundRC3, 

XmNarrowDirection, XmARROW UP, NULL); 

// Add the callback for Up Arrow 

XtAddCallback (UpArrow, XmNarmCallback, ChangeWindSpeedCB, (XtPointer) 1); 
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// Here we create an arrow to decrement the wind velocity.. 

Widget DownArrow = XtVaCreateManagedWidget("DownArrow", 
xmArrowButtonGadgetClass,BackgroundRC3, 

XmNarrowDirection, XmARROW_DOWN, NULL); 

// Add the callback for Down Arrow 

XtAddCallback (DownArrow, XmNarmCallback, ChangeWindSpeedCB, (XtPointer) -1); 

II Here manage the widget. 

XtManageChild (BackgroundRC3); 

// Here define a rowcolumn widget for wind direction. 

n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH WIDGET); n++; 

XtSetArg(wargs[n],XmNleftAttachment, XmATTACHFORM); n++; 

XtSetArg(wargs[n],XmNrightAttachment, XmATTACH FORM); n++; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg(wargs[n],XmNtopWidget, BackgroundRC3); n++; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg(wargs[n],XmNentryVerticalAlignment, XmALIGNMENT_CONTENTS_TOP);n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK TIGHT); n++; 

Widget BackgroundRC4 = XmCreateRowColumn(Backgroundform2, 

"BackgroundRC4", wargs, n); 

// This is the Wind Direction gadget. 

XtVaCreateManagedWidget("Wind Dir :",xmLabelGadgetClass, 

BackgroundRC4,XmNorientation, XmHORIZONTAL, 

NULL); 

// Create the text widget for the WindDirection. 

Widget WindDirection = XtVaCreateManagedWidget("WindDirection", 
xmTextWidgetClass ,BackgroundRC4 ,NULL); 

// Add Callbacks to the WindDirection widget. 

XtAddCallback(WindDirection,XmNmodifyVerifyCallback,WindDirectionCB, 

NULL); 

// This is the Aeresol gadget. 

XtVaCreateManagedWidget(" Aeresol :",xmLabelGadgetClass, 

BackgroundRC4,XmNorientation, XmHORIZONTAL, 

NULL); 

// Create the text widget for the Aeresol. 

Widget Aeresol = XtVaCreateManagedWidget("Aeresol", 

xmTextWidgetClass ,BackgroundRC4 ,NULL); 

// Add Callbacks to the Aeresol widget. 

XtAddCallback( Aeresol, XmNmodifyVerifyCallback,AeresolCB, 

NULL); 

// Now manage the parent RowColumn Widget. 

XtManageChild(BackgroundRC4); 


// Here define a rowcolumn widget for wind direction. 



n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n],XmNlefLAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg(wargs[n],XmNtopWidget, BackgroundRC4); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg(wargs[n] ,XmNentry Vertical Alignment, XmALIGNMENT_CONTENTS_TOP);n++; 
Widget BackgroundRC41 = XmCreateRowColumn(Backgroundform2, 

"BackgroundRC41", wargs, n); 


// This is the Lat/Long gadget. 

XtVaCreateManagedWidget("Latitude :",xmLabelGadgetClass, 

BackgroundRC41,XmNorientation, XmHORIZONTAL, 
NULL); 

// Create the text widget for the Latitude. 

Widget Latitude = XtVaCreateManagedWidget("Latitude", 

xmTextWidgetClass ,BackgroundRC41 ,NULL); 

// Add Callbacks to the Latitude widget. 
XtAddCallback(Latitude,XmNmodifyVerifyCallback,LatitudeCB, 
NULL); 


// This is the Lat/Long gadget. 

XtVaCreateManagedWidget(" Longitude: ",xmLabelGadgetClass, 

BackgroundRC41,XmNorientation, XmHORIZONTAL, 
NULL); 

// Create the text widget for the Longitude. 

Widget Longitude = XtVaCreateManagedWidgetC'Longitude", 
xmTextWidgetClass ,BackgroundRC41 ,NULL); 

// Add Callbacks to the Latitude widget. 
XtAddCallback(Longitude,XmNmodifyVerifyCallback,LongitudeCB, 
NULL); 


// Here manage the widget. 
XtManageChild (BackgroundRC41); 


// Create the bulletin board 1 widget 
n=0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_WIDGET); n++; 
XtSetArg(wargs[n],XmNleflAttachment, XmATTACH FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACHFORM); n++; 
XtSetArg(wargsln],XmNtopWidget, BackgroundRC41); n++; 

Widget BB= XmCreateBulletinBoard(Backgroundform2, "BB", wargs, n); 
XtManageChild(BB); 
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// Create a Frame widget. 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2); n++; 

XtSetArg (wargsfn], XmNshadowType,XmSHADOW_OUT); n++; 

Widget BackgroundRC6frame = XmCreateFrame (BB, "BackgroundRC6frame", wargs, n); 
XtManageChild (BackgroundRC6frame); 


// This the rowcolumn that holds Background widgets under it. 

n = 0; 

XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACK_TIGHT); n++; 

Widget BackgroundRC6 = XmCreateRowColumn(BackgroundRC6frame, 

"BackgroundRC6”, wargs, n); 

// This is the Time Of The Day label gadget. 

XtVaCreateManagedWidget("Time Of The Day :",xmLabelGadgetClass, 

BackgroundRC6,XmNorientation, XmHORIZONTAL,NULL); 

// This is the Hour label gadget. 

XtVaCreateManagedWidget("Hour :",xmLabelGadgetClass, 

BackgroundRC6,XmNorientation, XmHORIZONTAL,NULL); 

// This is the label widget which holds the hourlabel. 

TimeLabelHour=XtVaCreateManagedWidget("TimeLabelHour", 
xmLabelGadgetClass,BackgroundRC6, 
XtVaTypedArg,XmNlabelString,XmRString," 1" ,2, 

XmNuserData, &Time_of_day, 

NULL); 

// Here we create an arrow to increment the hour.. 

Widget HourUpArrow = XtVaCreateManagedWidget("HourUpArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROW UP, NULL); 

// Add the callback for Up Arrow 

XtAddCallback (HourUpArrow, XmNarmCallback, ChangeTimeHourCB, (XtPointer) 1); 

// Here we create an arrow to decrement the hour. 

Widget HourDownArrow = XtVaCreateManagedWidgetfHourDownArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROW_DOWN, NULL); 

// Add the callback for Down Arrow 

XtAddCallback (HourDownArrow, XmNarmCallback, ChangeTimeHourCB, (XtPointer) -1); 

// This is the Month label gadget. 

XtVaCreateManagedWidget(" Month :",xmLabelGadgetClass, 

BackgroundRC6,XmNorientation, XmHORIZONTAL,NULL); 

// This is the label widget which holds the Month label. 

TimeLabelMonth=XtVaCreateManagedWidget("TimeLabelMonth", 

xmLabelGadgetClass,BackgroundRC6, 
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XtVaTypedArg,XmNlabelString,XmRString," 1 ",2, 

XmNuserData, &Time_of_day, 

NULL); 

// Here we create an arrow to increment the Month 
Widget MonthUpArrow = XtVaCreateManagedWidgetf'MonthUpArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROWUP, NULL); 

// Add the callback for Up Arrow 

XtAddCallback (MonthUpArrow, XmNarmCallback, ChangeTimeMonthCB, (XtPointer) 1); 

// Here we create an arrow to deccrement theMonth 
Widget MonthDownArrow = XtVaCreateManagedWidget("MonthDownArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROW_DOWN, NULL); 

// Add the callback for Down Arrow 

XtAddCallback (MonthDownArrow, XmNarmCallback, ChangeTimeMonthCB, (XtPointer)-!); 


// This is the Year label gadget. 

XtVaCreateManagedWidget(" Year :",xmLabelGadgetClass, 

BackgroundRC6,XmNorientation, XmHORIZONTAL,NULL); 

// This is the label widget which holds the wind speed and it is modified 

// when the wind speed is changed. 

TimeLabelYear =XtVaCreateManagedWidget("TimeLabelYear", 
xmLabelGadgetClass,BackgroundRC6, 
XtVaTypedArg,XmNlabelString,XmRString," 1994",4, 

XmNuserData, &Time_of_day, 

NULL); 

// Here we create an arrow to increment the wind velocity.. 

Widget YearUpArrow = XtVaCreateManagedWidget("YearUpArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROW_UP, NULL); 

// Add the callback for Up Arrow 

XtAddCallback (YearUpArrow, XmNarmCallback, ChangeTimeYearCB, (XtPointer) 1); 

// Here we create an arrow to deccrement the Month. 

Widget YearDownArrow = XtVaCreateManagedWidget("MonthDownArrow", 
xmArrowButtonGadgetClass,BackgroundRC6, 

XmNarrowDirection, XmARROW DOWN, NULL); 

// Add the callback for Down Arrow 

XtAddCallback (YearDownArrow, XmNarmCallback, ChangeTimeYearCB, (XtPointer) -1); 

// Here manage the widget. 

XtManageChild (BackgroundRC6); 

// This the rowcolumn that holds Pull down Menu widget under it. 

n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_WIDGET); n++; 

XtSetArg(wargs[n],XmNlefLAttachment, XmATTACH FORM); n++; 
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XtSetArg(wargs[n],XmNrightAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNtopWidget, BackgroundRC6); n++; 
XtSetArg(wargs[n],XmNorientation, XmHORIZONTAL); n++; 
XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNpacking, XmPACKTIGHT); n++; 

Widget BackgroundRC7 = XmCreateRowColumn(Backgroundform2, 

"BackgroundRC7", wargs, n); 


// Add another pulldown to the RowColumn cloudiness sellection button 
n = 0; 

Widget pulldownlO = XmCreatePulldownMenu (BackgroundRC7, 

"pulldownlO", wargs, n); 


// Add the entries to the pulldownlO. 

make_pulldown_entry(pulldownlO, "OVERCAST ", CloudinessCB, (XtPointer) & SOIL); 
make_pulldown_entry(pulldownlO, "CLOUDY ", CloudinessCB, (XtPointer) & GRASS); 
make_pulldown_entry(pulldownlO, "CLEAR ", CloudinessCB, (XtPointer) & SEA); 


// Create the option menu that has the pulldown under it. 
n = 0; 

label_string = XmStringCreateSimple("CLOUDINESS:"); 

XtSetArg(wargs[n], XmNlabelString, label_string); n++; 

XtSetArg(wargs[n], XmNsubMenuId, pulldownlO); n++; 

Widget optionlO = XmCreateOptionMenu (BackgroundRC7, "optionMenulO", 

wargs, n); 


XmStringFree(label_string); 


// Here manage the widget. 
XtManageChild (optionlO); 


// This is the FOG DENSITY gadget. 

XtVaCreateManagedWidget("FOG DENSITY: ".xmLabelGadgetClass, 

BackgroundRC7,XmNorientation, XmHORIZONTAL,NULL); 


// Create a Frame widget 
n = 0; 

XtSetArg(wargs[n],XmNtopAttachment, XmATTACH_FORM); n++; 
XtSetArg(wargs[n],XmNrightAttachment, XmATTACHFORM); n++; 
XtSetArg(wargs[n],XmNbottomAttachment, XmATTACH FORM); n++; 
XtSetArg (wargs[n], XmNshadowThickness,2); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_OUT); n++; 

Widget FogFrame = XmCreateFrame (BackgroundRC7, "FogFrame", wargs, n); 
XtManageChild (FogFrame); 

// Crate the fog scale 
n = 0; 

XtSetArg (wargs[n], XmNshowValue, True); n++; 

XtSetArg (wargs[n], XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNminimum, 0); n++; 

XtSetArg (wargs[n], XmNmaximum, 100); n++; 

XtSetArg (wargs[n], XmNprocessingDirection, XmMAX ON RIGHT); n++; 
Widget Fogscale = XmCreateScale(FogFrame, "Fogscale", wargs, n); 
XtManageChild (Fogscale); 
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// Add the callbacks 

XtAddCallback (Fogscale, XmNvalueChangedCallback, FogDensityCB , (XtPointer) NULL); 
XtAddCallback (Fogscale, XmNdragCallback, FogDensityCB, (XtPointer) NULL); 

// Here manage the widget. 

XtManageChild (BackgroundRC7); 


}//Endof create_background_screen_widgets(. ).... 
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/************************************************************************* 

* make_pulldown_entry - * 

* The goal of this function is to simplify the construction of individual entries in the * 

* pulldown menu. The routine assumes that a menubar has already been created and a * 

* menupane created for that menubar. * 

************************************************************************/ 

void make_pulldown_entry( 

Widget menupane, // A menupane of the menubar. 

char *entrytext, // Display text for the pulldown entry. 

XtCallbackProc callback, // Name of procedure to call 

// when this menu entry is selected. 

XtPointer userdata) // Data to be sent the callback. 

{ 

Arg wargs[ 10]; // Same old args stuff. 

Widget button; // Temp to hold the widget (not managed). 

// Specify the character set (set up XmStrings). 

XmStringCharSet charset = (XmStringCharSet) XmSTRING_DEFAULT_CHARSET; 


// Create this particular entry in the pulldown menu, 
n = 0; 

XtSetArg(wargs[n], XmNlabelString,XmStringCreateLtoR(entrytext, charset)); n++; 
button = XmCreatePushButton (menupane, entrytext, wargs, n); 

// We are assuming no data structure needs to be passed in to the callback. 
XtAddCallback (button, XmNactivateCallback, callback, user data); 

// Manage the entry in the menupane. 

XtManageChild (button); 

} 
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void BackgroundlCB(Widget, XtPointer, XtPointer) 

{ 

// Call the display list for soil texturing 
BACKGROUND = Index + 3 ; 
glCallList(BACKGROUND); 


void Background2CB(Widget, XtPointer, XtPointer) 

{ 

// Call the display list for grass texturing 
BACKGROUND = Index + 1 ; 
glCallList(BACKGROUND); 


} 

void Backgrounds CB(Widget, XtPointer, XtPointer) 

{ 

// Call the display list for sea texturing 
BACKGROUND = Index ; 
glCallList(BACKGROUND); 

} 


void Background4CB(Widget, XtPointer user data, XtPointer) 

{ 


// Call the display list for asphalt texturing 
BACKGROUND = Index + 2; 
glCallList(BACKGROUND); 

} 
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void SensorCB(Widget, XtPointer userdata, XtPointer) 

{ 

float Temp,mu,total_radiance; 

// First Load the Sensor ID to SENSOR variable, 
int * SENSOR J'JO = (int *)user_data ; 

SENSOR = *SENSOR_NO; 

// If the sensor is IR then we must modify the color table. 
if(SENSOR == _IR) 

{ 

// Manage IR_scale widget so color sellection button will be visible.. 
XtManageChild(SensorRCl 1); 

// First Find the extinction cooficient value. 

mu = CalculateJExtinctionO; 

// Now modify the color table with current Firing Distance. 

for(int index = 0 ; index < VERTICE NUM ; index++) 

{ 

Temp = OriginalJempfindexHO] ; 
total_radiance = Find_Total_Radiance(Temp); 

Temp = Temp_At_Sensor(total_radiance,mu); 

Temperature[index][0] = Temp ; 
Fill_Up_Color_Table(index,MODE); 

}//End of for loop.. 

} // End of if clause.. 

if(SENSOR == _TV) 

{ 

// Manage IR_scale widget so it will be not displayed anymore 
XtUnmanageChild(SensorRCl 1); 

// Modify the environment color to dark gray. 
for(int j = 0 ; j < 2500 ; j++) 

{ 

Color [j][0] = 0.7 ; 

Color [j][l] = 0.7; 

Color lj][2] = 0.7 ; 

Color [j][3] = 1.0 ; 

} //End of for loop. 

} // End of if(TV) 

if(SENSOR = _LASER){// If color scale was displayed make it unvisible.. 

XtUnmanageChild(SensorRCll);} 

} //End of SensorCB. 
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// This callback is called when the target option menu is called. 

void targetCB(Widget, XtPointer user data, XtPointer) 

{ 

float Temp,mu,total_radiance; 

//Load the target ID. 

int *TARGET_NO = (int *)user_data ; 

Target = *TARGET_NO; 


// If the current sensor is IR then we must modify the color table. 
if(SENSOR = JR) 

{ 


// First Find the mu value. 

mu = Calculate_Extinction(); 

// Now modify the color table with previous Firing Distance. 

for(int index = 0 ; index < VERTICE_NUM; index++) 

{ 

Temp = OriginalJTemp[index][0] ; 
totalradiance = Find_Total_Radiance(Temp); 
Temp = Temp_At_Sensor(total_radiance,mu); 
Temperature[index][0] = Temp ; 

Fill Up_Color_Table(index,MODE); 

}//End of for loop.. 

} // End of if clause.. 

if(SENSOR == TV) 

{ 

for(int j = 0 ; j < 2500 ; j++) 

{ 

Color [j][0] = 0.7 ; 

Color [j][l]= 0.7; 

Color [j][2] = 0.7 ; 

Color [j] [3] = 1.0; 

} //End of for loop. 


> 
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// Here we are going to open the related target file with according to the 
// target selected by the user, and call the initilization function. 
switch(Target) 

{ 

case TANK: 

FillInInitialTargetArrays("TankInitialFile.dat"); 

XtManageChild(TargetRowColumn3); 

XtManageChild(TargetRowColumn4); 

XtManageChild(TargetBB5); 

XtManageChild(TargetBB6); 

XtUnmanageChild(TargetBB7); 

BACKGROUND = Index + 3 ; 
glCallList(BACKGROUND); 
break; 

case TRUCK: //This is not yet implemented... 

//FillInInitialTargetArrays("TruckInitialFile.dat"); 

XtManageChild(TargetRowColumn3); 

XtManageChild(TargetRowColumn4); 

XtManageChild(TargetBB5); 

XtManageChild(TargetBB6); 

XtUiunanageChild(TargetBB7); 

BACKGROUND = Index + 3 ; 
glCallList(BACKGROUND); 
break; 

case SHIP: 

FillInInitialTargetArrays("ShipInitialFile.dat"); 

XtManageChild(TargetRowColumn3); 

XtManageChild(TargetRowColumn4); 

XtManageChild(TargetBB5); 

XtManageChild(TargetBB6); 

XtUnmanageChild(TargetBB7); 

BACKGROUND = Index ; 
gICallList(BACKGROUND); 
break; 

case HELO: // This target is also not yet implemented. 

//FillInInitialTargetArrays("HeloInitialFile.dat"); 

XtManageChild(TargetRowColumn3); 

XtManageChild(TargetRowColumn4); 

XtManageChild(TargetBB5); 

XtManageChild(TargetBB6); 

XtUnmanageChild(TargetBB7); 

BACKGROUND = Index + 3 ; 
glCallList(BACKGROUND); 
break; 


case HOUSE: 

FillInInitialTargetArrays("BuildingInitialFile.dat"); 

XtUnmanageChild(TargetRowColumn3); 

XtUnmanageChild(TargetRowColumn4); 

XtUnmanageChild(TargetBB5); 
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XtUnmanageChild(TargetBB6); 
XtManageChild(T argetBB7); 
BACKGROUND = Index + 3 ; 
glCalIList(BACKGROUND); 
break; 


case PLANE: 

FillInInitialTargetArrays(" AircraftInitialFile.dat"); 

XtManageChild(TargetRowColumn3); 

XtManageChild(TargetRowColumn4); 

XtManageChild(TargetBB5); 

XtManageChild(TargetBB6); 

XtUnmanageChild(TargetBB7); 

BACKGROUND = Index + 10 ; 
glCallList(BACKGROUND); 
break; 


default: 
break; 


}//End of targetCB.... 
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/************************************************************************ 

* This is the Display Callback.This callback activates the display and * 

* ControlWindow widgets, and Unmanage the InputWindow widget. * 

************************************************************************/ 


void DisplayCB(Widget, XtPointer, XtPointer) 

{ 


float Temp,mu,total_radiance; 

// Remove the InputWindow widget so it will be unvisible 
XtUnmanageChild (InputWindow); 

// Manage DisplayWindow Widget Make it visible. 

XtManageChild (DisplayWindow); 

// Manage Main Window Widget. Make this window visible too. 
XtManageChild(ControlWindow); 

// Store the targets vertice number in VERTICE NUM so when we are modifying the color arrays 

// we don't try to modify the unnacessary part of the array. 

switch(Target) 

{ 

case TANK: 

VERTICENUM = Tankjndex; 
break; 

case TRUCK: 

//VERTICE_NUM = TruckJndex; 

break; 

case SHIP: 

VERTICENUM = Shipjndex; 

break; 

case HELO: 

//VERTICE_NUM = Helojndex; 
break; 

case HOUSE: 

VERTICE_NUM = Building_Index; 

break; 

case PLANE : 

VERTICE_NUM = Airplane_Index; 

break; 

}// End of switch statement 

// If the sensor is IR then we must modify the color table. 
if(SENSOR == JR) 

{ 

// First Find the Extinction coeficient value. 

mu = Calculate_Extinction(); 

// Now modify the color table with previous Firing Distance. 

for(int index = 0 ; index < VERTICE NUM; index++) 
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{ 


Temp = Original_Temp[index][0] ; 
total_radiance = Find_Total_Radiance(Temp); 
Temp = Temp_At_Sensor(total_radiance,mu); 
Temperature[index][0] = Temp ; 
Fill_Up_Color_Table(index,MODE); 

}//Endoffor loop.. 

} //End of if clause.. 

// Here reset the state control variables. 

DISPLAYENABLED = 1; 

SEARCHED ONES = 0; 

READONES = 0; 

OLDMAXX = 1; 

OLDMAXY = 1; 

MAX_X = 0; 

MAXY = 0; 

MAXRED = 0.0; 

ROW = 0; 

COLOMN = 900; 

Max_Contrast_Of_Scene = 0; 


}// End of Display CallBack.. 
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/************************************************************************ 
♦These following callbacks are not used in the program but since the interface has already 

* been designed the necessary callbacs are defined but no action has been taken.They kept 

* for future use. 

************************************************************************/ 

static void materialCB(Widget, XtPointer, XtPointer){;} 

static void TargetAreaCB(Widget, XtPointer, XtPointer){;} 

static void TargetHeadingCB(Widget, XtPointer, XtPointer){;} 

static void TargetSpeedCB(Widget, XtPointer, XtPointer){;> 

static void EngineStatusCB( Widget, XtPointer, XtPointer){;} 

static void FireStatusCB(Widget, XtPointer, XtPointer){;} 

static void IsolationStatusCB(Widget, XtPointer, XtPointer){;} 

static void RelativeApertureCB(Widget, XtPointer, XtPointer){;} 

static void ApertureDramCB(Widget, XtPointer, XtPointer){;} 

static void MagnificationCB(Widget, XtPointer, XtPointer){;} 

static void OpticalTransmittanceCB(Widget, XtPointer, XtPointer){;} 

static void ElectronicBandwidthCB(Widget, XtPointer, XtPointer){;} 

static void SensorHeightCB(Widget, XtPointer, XtPointer){;} 

static void DetectivityCB(Widget, XtPointer, XtPointer){;} 

static void DetectorElementsCB(Widget, XtPointer, XtPointer){;} 

static void FiringAngleCB(Widget, XtPointer, XtPointer){;} 

static void FiringDistanceCB(Widget, XtPointer, XtPointer){;} 

static void WindDirectionCB(Widget, XtPointer , XtPointer){;} 

static void AeresolCB(Widget, XtPointer, XtPointer){;} 

static void LongitudeCB(Widget, XtPointer, XtPointer){;} 

static void LatitudeCB(Widget, XtPointer, XtPointer){;} 

static void CloudinessCB(Widget, XtPointer, XtPointer){;} 


/***************♦******************************************************** 
* This is the Change Wind Speed Callback. 

************************************************************************/ 


static void ChangeWindSpeedCB(Widget, XtPointer userdata, XtPointer call_data) 

{ 

int incr = (int) userdata; 
winddata *Wind; 

char buff 8]; 

XtVaGetValues (WindSpeedLabel,XmNuserData,&Wind,NULL); 
Wind->value += incr; 
if(Wind->value < 0) Wind->value = 0 ; 
sprintf(buf,"%d",Wind->value); 

XtVaSetValues(WindSpeedLabel, 

XtVaTypedArg,XmNlabelString,XmRString,buf,strlen(buf), 

NULL); 


} 


/************************************************************************ 


♦ This is the Change Change Time Hour Callback. 
************************************************************************/ 


static void ChangeTimeHourCB(Widget, XtPointer user data, XtPointer) 

{ 


int incr = (int) user data ; 
timedata*TIME; 
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charbuf[8]; 

XtVaGetValues (TimeLabelHour,XmNuserData,&TIME,NULL); 

TIME->hour += incr; 

if(TIME->hour < 1) TIME->hour = 1 ; 

if(TIME->hour > 24) TIME->hour = 24 ; 

sprintf(buf,"%d",TIME->hour); 

XtVaSetValues(TimeLabelHour, 

XtVaTypedArg,XmNlabelString,XmRString,buf,2, 

NULL); 

HOUR = TIME->hour; 


} 


/♦♦♦I********************************************************************* 

* This is the Change Time Month Callback. 


* * * * * * * * * * * * 4c * ** * * * * * * * * * * * * * * * * * * * * * * * * * * * *** * * * * * * * * ** * * * * * * * * * * * * * * * */ 


static void ChangeTimeMonthCB(Widget, XtPointer userdata, XtPointer) 

{ 


int incr = (int) user data ; 
timedata*TIME; 
char buf[8]; 

XtVaGetValues (TimeLabelMonth,XmNuserData,&TIME,NULL); 

TEME->month += incr; 

if(TIME->month < 1) TIME->month = 1; 

if(TIME->month >31) TIME->month = 31 ; 

sprintf(buf,"%d",TIME->month); 

XtVaSetValues(TimeLabelMonth, 

XtVaTypedArg,XmNlabelString,XmRString,buf,2, 

NULL); 


} 


/**:M***********>M********>M********************************************* 


* This is the Change Time Year Callback. 

* * * * * * ** * * * ★* * * * * * ** * * * * * * * 4c ** * * * * * * * * * * 4c * * 4c * * * 4c * * ** * * * * 4c * * * * * * ** * * * * * * j 


static void ChangeTimeYearCB(Widget, XtPointer user data, XtPointer) 

{ 


int incr = (int) user data; 
timedata*TIME ; 
char bufI8]; 

XtVaGetValues (TimeLabelYear,XmNuserData,&TIME,NULL); 
TIME->year += incr; 

if(TIME->year < 1994) TIME->year = 1994 ; 
if(TIME->year > 2100) TIME->year = 2100 ; 
sprintf(buf,"%d",TIME->year); 

XtVaSetValues(TimeLabelYear, 

XtVaTypedArg,XmNlabelString,XmRString,buf,4, 

NULL); 


} 
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******************************************************************* 

* This is the Fog Density Callback. The scale value is divided to 30000 to make the 

* density smaller so the fog color doesn't cover the whole screen. 
************************* **************************************** *******/ 


static void FogDensityCB(Widget, XtPointer, XtPointer call data) 

{ 

XmScaleCallbackStruct * call_value = (XmScaleCallbackStruct *) call_data; 
density = (call value -> value)/30000.0 ; 

} 


/************************************************************************ 


* This is the IR Scale Callback. With the help of "Mode" variable the color that is going to 

* loaded onto the targets are desided. 

************************************************************************/ 


static void IR_ScaleCB(Widget, XtPointer user_data, XtPointer) 

{ 

int *IR_SCALE_NO = (int *)user_data ; 

MODE = <TR_SCALE_NO; 

} 
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/************************************************************************* 

* Here we CreateTheScales () for the display screen.There are three scales * 

* created. * 

************************************* *************** ********************* J 

void Create_The_Scales(Widget parent) 


{ 

Widget pulldown, // Row coloumn widget 

option , 
button ; 


Arg wargs[ 16]; // Args used with XtSetArg below. 

XmFontList fontlist; // Font List 


XmString titleStringl, //Titles for scales.. 
titleString2, 
titleString3 ; 

XFontStruct *font; // Ptr to a font structure. 


XmString label string; 


// Create a compound string for the scale title 
font = XLoadQueryFont (XtDisplay (toplevel), "spcl2xl2c"); 
fontlist = XmStringCreateFontList (font, charset); 
titleStringl = XmStringCreateLtoR ("Side View Angle ", charset); 
titleString2 = XmStringCreateLtoR ("Top View Angle ", charset); 
titleString3 = XmStringCreateLtoR ("Viewing Distance(km.)", charset); 

// Create the bulletin board widget 
n=0; 

XtSetArg (wargs[n], XmNbottomAttachment, XmATTACH FORM); n++; 
XtSetArg (wargs[n], XmNleftAttachment, XmATTACH FORM); n++; 
XtSetArg (wargs[n], XmNbottomOfifset, 2); n++; 
bulletin_boardl=XmCreateBulletinBoard(parent, "bulletinl", wargs, n); 
XtManageChild(bulletin_boardl); 


// Create a Frame widget to make it so we can see the size of the window, 
n = 0; 

XtSetArg (wargs[n], XmNshadowThickness,2 ); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_IN); n++; 

Widget ffamel = XmCreateFrame (bulletin boardl, "framel", wargs, n); 
XtManageChild (framel); 

// Set up arglist and create the scale 1 
n = 0; 

XtSetArg (wargs[n], XmNfontList, fontlist); n++; 

XtSetArg (wargs[n], XmNshowValue, True); n++; 
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XtSetArg (wargs[n], XmNtitleString, titleStringl); n++; 

XtSetArg (wargsfn], XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[nj, XmNmaximum, 360); n++; 

XtSetArg (wargsfn], XmNprocessingDirection, XmMAXONRIGHT); n++; 
scalel = XmCreateScale(framel, "scalel", wargs, n); 

XtManageChild (scalel); 

// Add the callbacks 

XtAddCallback (scalel, XmNvalueChangedCallback, rotateyCB, (XtPointer) NULL); 
XtAddCallback (scalel, XmNdragCallback, rotateyCB , (XtPointer) NULL); 


// Add a separator 
n = 0; 

XtSetArg (wargsfn], XmNorientation, XmVERTICAL); n++; 
sep = XmCreateSeparator(parent,"separator",wargs,n); 

XtManageChild(sep); 

// Create the bulletin board2 widget 
n=0; 

XtSetArg (wargsfn], XmNleftAttachment, Xm ATT ACH_ WIDGET); n++; 

XtSetArg (wargsfn], XmNleftWidget,bulletin_boardl); n++; 

XtSetArg (wargsfn], XmNbottomAttachment, XmATTACH FORM); n++; 

XtSetArg (wargsfn], XmNbottomOffset, 2); n++; 

bulletin_board2= XmCreateBulletinBoard(parent, "bulletin2", wargs, n); 

XtManageChild(bulletin_board2); 

// Create a Frame widget 
n = 0; 

XtSetArg (wargsfn], XmNshadowThickness,2); n++; 

XtSetArg (wargsfn], XmNshadowType,XmSHADOW_IN); n++; 
frame2 = XmCreateFrame (bulletin_board2, "frame2", wargs, n); 

XtManageChild (frame2); 

// Set up arglist and create the scale2 
n = 0; 

XtSetArg (wargsfn], XmNfontList, fontlist); n++; 

XtSetArg (wargsfn], XmNshowValue, True); n++; 

XtSetArg (wargsfn], XmNtitleString, titleString2); n++; 

XtSetArg (wargsfn], XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargsfn], XmNmaximum, 90); n++; 

XtSetArg (wargsfn], XmNminimum, 0); n++; 

XtSetArg (wargsfn], XmNprocessingDirection, XmMAX_ON_RIGHT); n++; 
scale2 = XmCreateScale(frame2, "scale2", wargs, n); 

XtManageChild (scale2); 

// Here we set the initial value of this scale since it is different than 
// minvalue of ths scale. 

XmScaleSetValue(scale2 ,int(Firing_Angle)); 

// Add the callbacks 

XtAddCallback (scale2, XmNvalueChangedCallback, rotatexCB, (XtPointer) NULL); 
XtAddCallback (scale2, XmNdragCallback, rotatexCB, (XtPointer) NULL); 
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// Add a separator 
n = 0; 

XtSetArg (wargs[n], XmNorientation, XmVERTICAL); n++; 
sep = XmCreateSeparator(parent,"separator",wargs,n); 
XtManageChild(sep); 


// Create the bulletin board3 widget 
n=0; 

XtSetArg (wargs[n], XmNleftAttachment, XmATTACHWIDGET); n++; 
XtSetArg (wargs[n], XmNleftWidget,bulletin_board2); n++; 

XtSetArg (wargs[n], XmNbottomAttachment, XmATTACH FORM); n++; 
XtSetArg (wargs[n], XmNbottomOffset, 2); n++; 
bulletin_board3=XmCreateBulletinBoard(parent, "bulletin3", wargs, n); 
XtManageChild(bulletin_board3); 


// Create a Frame widget to make it so we can resize the window, 
n = 0; 

XtSetArg (wargsfn], XmNshadowThickness,2 ); n++; 

XtSetArg (wargs[n], XmNshadowType,XmSHADOW_IN); n++; 

Widget frame3 = XmCreateFrame (bulletin_board3, "frame3", wargs, n); 
XtManageChild (frame3); 

// Set up arglist and create the scale3 
n = 0; 

XtSetArg (wargs[n], XmNfontList, fontlist); n++; 

XtSetArg (wargs[n], XmNshowValue, True); n++; 

XtSetArg (wargs[n], XmNtitleString, tit!eString3); n++; 

XtSetArg (wargs[n], XmNorientation, XmHORIZONTAL); n++; 

XtSetArg (wargs[n], XmNminimum, 0); n++; 

XtSetArg (wargs[nj, XmNvalue, 400); n++; 

XtSetArg (wargs[n], XmNmaximum, 400); n++; 

XtSetArg (wargs[n], XmNdecimalPoints, 2); n++; 

XtSetArg (wargs[n], XmNprocessingDirection, XmMAX_ON_RIGHT); n++; 
scale3 = XmCreateScale(frame3, "scale3", wargs, n); 

XtManageChild (scale3); 

// Add the callbacks 

XtAddCallback (scale3, XmNvalueChangedCallback, translatezCB, (XtPointer) NULL); 
XtAddCallback (scale3, XmNdragCallback, translatezCB, (XtPointer) NULL); 

// Here we set the initial value of this scale since it is different than minvalue of this scale. 
XmScaleSetValue(scale3 ,400); 

// Get rid of the fontlists 
XmFontListFree (fontlist); 

// Get rid of the Stringlists. 

XmStringFree (titleStringl); 

XmStringFree (titleString2); 

XmStringFree (titleString3); 
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// Create widget to hold the push buttons.. 
n = 0; 

XtSetArg (wargs[n], XmNleftAttachment, XmATTACHWIDGET); n++; 
XtSetArg(wargs[n], XmNleftWidget,bulletin_board3); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNorientation, XmVERTICAL); n++; 
rc2 = XmCreateRowColumn(parent,"degree",wargs,n); 
XtManageChild(rc2); 


// Create the bulletin board widget for INPUTWINDOW push button. 
n=0; 

BboardForInput= XmCreateBulletinBoard(rc2, "BboardForlnput", wargs, n); 
XtManageChild(BboardForInput); 

// Now add an Input PushButton. 

label string = XmStringCreateLtoR("INPUTWINDOW", charset); 
n = 0; 

XtSetArg(wargs[n], XmNalignment, XmALIGNMENT CENTER); n++; 

XtSetArg(wargs[n], XmNlabelString, label string); n++; 

button = XmCreatePushButton(BboardForInput, "InputScreen", wargs, n); 

XtManageChild(button); 

XmStringFree(label_string); 

// Now add an Exit callback for the PushButton. 

XtAddCallback(button, XmNarmCallback, BackToInputScreenCB, (XtPointer) NULL); 

// Create widget to hold the push buttons., 
n = 0; 

XtSetArg (wargs[n], XmNleftAttachment, XmATTACHWIDGET); n++; 

XtSetArg (wargs[n], XmNleftWidget,rc2); n++; 

XtSetArg (wargs[n], XmNnumColumns, 1); n++; 

XtSetArg (wargs[n], XmNorientation, Xm VERTICAL); n++; 

Widget rc3 = XmCreateRowColumn(parent,"degree",wargs,n); 

XtManageChild(rc3); 

// Create the bulletin board widget for Exit push button. 
n=0; 

BboardForExit= XmCreateBulletinBoard(rc3, "BboardForExit", wargs, n); 
XtManageChild(BboardForExit); 

// Now add an Exit PushButton. 
label string = XmStringCreateLtoR("EX3T ", charset); 
n = 0; 

XtSetArg(wargs[n], XmNalignment, XmALIGNMENT CENTER); n++; 
XtSetArg(wargs[n], XmNlabelString, label_string); n++; 
button = XmCreatePushButton(BboardForExit, "exit", wargs, n); 
XtManageChild(button); 

XmStringFree(labelstring); 

// Now add an Exit callback for the PushButton. 

XtAddCallback(button, XmNarmCallback, quitCB, (XtPointer) NULL);} 

// End of the create the scales. 
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/****************************************************************************** 

* This initial calback where it iscalled when the program reaches the definition of the 
*"glw" display widget. What this callback does is it initializes the window routines 

* display lists for texturing and the cube which is used for sky polygons. It also loads the 

* initial values of some GL and user defined variables. 

*********************************************************************** ******/ 
static void initCB (Widget w, XtPointer, XtPointer calldata) 

{ 

// Get us a GLwDrawingAreaCallbackStruct ptr to the call data. 

GLwDrawingAreaCallbackStruct ♦glptr = (GLwDrawingAreaCallbackStruct *) call_data; 

Arg wargs[ 1 ]; // Arg temp. 

XVisuallnfo *vi; // Pointer to XVisuallnfo. 

GLint gdtmp; // Used in the get capability stuffbelow. 


// Get the visual info... 

XtSetArg(wargs[0], GLwNvisuallnfo, &vi); 

XtGetValues(w, wargs, 1); 

// Create a new GLX rendering context. 

// GLTRUE -> Specify direct connection to graphics system (if possible), 
glxcontext = glXCreateContext(XtDisplay(w), vi, 0, GL TRUE); 

// Set the global window and display, 
globaldisplay = XtDisplay(w); 
global_window = XtWindow(w); 

// Make this drawing area the current one. 
GLwDrawingAreaMakeCurrent(w, glxcontext); 

// Test to see if this machine has a z-buffer. 
if((glGetIntegerv(GL_DEPTH_BITS, &gdtmp), gdtmp) == 0) 

{ 

cerr « "This machine does not have a hardware zbufifer" « endl; 
exit(0); 

} 

// Turn on z-buflfering. 

glEnable(GL_DEPTH_TEST); 


// Turn on Gouraud shading. 
glShadeModel(GL_SMOOTH); 


// Get a font we can use to display our messages. We usebthe first font to display "LOCK ON" 

// and the second one to display numbers on the color scale... 

fontHandlel = makeRasterFont(w, "-*-Umes-mediuiri-r-normal~100-*-*-*-p-84-iso8859-r'); 
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fontHandle2 = makeRasterFont(w, "-*-times-medium-r-normal~14-*-*-*-p-84-iso8859-l"); 


// Set the Viewport for the first draw of the window. 
gIViewport (0,0, glptr->width, glptr->height); 

// Set the initial viewing parameters. 
set_initial_viewing_values(); 

// Load the default target arrays here 
FillInInitialTargetArrays("TankInitialFile.dat''); 

// Call the display list to form the texturing list 
CallDisplayListForTexturingO; 


// Make a display List for Cube. 
CallDisplayListForCube(); 

// Load the default initial values to the variables 
MOUNTAIN = Index + 6 ; 

CLOUDS = Index + 5 ; 

BACKGROUND = Index+3 ; 


// Set the texture environment we want. 
setGLTextureEnvironmentO; 

// Unmanage the info screen so it wnt be visiuble any more... 
XtUnmanageChild (info); 

} // end of initCBO 
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/************************************************************************* 

* exposeCB - * 

* This routine is called whenever the window is uncovered. * 

*************************************************************************/ 
static void exposeCB(Widget w, XtPointer, XtPointer) 

{ 

// Set the window into which GL drawing should be done. 

GLwDrawingAreaMakeCurrent(w, glxcontext); 

if(DISPLAY_ENABLED){ draw_the_scene(); } // otherwise we are in the input screen so 

// don't spent time by calling display window. 

} 


/************************************************************************* 

* resizeCB - * 

* This routine is called whenever the window is moved or resized. * 

*********************** j^***>M*******************************************/ 

static void resizeCB (Widget w, XtPointer, XtPointer call_data) 

{ 

// Get us a GlwDrawingAreaCallbackStruct ptr to the call_data. 

GLwDrawingAreaCallbackStruct *glptr = (GLwDrawingAreaCallbackStruct *) calldata; 


// Set the window into which GL drawing should be done. 
GLwDrawingAreaMakeCurrent(w, glxcontext); 

// Set the viewport using the window size currently set. 
glViewport (0, 0, (GLsizei) glptr->width, (GLsizei) glptr->height); 


if(DISPLAY_ENABLED){ draw_the_scene(); } // otherwise we are in the input screen so 

// don't spent time by calling display window. 


} 
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/*********************S^**!M*********************************************** 


* inputCB - 

* This routine handles all types of input from the GL widget. ' 

* In this particular example, the KeyRelease handles the 

* Escape-Key so that its release exits the program. 1 

* * * * * * * * * * * * ** * * * * * * * * * * 4c * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * % * * * * * % * 4c * j|y 


static void inputCB (Widget, XtPointer, XtPointer calldata) 

{ 


// Get us a GLwDrawingAreaCallbackStruct ptr to the call data. 
GLwDrawingAreaCallbackStruct *glptr = (GLwDrawingAreaCallbackStruct *) call data; 

char ascii[l]; // Hold for the ascii chars returned from XLookupString. 

int nchars; // Number of characters returned from XLookupString. 

KeySym keysym; // The X version of the returned character. 

XKeyEvent *ptr; // ptr to the Key Event structure. 


// We look at the type of event to see if we should pay attention... 
// In this example, we only pay attention to key release events. 
switch(glptr->event->type) 

{ 


case KeyRelease: 

// We must convert the keycode to a KeySym before it is possible 

// to check if it is an escape. We also dump the Ascii into the 

// array ascii. The return value from XLookupString is the 

// number of characters dumped into array ascii. 

ptr = (XKeyEvent *) glptr->event; 

nchars = XLookupString(ptr, ascii, 1, &keysym, NULL); 

if(nchars == 1 && keysym — (KeySym) XKEscape) 

{ 

// We have an escape. Time to exit. 
exit(O); 

} 

break; 


// This part of the program kept in the program forseeing that in the future some functions will be 
// assigned to the mouse buttons, 
case ButtonPress: 

// We have a button press from the mouse. 

// Which mouse button was it? 
switch(glptr->event->xbutton.button) 

{ 

case Button2: // The middle button was pressed. 
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break; 

case Button 1: // The left button was pressed, 
break; 

case Button3: // The right button was pressed, 
break; 

default: 

break; 

} // end switch on which mouse button it was. 
break; 

case ButtonRelease: 

// We have a button release from the mouse. 

// Which mouse button was it? 
switch(glptr->event->xbutton.button) 

{ 

case Button2: // The middle button was pressed, 
break; 

case Button 1: // The left button was pressed, 
break; 

case Button3: // The right button was pressed, 
break; 


default: 

break; 

} // end switch on which mouse button it was. 
break; 

case MotionNotify: 

// We have motion on the mouse. We only get notified of motion when a mouse button is 
// pressed, 
break; 
default: 
break; 

} // end of switch on event->type. 

} // end of inputCB. 

/****************************************************************************** 

* quitCB - * 

* This function is called whenever the "Quit" menu option is selected. * 

*****************im***********************************************************/ 

static void quitCB (Widget w, XtPointer, XtPointer) 

{ 

//Bye, bye... 

XtCloseDisplay(XtDisplay(w)); 

exit(O); 

} // End of quitCB 

/**************************************jM*************************************** 

* rotatexCB * 

* This function is called whenever the Top view scale is moved. * 

******************************************************************************/ 



static void rotatexCB(Widget, XtPointer, XtPointer calldata) 

{ 


float Temp,mu,total_radiance,Extincted_temp; 

XmScaleCallbackStruct * callvalue = (XmScaleCallbackStruct *) call data; 

pitch = call_value -> value ; 

FiringAngle = pitch; 

if(SENSOR== JR) 

{ 


// First Find the mu value. 

mu = Calculate_Extinction(); 


// Now modify the color table with previous Firing Angle. 
for(int index = 0 ; index < VERTICE NUM; index++) 

{ 

Temp = OriginalJemp[index][0] ; 
total_radiance = Find_Total_Radiance(Temp); 
Temp = Temp_At_Sensor(total_radiance,mu); 
Temperature[index][0] = Temp ; 

Fill JJp_Color_Table(index,MODE); 

}//End of for loop.. 

} // End of if clause.. 

// Modify the statre control variables. 

SEARCHED ONES = 0; 

READONES = 0; 

ROW = 0; 

COLOMN = 900 ; 

OLD MAX X = MAX_X; 

OLD_MAX_Y = MAX_Y; 

MAXRED = 0.0; 

Max Contrast Of Scene = 0 ; 

} // End of rotatexCB.... 


/******^***************************************************************** 

* rotateyCB 

* This function is called whenever the Side view scale is moved. 
************************************************************************/ 
static void rotateyCB(Widget, XtPointer, XtPointer call_data) 

{ 


XmScaleCallbackStruct * call_value = (XmScaleCallbackStruct *) call_data; 
pitch = call_value -> value; 


} 
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/******************************************************************************* 

* translatezCB * 

* This function is called whenever the Distance scale is moved. * 

>******************************************************************************/ 


static void translatezCB(Widget, XtPointer, XtPointer calldata) 

{ 


float Temp,mu,kt,total_radiance; 

XmScaleCallbackStruct * call_value = (XmScaleCallbackStruct *) call data; 

tz = callvalue -> value; 

FiringDistance = tz; 

// If the sensor is IR then we must modify the color table. 
if(SENSOR — _IR) { 

// First Find the mu value. 

mu = Calculate_Extinction(); 

// Now modify the color table with previous Firing Distance, 
forfint index = 0 ; index < VERTICE NUM ; index++) 

{ 

Temp = Original_Temp[index][0] ; 
totalradiance = Find_Total_Radiance(Temp); 

Temp = Temp_At_Sensor(total_radiance,mu); 
Temperature[index][0] = Temp ; 
Fill_Up_Color_Table(index,MODE); 

}//End of for loop.. 

} //End of if clause.. 


// Modify the state control variables.... 

SEARCHED ONES = 0; 

READ_ONES = 0; 

ROW = 0; 

COLOMN = 900 ; 

OLD_MAX_X = MAX_X; 

OLD MAX Y = MAX Y; 

MAXRED = 0.0; 

Max_Contrast_Of_Scene = 0; 

}//End of translatezCB... 

mt+m+m************* ************ ********************************** ***** **** 

* drawWP - * 

* This function is called by the work procedure. It is called repeatedly whenever there are * 

* no events to process. This function returns FALSE so that the work procedure does NOT * 

* stop calling it. * 

******+*+*+***++*++***+***+**************+***++********************************/ 
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GLboolean drawWP() 

{ 

if(DISPLAY_ENABLED) { draw_the_scene(); } 

// If we do not return FALSE, the work procedure will stop calling this function. 
retum(FALSE); 


} 
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y******************************************************************************* 

* drawthescene * 

* This function draws the picture. * 

******************************************************************************/ 


void draw_the_scene() 

{ 


// Turn on z-buffering! 

glEnable(GL_DEPTH_TEST); 


// Set the background color to cyan (RGBA) and clear the z-buffer. 
glClearColor(0.25, 0.0, 1.0, 1.0); 

glClear(GL_DEPTH_BUFFER_BIT | GLCOLORBUFFERBIT); 


//Set the projection matrix. 

// The near and far values are distances from the viewer 
// to the near (1.0) and far (10000.0) clipping planes. 
glMatrixMode(GLPROJECTION); 
glLoadldentityO; 

gluPerspective(45.0, 1.25, 1.0, 18000.0); 


// Load the ModelView matrix with a unit matrix. 

glMatrixMode(GLMODELVIEW); 

glLoaidldentityO; 


if(SENSOR == _TV) 

{ 

// Turn on the Light Model... 
tumOnTheLightModel(); 

// Turn on the Lights ... 
tumOnTheLights(); 

//Enable the fog for the target. 
EnableFog(density); 


} 


//Here we check the position of the viewer.. 
load_the_viewpoint(); 

// We are at the viewpoint and looking towards 
// the reference point of the object. 

// Up vector is the vector orienting our view volume around 
// the line of sight. 

gluLookAt(Viewpoint[0], Viewpoint^], Viewpoint^], // View point. 

Refpoint[0], Refjx>int[l], Refjx>int[2], // Ref point, point we are looking towards. 

0.0, 1.0, 0.0); // Up vector. 
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// Concatenate onto the ModelView matrix the results of our 
// accumulative transformations. 

// The coordinate passed in is the point about which 
// we would like our operations (rotations/scales) to occur. 

// Think of this point as the center of the displayed object. 

concatenate_accumulative_matrices(0.0,0.0,0.0); 


//First draw the target.That is the nearest object to the user. 
drawTarget(Target,0.0,-2.0,0.0); 

if(SENSOR == _TV) 

{ 

// Draw the target shadow .. 
drawTargetShadow(Target,0.0,-2.0,0.0); 

// Turn on texturing... 
glEnable(GL_TEXTURE_2D); 


} 


// Here we are drawing the ground and calling the ground texturing. 
glCallList(BACKGROUND); 
drawGround(0.0, -2.0, 0.0,BACKGROUND,MODE); 

// Here we are drawing the mountains and calling the mountain texturing. 
glCallList(MOUNTAIN); 

DrawTheMoutains(MODE); 


//Call the display list for clouds texture and Cube. 
glCallList(CLOUDS); 
glCallList(CUBE); 

//Here we draw the sun.The position of the sun is time dependent. 
DrawTheSun(HOUR); 


if(SENSOR — _TV) 

{ 

//Turn off everything being turned on. (Save energy ) 

glDisable(GL_TEXTURE_2D); 

DisableFogO; 

tumOffTheLights(); 

turnOffTheLightModelO; 

Call_TV_Functions(); 

}//End of if(TV).. 


Ill 




if(SENSOR == _IR) { Call_IR_Functions();} 

// Swap the buffers as we are doing double buffering. 
glXSwapBuffers(globaI_display, global window); 

} // end of draw the scene 
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/******************************************************************************* 

* The following routine actually does the matrix multiplication to compute * 

* the following formula. * 

* * 

* P' = P . T(to origin). S(acc). R(x acc). R(y acc). R(z acc) * 

* . T(to acc. loc). T(back to specified center).. lookat() ,perspective() * 

* * 

* The input coordinate is the point about which operations should be done. * 

* Think of that point as the point translated to the origin. * 

******************************************************************************/ 
void concatenate_accumulative_matrices(float refx, float refy, float refz) 

{ 

// Since the operation is done by multiplying the matrices,the order 
// of operation is important. 

// We are going to use the ModelView stack. 
glMatrixMode(GLMODELVIEW); 

// translate center of object back to original location. 
glTranslateflrefx, refy, refz); 

glRotatef((pitch + 1.0), 1.0,0.0,0.0); // Rotate to look from the top 

glRotatef(heading, 0.0, 1.0, 0.0); // Rotate to lok from each side.. 

// translate center of object to the origin. 
glTranslatef(-refx, -refy, -refz); 


} 


/******>m*********>i^****************************************************** 

* The following routine loads a unit matrix into the input array. 
************************************************************************/ 

void unit(GLfloat *m) 

{ 

static GLfloat un[4][4] = { 1.0,0.0,0.0,0.0, 

0 . 0 , 1 . 0 , 0 . 0 , 0 . 0 , 

0 . 0 , 0 . 0 , 1 . 0 , 0 . 0 , 

0 . 0 , 0 . 0 , 0 . 0 , 1.0 }; 

long i j; 

GLfloat *ptr; // Get a temp ptr. 


ptr = m; // Get a ptr to the start of m. 

// copy the matrix elements... 
for(i=0; i < 4; i=i+l) 

{ 

for(j=0; j <4; j=j+l) 

{ 


*ptr++ = un[j][i]; 




} 


}// End oj unit(..) 


/************************************************************************ 

* This function sets the initial viewing parameters. * 

*************>m*********************************************************/ 


void set_initial_viewing_values() 

{ 


// Set the initial angles, 
heading = 0.0; 
pitch = Firing_Angle; 
tz = FiringDistance; 

roll = 0.0; //This value is not used but saved for future implementations which can 
// be modified in rotatezCB(....) 
density = 0.0 ;//This value is used for fog implimentation 

Viewpoint[0] = 0.0; 

Viewpoint 1] = 0.0 ; 

Viewpoint^] = tz+5; // We don't want get in the targets which is at the origin.. 

// We are loo,ing to the center.. 

Refpoint[0] = 0.0; 

Refpoint[ 1 ] = 0.0 ; 

Refpoint[2] = 0.0 ; 


} 


/************************************************************************* 
* This function sets the initial viewing parameters. 

f***^**********^************^****************************************/ 


void load_the_viewpoint() 

{ Viewpoint[0] = 0.0 ; Viewpoint 1] = 0.0 ; Viewpoint[2] = tz+5;} 


* 


/* 4 ^********************************************************************* 

* This callback is used to return back to input screen from the display 

* screen. This callback manages the InputScreen widget, and Unmanage the * 

* DisplayScreen and ControlScreen widgets. 

*************4^****4^***************:t***********************************/ 

static void BackToInputScreenCB(Widget, XtPointer, XtPointer) 

{ 


// Remove the DisplayWindow widget 
XtUnmanageChild (DisplayWindow); 

// Remove the ControlWindow widget 
XtUnmanageChild (ControlWindow); 


// Manage Main Window Widget 
XtManageChild(InputWindow); 
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// Here set the global DISPLAYENABLED variable to zero so when draw the 
// scene is called no function is executed.. 

DISPLAYENABLED = 0; 


} 


/************************************************************************ 

* This function reads-in the target data files and loads the vertex temperature and color * 

* arrays. Function return type is integer for error checking purposes. * 

************************************************************************/ 


int FillInInitiaITargetArrays(char * InputFile) 

{ 


int Verticeno; 
float Vertexx, 

Vertex_y, 

Vertexz, 

Temp , 

kt, 

mu, 

total_radiance = 0.0; 
const int size = 100 ; 

char buffer[size]; // Used to store unnecassaiy characters. 

ifstream from (InputFile); // Input Data File 

if ((from) 

{ 

cout«" I can't open"«InputFile «"File"« endl; 
return 0; 

} 

// Since this function is called when a target is initialized find the "mu" value according to either 
// default or sellected values of firin angle, 
mu = Calculate_Extinction(); 

while (Ifrom.eofO) 

{ 

from » Verticeno; 
from » Vertexx; 
from » Vertex_y; 
from » Vertex_z; 
from » Temp ; 
from. getline(buffer,size,V); 

Temp = Temp ; 

Original_Temp[ Verticeno] [0] 

Vertex[ Verticeno] [0] 

Vertex! Vertice_no] [ 1 ] 

Vertex[Vertice_no]f2] 

// We read the temperature values then calculate theradiance of the target for this temperature 
// value assuming the target is a black body radiator. 


// Go to next line 

= Temp ; 

= Vertexx; 

= Vertex_y; 

= Vertexz; 
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kt=Boltzman_Constant*(Temp); 
totalradiance = FindTotalRadiance(Temp); 


// Now with this radiance and mu extinction coef. find the equilant temp on the sensor. This 
// function applies necessary extinction to the signal and after finding the new level of the signal 
// it transforms it to temperature value again.This calculated temp, is loaded to the current temp. 
// array. 

Temp = Temp_At_Sensor(total_radiance,mu); 

Temperature[Vertice_no][0]= Temp; 

// Here the temperature values transformed to the equilant color values, and color array is loaded. 
Fill_Up_Color_Table(Vertice_no,MODE); 


} // end while... 

from.close(); 
return 1; 

}// End of function... 


f****************************************************************************** 

* This function calculte the extinction coeficient with respect to the angle that comes from * 

* " Top View Angle " scale. We used the complimentary angle for this calculation * 

* because the top view is 90 degree for scale but 0 degree for LOWTRAN (LT.)codes * 

* This extinction values can be found directly running the LT. itself but it is very * 

* time consuming so we couldn't afford to spend that much time for a simulation that the * 

* time is very critical . Instead we calculated the extinction coefficient by running LT* * 

* with five degree increments and used those values here. * 

**********************^****sm************************************************/ 


float CalculateJExtinctionO 

{ 


float App_Angle= 90.0 - Firing_Angle; 

if (AppAngle <= 90 && AppAngle > 85) return 1.1660 
if (App_Angle <= 85 && App Angle > 80) return 0.8099 
if (App Angle <= 80 && App Angle > 75) return 0.6134 
if (App Angle <= 75 && App Angle > 70) return 0.4650 
if (App Angle <= 70 && App Angle > 65) return 0.3518 
if (App Angle <= 65 && App Angle > 60) return 0.2702 
if (App Angle <= 60 && App Angle > 55) return 0.2123 
if (App Angle <= 55 && App Angle > 50) return 0.1251 
if (App_Angle <= 50 && App Angle > 45) return 0.1239 
if (App Angle <= 45 && App_Angle > 40) return 0.1217 
if (App_Angle <= 40 && App_Angle > 35) return 0.1068 
if (App Angle <= 35 && App Angle > 30) return 0.0960 
if (App Angle <= 30 && App Angle > 25) return 0.0877 
if (App_Angle <= 25 && App_Angle > 20) return 0.0816 
if (App Angle <= 20 && App Angle > 15) return 0.0739 
if (App Angle <= 15 && App Angle > 10) return 0.0739 
if (App Angle <= 10 && App Angle > 5) return 0.0716 ; 
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if (App_Angle <= 5 && AppAngle > 0) return 0.0703 ; 
if (App Angle == 0) return 0.0703 ; 


} 


J* *************************** *********************** ********************* 

* This function finds the temperature equilant color for three different coloring schemes* 

* and loads the color array. * 

************************************************************************/ 


int Fill_Up_Color_Table(int Vertice_no,int COLOR_MODE) 

{ 

double red; 

if(Temperature[Vertice_no][0] <= 10) //Dark Blue 

{ 

red = 0.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);) 
else { 

Color[Vertice_no][0]= 0.0; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 1.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1; 

} 

if(Temperature[Vertice_no][0] >10 && Temperature[Vertice_no][0] <= 15)// Dirty_Green 

{ 

red = 1.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Amy_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.0625; 

Color[Vertice_no][l]= 0.3137; 

Color[Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertice_no][0] >15 && Temperature[Vertice_no][0] <= 20)//Light_Green 


red = 2.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 
{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.1250; 
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Color[Vertice_no][l]= 0.8784; 

Color[Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

> 

return 1 ; 

} 

if(Temperature[Vertice_no][0] >20 && Temperature[Vertice_no][0] <= 25)//Light_Blue 

{ 

red = 3.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

(ModifyColorArrayWithAGiven_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.1875; 

Color[Vertice_no][l]= 1.0; 

Color[Vertice_no][2]= 1.0; 

ColortVertice_no][3]= 1.0; 

} 

return 1; 

} 

if(Temperature[Vertice no][0] >25 && Temperature[Vertice_no][0] <= 30)//Brown 

{ 

red - 4.0/16.0; 

if( COLOR MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.25; 

Color[Vertice_no][l]= 0.25; 

Color[Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertice_no][0] >30 && Temperature[Vertice_no][0] <= 35)// Purplish 

{ 

red = 5.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_AiTay_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.3125; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 1.0; 

Color[Vertice_no][3]= 1.0; 

} 


return 1; 
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} 


if(Temperature[Vertice_no][0] >35 && Temperature[Vertice_no][0] <= 40) //Gray 

{ 

red = 6.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.3750; 

Color[Vertice_no][l]= 0.3750; 

Color[Vertice_no][2]= 0.3750; 

Color[Vertice_no][3]= 1.0; 

} 

return 1; 

} 

if(Temperature[Vertice_no][0] >40 && Temperature!Vertice_no][0] <= 45)//Dark_Brown 

{ 

red = 7.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[ Vertice_no][0]= 0.4375; 

Color[Vertice_no] [ 1]= 0.5843; 

Color! Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1; 

} 

if(Temperature[Vertice_no][0] >45 && Temperature[Vertice_no][0] <= 50) // Cherry 

{ 

red = 8.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no] [0]= 0.5; 

Color[Vertice_no][l]= 0.0; 

Color! Vertice_no][2]= 0.5; 

Color! Vertice_no] [3]= 1.0; 

} 


return 1; 

} 

if(Temperature!Vertice_no][0] >50 && Temperature[Vertice_no][0] <= 55) //Fading_Red 

{ 

red = 9.0/16.0; 
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if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color} Vertice_no]}0]= 0.5625; 

Color[Vertice_no][l]= 0.1568; 

Color[Vertice_no] [2]= 0.5607; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertice_no][0] >55 &&Temperature[Vertice_no][0] <=60) // Purple 

{ 

red =10.0/16.0; 

if( COLOR MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.6250; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 1.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertice_no]}0] >60 && Temperature[Vertice_no]}0] <= 65)//Dark_Gold 

{ 

red =11.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color} Verticeno] [0]= 0.6875; 

Color[Vertice_no][l]= 0.8156; 

Color} Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature}Vertice_no][0] >65 && Temperature[Vertice_no][0] <= 70)//Dirty_Red 

{ 

red= 12.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color} Vertice_no]}0]= 0.75; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 0.3960; 

Color[Vertice_no][3]= 1.0; 

} 
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return 1 ; 

} 


if(Temperature[Vertice_no][0] >70 && Temperature[Vertice_no][0] <= 75)//Dirty_Pink 

{ 

red = 13.0/16.0; 

if( COLOR MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.8125; 

Color[Vertice_no][l]= 0.3725; 

Color[Vertice_no][2]= 1.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertiee_no][0] >75 && Temperature[Vertice_no][0] <= 80)//Pinkish 

{ 

red = 14.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modiiy_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.8750; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 1.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(Temperature[Vertice_no][0] >80 && Temperature!Vertice_no][0] <= 85) // Gold 

{ 

red = 15.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
else { 

Color[Vertice_no][0]= 0.9375; 

Color[Vertice_no][l]= 0.9490; 

Colorf Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

if(TemperaturelVertice_no][0] >85)// 85<T Degrees Dark_Red 

{ 

red = 16.0/16.0; 

if( COLOR_MODE != SYCLIC_SCALE) 

{Modify_Color_Array_With_A_Given_Mode(Vertice_no,COLOR_MODE, red);} 
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else { 

Color[Vertice_no][0]= 1.0; 

Color[Vertice_no][l]= 0.0; 

Color[Vertice_no][2]= 0.0; 

Color[Vertice_no][3]= 1.0; 

} 

return 1 ; 

} 

return 1; 

}// End of function Fill_Up_Color_Table 

/****************************************************************************** 

* This function reads the display screen pixels from the frame buffer with a 100 by 100 

* square and finds the avarage red value of square and finally after reading the last square 

* of the screen it finds the maximum red value of the screen . 
******************************************************************************/ 
void Read_Screen() 

{ 

float sum_of_red = 0.0,avr_red = 0.0; 
int index = 0, deb = 0; 


for (int y = 900 ;y>=0;y=y-100) 

{ 

for(int x = 0;x <=1000 ; x= x + 100) 

{ 

glReadPixels (x,y, 100,100,GL_RED,GL_UNSIGNED_BYTE,bufferl); 

for (int counter = 0;counter < 10000 ;counter++) 

{ 

// We defined the bufferl array as character to make it occupy less memory space. 
sum_of_red=sum_of_red+int(bufferl[counter]); 

}//End of for(counter) 
avr_red=sum_of_red/10000; 
avarage_red_datas[index] [0]=avr_red; 
avarage_red_datas[index] [ 1 ]=x; 
avarage_red_datas[index] [2]=y; 
index=index+l; 

sum_of_red = 0.0; 

}//End of for(x).. 


} // End of for(y).. 


}// End of Read Buffer... 
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/****************************************************************************** 

* This function reads the display screen pixels from the frame buffer with a 100 by 100 * 

* square and finds the maximum and minimum illuminated pixel in that square and finds * 

* a contrast value for that square by subtracting the max from min illumination value 11 

* After finding the contrast values for all squares it picks up the max. contrast of the * 

* screen * 

******************************************************************************/ 

void Read_TV_Screen() 

{ 

int maxcontrast = 0,min_contrast = 255; 

int max_contrast_x ,max_contrast_y, 

min_contrast_x ,min_contrast_y; 
float Contrast! 110]; 


int index = 0; 


for (int y = 800 ;y>= 100;y=y-100) 

{ 

for(int x =100 ;x <= 900 ; x= x + 100) 

{ 

glReadPixels (x,y, 100,100,GL_LUMINANCE,GL_UNSIGNED_BYTE,buflfer2); 

for (int counter = 0;counter < 10000 ;counter++) 

{ 

if(int(buffer2[counter]) > max_contrast) 

{ 

maxcontrast = int(buffer2[counter]); 
maxcontrastx = x; 
max_contrast_y = y; 

} 

if(int(buffer2 [counter]) < mincontrast) 

{ 

mincontrast = int(buffer2[counter]); 
min_contrast_x = x; 
min_contrast_y = y; 

} 

}//End of for(counter) 


Contrastfindex] = maxcontrast- mincontrast; 
if (Contrastfindex] > Max_Contrast_Of_Scene){ 

Max_Contrast_Of_Scene = Contrastfindex]; 

MaxContrastX = maxcontrastx; 

MaxContrastY = max_contrast_y; 

}//End of if(...). 

// Initialize the values for comparison purposes. If they are not initialized we never find 
// the correct min and max values, 
maxcontrast = 0; 
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mincontrast = 255 ; 

index=index+l; 

} //Endoffor(x).. 

} // End of for(y).. 

}// End of Read_TV_Screen.... 


/****************************************************************************** 

* This is a convenience routine for finding the index value of global "avarage red datas" * 

* array which has the max red value in it. * 


int Find_Index_For_Max_Value() 

{ 

int index; 


for(int i = 0; i < 110 ;i++) 

{ 

if(avarage_red_datas[i][0] > MAXRED) 

{ 

MAXRED = avarage_red_datas[i][0]; 
index = i; 

} 


} 

return index; 

}//End of Find_Max_Value. 
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/********************************************* 5 ^************************* 

* This function prints LOCK ON to the screen with big fonts. * 

************************************************************************/ 


void Print Lock On (int xpos , int ypos) 

{ 

// Since we are going to write on a 2D screen we must use the Orthographic projection. 
glMatrixMode(GLPROJECTION); 
glPushMatrix(); 
glLoadldentityO; 

gluOrtho2D(0.0,1400.0,0.0,800.0); 


glMatrixMode(GLMODELVIEW); 

glPushMatrixO; 

glLoadldentityO; 

glPushAttrib(GL_COLORJBUFFER_BIT); 
glColor3f(0.1, 1.0, 0.5); 

// Do a setpoint for the lower lefthand coordinate of the text string. 
glRasterPos3f(xpos, ypos, 0.0); 
drawCharstringffontHandlel, "LOCKED ON!"); 
glPopAttribO; 

glMatrixMode(GLMODELVIEW); 

glPopMatrixf); 

glMatrixMode(GLPROJECTION); 

glPopMatrixf); 

} 


/************************************************************************ 


* This function prints the temperature equilant values of the colors on each color box on 

* bottom of the screen. 

^*>|^**********>l^>M*****************************************************/ 


void Print_Numbers_On_Color_Scale() 

{ 

glMatrixMode(GL_PROJECTION); 
glPushMatrixf); 
glLoadldentityO; 

gluOrtho2D(0.0,1400.0,0.0,800.0); 


glMatrixMode(GLMODELVIEW); 
glPushMatrixf); 
glLoadldentityO; 
glColor3f(1.0, 1.0, 1.0); 

glRasterPos3f( 175.0,63.0,0.0); 
drawCharstring(fontHandle2," 10"); 
glRasterPos3f(225.0,63.0,0.0); 
drawCharstring(fontHandle2, "15"); 
glRasterPos3f(275.0,63.0,0.0); 
drawCharstring(fontHandle2, "20"); 
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glRasterPos3f(325.0,63.0,0.0); 
drawCharstring(fontHandle2, "25"); 
glRasterPos3f(375.0,63.0,0.0); 
drawCharstring(fontHandle2, "30"); 
glRasterPos3f(425.0,63.0,0.0); 
drawCharstring(fontHandle2, "35"); 
glRasterPos3f(475.0,63.0,0.0); 
drawCharstring(fontHandle2, "40"); 
glRasterPos3f(525.0,63.0,0.0); 
drawCharstring(fontHandle2, "45"); 
glRasterPos3f(575.0,63.0,0.0); 
drawCharstring(fontHandle2, "50"); 
glRasterPos3f(625.0,63.0,0.0); 
drawCharstring(fontHandle2, "55"); 
glRasterPos3f(675.0,63.0,0.0); 
drawCharstring(fontHandle2, "60"); 
glRasterPos3f(725.0,63.0,0.0); 
drawCharstring(fontHandle2, "65"); 
glRasterPos3f(775.0,63.0,0.0); 
drawCharstring(fontHandle2, "70"); 
glRasterPos3f(825.0,63.0,0.0); 
drawCharstring(fontHandle2, "75"); 
glRasterPos3f(875.0,63.0,0.0); 
drawCharstring(fontHandle2, "80"); 
glRasterPos3f(925.0,63.0,0.0); 
drawCharstring(fontHandle2, "85"); 
glRasterPos3f(975.0,63.0,0.0); 
drawCharstring(fontHandle2, "90"); 


glMatrixMode(GL_MODELVIEW); 

glPopMatrix(); 

glMatrixMode(GLPROJECTION); 

gIPopMatrix(); 


} 
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/************************************************************************ 
* This function finds the total radiance for a given temperature. * 

************************************************************************/ 
float Find_Total_Radiance(float Temp) 

{ 

float total_radiance = 0.0,e = 1.0 ; 

total_radiance=e *sigma*Temp*Temp*Temp*Temp; 
return total_radiance; 

} 


/*********************+************************************************** 


* This function apply the atmospheric extinction to the signal and finds the equilant 

* temp, value at the sensor. 


*******«****************************************************************/ 


float Temp_At_Sensor(float total_radiance, float mu) 

{ 

float pas, Temp ,R = (Firing_Distance/100),e = 1.0 ; 


pas=(total_radiance/PI)*exp(-mu*R); 
Temp=sqrt(sqrt(pas/(e*sigma))); 
return Temp; 


/*****>^>|^*:M**<M******************************************************** 

* This function modifies the position of the box that we draw to simulate the scanning * 

* function of the sensor. 

>|^**************j|^******************************************************/ 

void Update_ROW_COLOMN(int Row, int Colomn) 

{ 

if(Row > 1000) 

{ 

ROW = 0; 

COLOMN = COLOMN -100 ; 
if(Colomn <= 0) 

{ 

COLOMN = 900; 

SEARCHEDONES = 1; 

}// End of iffColomn). 

}//End of if(Row). 

else ROW = ROW + 100 ; 

}// End of Update_ROW_COLOMN(..) 
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ft************************ **************** ******************************* 

* This function is called from "draw_the_scene" function if the sensor is IR type. * 
**************************************>m********************************/ 


void Call_IR_Functions() 

{ 


if(!SEARCHED_ONES) // Then the box we are drawing to simulate the scanning doesn't 

// completly scanned the whole display screen. So go ahead and 
// update the position of the box for the next frame. 

{ 


draw_a_rectangle(RO W, COLOMN); 
Update_ROW_COLOMN(ROW,COLOMN); 

>// End of if(! SEARCHED ONES). 


else // Yes we completed scanning the sceene.. 

{ 

if(!READ_ONES) // Then we haven't read the screen previously. 

{ 

Read_Screen (); 

INDEXJFOR_MAX_RED = Find_Index_For_Max_Value(); 

READONES = 1; 

} 

else // It means this screen reading is at least second reading. 

{ 

if(MAX_RED > 55.0) // 55 is our treshhold value for red pixels.lt is approximately 

// equial to 0.2 red value of the pixel where 1.0 being the max. 


{ 

if(MAX_X == OLD MAX X && MAX Y — OLD MAX Y) 

{ 

Print Lock On (400,150); 

} 

if(MAX_X <= 1000 && MAX_Y < 1000) // We are inside the display screen. 

{ 

draw_a_rectangle(int(avarage_red_datas[INDEX_FOR_MAX_RED][l]), 
int(avarage_red_datas[INDEX_FOR_MAX_RED] [2])+50); 


} 


}// End of if(MAX_RED > 55.0). 

}//End of else.... 

}//End of else. 

// Here we draw the color scale in two dimensional screen... 
if(MODE == SYCLIC SCALE) (DrawTheSyclicColorScale();} 
else {DrawTheColorScale(MODE); } 

Print_Numbers_On_Color_Scale(); 
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}//End of Call_IR_Functions() 


/************************************************************************ 
* This function is called from "draw_the_scene" function if the sensor is IR type. * 
************************************************************************/ 


void CallTVFunctionsf) 

{ 


if(!SEARCHED_ONES) // Then the box we are drawing to simulate the scanning doesn't 

// completly scanned the whole display screen. So go ahead and 
// update the position of the box for the next frame. 


draw_a_rectangle(ROW,COLOMN); 

Update_ROW_COLOMN(ROW,COLOMN); 

}// End of iff! SE ARCHED ONES). 

else // Yes we completed scanning the sceene.. 

{ 

if(!READ_ONES) // Then we haven't read the screen previously. 

{ 

Read_TV_Screen (); 
if(Max_Contrast_Of_Scene > 55.0) 

{ 

draw_a_rectangle(Max_Contrast_X,Max_Contrast_Y + 50); 
}// End of Max_Contrast_Of_Scene > 55.0). 


READONES = 1; 

} 

else 

{ 

if(Max_Contrast_Of_Scene > 55.0) 

{ 

PrintLockOn (400,150); 

draw_a_rectangle(Max_Contrast_X,Max_Contrast_Y + 50); 
)// End of Max_Contrast_Of_Scene > 55.0). 

}//End of else.... 

}//End of else. 

}//End of Call_IR_Functions() 


y*<i********************************************************************** 
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* This function modifies the color array with respect to the color mode sellected. * 
************************************************************************/ 
void Modify_Color_Array_With_A_Given_Mode(int Verticeno, int SELECTED MODE, double redvalue) 
{ 

switch(SELECTEDMODE) 

{ 

case GRAY SCALE: 

Color[Vertice_no][0]= red_value; 

Color[Vertice_no][l]= redvalue; 

Color[Vertice_no]t2]= red_value; 

Color[Vertice_no][3]= 1.0; 

break; 

case TEMPERATURE_SCALE: 

Color[Vertice_no][0]= redvalue; 

Colorf Verticeno] [ 1 ]= 0.0; 

Color[Vertice_no] [2]= (1 -redvalue); 

Color[Vertice_no][3]= 1.0; 

break; 


default: 

break; 

}//End of switch.. 

}//End of Modify_Color_Table_With. 
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/*>M*********************************************************************************** 

* This is materialsupport.C * 

* Main body of this ".C" file is taken from CS4202 course notes and some material * 

* definitions are added. For example originally in this file "emission type” was not * 

* included in definitions of the properties of the materials. Since sun has it's own emission * 

* We seperatly defined a function call for "SUN" material. A function to enable and disable * 

* the fog is also defined for convience. * 

*************************************************************************************/ 
// These functions are self-contained and can be used by you in your applications. 


#include <Xm/Xm.h> // Get the Motif stuff... 

//include <GL/GLwMDrawA.h> // We are going to use an OpenGL Motif Draw widget 

//include <GL/gl.h> // Get the OpenGL required includes. 

#include <GL/glu.h> 

//include <GL/glx.h> 

//include <iostream.h> // C++ I/O subsystem, 
include <math.h> 

//include <stdlib.h> // Get exit() function. 

//include "materialsupport.h" // Get the names of the available materials. 

#include "materialsupport_funcs.h" // Get material support function names. 


/************************************************************************* 
* The following function turns on the Light Model. 
************************************************************************/ 


void tumOnTheLightModelf) 

{ 


// A dim white for the lighting model ambient color. 
GLfloat lmodel_ambient[] = { 0.2, 0.2, 0.2, 1.0 }; 


// Set the global ambient light intensity. 

glLightModelfv(GL_LIGHT_MODEL_AMBIENT, lmodel ambient); 

// Set whether the viewpoint position is local to the scene or 
// whether it should be considered to be an infinite distance away. 

// We have selected an infinite viewpoint. 

glLightModeli(GL_LIGHT_MODEL_LOCAL_VIEWER, GL FALSE); 

// Say that we want two-sided lighting. 
glLightModeli(GL_LIGHT_MODEL_TWO_SIDE, GL TRUE); 

// Enable the lighting model (there is a glDisable()). 
glEnable(GLLIGHTING); 


} 
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/************************************************************************* 
* The following function turns off the Light Model. 
************************************************************************/ 


void tumOffTheLightModel() 

{ 


glDisable(GLLIGHTING); 


} 


^************************************************************************* 


* The following function turns on the Sun Light 
************************************************************************/ 


void tumOnTheLights() 

{ 


GLfloat light_ambient[] = { 0.5, 0.5,0.5,1.0 }; 

GLfloat light_diffuse[] = { 0.7, 0.7, 0.7, 1.0 }; 

GLfloat light_diflusel[] = { 0.7, 0.7,0.7, 1.0 }; 

GLfloat light_specular[] = { 0.5,0.5, 0.5, 1.0 }; 

GLfloat light_specularl[] = { 1.0, 1.0, 1.0,1.0 }; 

// The ending 0.0 means directional light at infinity. 
GLfloat light_position[] = { 0.0, 100.0, -100.0, 0.0 }; 

GLfloat light2_position[] = { 0.0, 0.0, -100.0,0.0 }; 

// Specify the ambient rgba for the lights. 
glLightfv(GL_LIGHTO, GLAMBIENT, light_ambient); 


// Specify the diffuse rgba for the lights. 
glLightfV(GL_LIGHTO, GL DIFFUSE, light_difluse); 


// Specify the specular rgba for the lights. 
glLightfv(GL_LIGHTO, GL_SPECULAR, light_specular); 

// The position of the light is moving. So we didn't specify here.. 
// glLightPosition(GL_LIGHTO,.); 


// Turn the light on ... 
glEnable(GLLIGHTO); 


} 
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y************************************************************************* 

* The following function turns off the Sun Light 
*************************************************************************/ 


void tumOffTheLights() 

{ 

glDisable(GL_LIGHTO); 

} 


/**************************++***+***************************************** 
* The following function enables the fog. 

*************************************************************************/ 


//Here we are enabling the fog so the objects in the far distance look blury.. 
void EnableFogffloat density) 

{ 

//Apply Extinction to Electro-optical Signal 

GLint fogmode; 

glEnable(GL_FOG); 

{ 

GLfloatfogcolor[4]={0.8,0.894,0.815,0.5}; 

fogmode=GL_EXP; 

glFogi(GL_FOG_MODE,fogmode); 

glFogfv(GL_FOG_COLOR,fogcolor); 

glFogf(GL_FOG_DENSrTY,density); 

glHint(GL_FOG_HINT,GL_NICEST); 

glClearColor(0.3 921,0.3921,0.3921,0.5); 

} 

}//End of EnableFog. 

/************************************************************************* 
* The following function disables the fog 

*** * * * * *** 9ft * * * *** * ************* * ****** ** ** ** ** * * *** * ** * ****************** I 


//Disable the fog here., 
void DisableFogO 
{ 

glDisable(GLFOG); 

} 

y******************************iM*****>M********************************** 

* Turn on a Material for a particular face. I(Prof. Zyda) collected a number of materials 

* I found and made them easily accessible via define constant. 

* Options for whichFace = GLJFRONT, GLJ3ACK & others... 

*******>m****************************************************************/ 
void tumOnMaterial(GLenum whichFace, int whichMaterial) 

{ 


// Set the appropriate material type... 
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switch(whichMaterial) 

{ 


case BRASS: 

{ 

// Here are the brass material values... 

GLfloat brass_ambient[] = { 0.35, 0.25, 0.1, 1.0 }; 

GLfloat brass_difluse[] = { 0.65, 0.5,0.35, 1.0 }; 

GLfloat brassspecularQ = { 0.65,0.5, 0.35, 1.0 }; 

GLfloat brass_shininess[] = { 5.0 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace, 

brassambient, 

brassdifluse, 

brassspecular, 

brassshininess); 

} 

break; 

case SHINYBRASS: 

{ 

// Here are the shinybrass material values... 

GLfloat shinybrass_ambient[] = { 0.25, 0.15, 0.0,1.0 }; 

GLfloat shinybrass_difluse[] = { 0.65, 0.5,0.35, 1.0 }; 

GLfloat shinybrass_specular[] = { 0.9,0.6,0.0,1.0 }; 

GLfloat shinybrass_shininess[] = { 10.0 }; 

// Make the shinybrass material calls. 

makeGLMaterialCalls(whichFace, 

shinybrass_ambient, 

shinybrass_diffuse, 

shinybrass_specular, 

shinybrass_shininess); 

} 

break; 

case PEWTER: 

{ 

// Here are the pewter material values ... 

GLfloat pewter_ambient[] = { 0.0, 0.0, 0.0,1.0 }; 

GLfloat pewter_diffuse[] = { 0.6, 0.55, 0.65, 1.0 }; 

GLfloat pewter_specular[] = { 0.9,0.9,0.95,1.0 }; 

GLfloat pewter_shininess[] = { 10.0 }; 

// Make the pewter material calls. 

makeGLMaterialCalls(whichFace, 

pewterambient, 

pewterdifihise, 

pewterspecular, 

pewter_shininess); 

} 

break; 
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case SILVER: 

{ 

// Here are the silver material values ... 

GLfloat silver_ambient[] = {0.4,0.4,0.4, 1.0 }; 
GLfloat silver_diffuse[] = { 0.3,0.3, 0.3, 1.0 }; 
GLfloat silverspecularf] = { 0.9,0.9,0.95,1.0 }; 
GLfloat silver_shininess[] = { 30.0 }; 

// Make the silver material calls. 
makeGLMaterialCalls(whichFace, 
silverambient, 
silverdifluse, 
silver_specular, 
silvershininess); 

} 

break; 

case GOLD: 

{ 

// Here are the gold material values... 

GLfloat gold_ambient[] = { 0.4,0.2, 0.0, 1.0 }; 
GLfloat gold_diffuse[] = { 0.9, 0.5, 0.0, 1.0 }; 
GLfloat gold_specular[] = { 0.7,0.7,0.0, 1.0 }; 
GLfloat gold_shininess[] = { 10.0 }; 

// Make the gold material calls. 
makeGLMaterialCalls(whichFace, 
gold_ambient, 
gold_diffuse, 
gold_specular, 
gold_shininess); 

} 

break; 

case SHADOW: 

{ 

// Here are the shadow values... 

GLfloat Ambientf] = { 0.0,0.0,0.0, 0.4 ); 
GLfloat Diffuse[] = { 0.0, 0.0, 0.0, 0.4 }; 

GLfloat Specular[] = { 0.0,0.0,0.0,0.4 }; 
GLfloat Shininess[] = { 0.0 }; 


// Make the shinygold material calls. 
makeGLMaterialCalls(whichFace, 
Ambient, 

Diffuse, 

Specular, 

Shininess); 


} 

break; 


case PLASTER: 

{ 
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// Here are the plaster material values... 

GLfloat plasterambientfl = { 0.2, 0.2, 0.2, 1.0 }; 

GLfloat plaster_difluse[] = { 0.95, 0.95, 0.95,1.0 }; 

GLfloat plaster_specular[] = { 0.0, 0.0,0.0, 1.0 }; 

GLfloat plaster_shininess[] = { 1.0 }; 

// Make the plaster material calls. 

makeGLMaterialCalls(whichFace, 

plasterambient, 

plasterdifluse, 

plasterspecular, 

plaster_shininess); 

} 

break; 

case REDPLASTIC: 

{ 

// Here are the redplastic material values ... 

GLfloat redplastic_ambient[] = { 0.3, 0.1,0.1, 1.0 }; 

GLfloat redplastic_difluse[] = { 0.5, 0.1, 0.1, 1.0 }; 

GLfloat redplastic_specular[] = { 0.45,0.45, 0.45, 1.0 }; 

GLfloat redplastic_shininess[] = { 30.0 }; 

// Make the redplastic material calls. 

makeGLMaterialCalls(whichFace, 

redplastic_ambient, 

redplastic_difiuse, 

redplastic_specular, 

redplastic_shininess); 

} 

break; 

case GREENPLASTIC: 

{ 

// Here are the greenplastic material values... 

GLfloat greenplastic_ambient[] = { 0.1,0.3,0.1, 1.0 }; 

GLfloat greenplastic_difluse[] = { 0.1,0.5,0.1,1.0 }; 

GLfloat greenplastic_specular[] = { 0.45, 0.45, 0.45, 1.0 }; 

GLfloat greenplastic_shininess[] = { 30.0 }; 

// Make the greenplastic material calls. 

makeGLMaterialCalls(whichFace, 

greenplasticambient, 

greenplastic_difluse, 

greenplasticspecular, 

greenplasticshininess); 

} 

break; 

case BLUEPLASTIC: 

{ 

// Here are the blueplastic material values ... 

GLfloat blueplastic_ambient[] = { 0.1,0.1,0.3, 1.0 }; 

GLfloat blueplastic_difluse[] = { 0.1, 0.1, 0.5, 1.0 }; 
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GLfloat blueplastic_specular[] = { 0.45, 0.45, 0.45, 1.0 }; 

GLfloat blueplastic_shininess[] = { 30.0 }; 

// Make the blueplastic material calls. 

makeGLMaterialCalls(whichFace, 

blueplastic_ambient, 

blueplastic_diffuse, 

blueplastic_specular, 

blueplasticshininess); 

} 

break; 

case RED: 

{ 

// Here are the red material values... 

GLfloat red_ambient[] = { 0.2,0.0,0.0, 1.0 }; 

GLfloat red_diEFuse[] = { 1.0,0.0,0.0,1.0 }; 

GLfloat red_specular[] = { 1.0, 0.0,0.0,1.0 }; 

GLfloat red_shininess[] = { 1.0 }; 

// Make the red material calls. 

makeGLMaterialCalls(whichFace, 

redambient, 

reddifluse, 

red_specular, 

red_shininess); 

} 

break; 

case WHITE: 

{ 

// Here are the white material values... 

GLfloat white_ambient[] = { 0.2, 0.2, 0.2, 1.0 }; 

GLfloat white_difluse[] = { 1.0, 1.0,1.0, 1.0 }; 

GLfloat white_specular[] = { 1.0,1.0, 1.0,1.0 }; 

GLfloat white_shininess[] = { 1.0 }; 

// Make the red material calls. 

makeGLMaterialCalls(whichFace, 

whiteambient, 

white_diffuse, 

whitespecular, 

white_shininess); 

} 

break; 

case t72mat0: 

{ 

GLfloat t72mat0_ambient[] = { 0.039216,0.054902,0.039216}; 

GLfloat t72mat0_difluse[] = { 0.196078,0.274510,0.196078}; 

GLfloat t72mat0_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat0_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72matO_ambient ,t72mat0_diffiise ,t72mat0_specular 

,t72mat0_shininess); 
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} 

break; 

caset72matl : 

{ 

GLfloat t72matl_ambient[] = { 0.039216,0.078431,0.039216}; 

GLfloat t72matl_difiuse[] = { 0.196078,0.392157,0.196078}; 

GLfloat t72matl_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72matl_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72matl_ambient ,t72matl_diffuse ,t72matl_specular 
,t72matl_shininess); 

} 

break; 

case t72mat2 : 

{ 

GLfloat t72mat2_ambient[] = { 0.039216,0.039216,0.039216}; 

GLfloat t72mat2_diffiise[] = { 0.196078,0.196078,0.196078}; 

GLfloat t72mat2_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat2_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat2_ambient ,t72mat2_diffuse ,t72mat2_specular 
,t72mat2_shininess); 

} 

break; 

case t72mat3 : 

{ 

GLfloat t72mat3_ambient[] = { 0.039216,0.070588,0.039216}; 

GLfloat t72mat3_diffiise[] = { 0.196078,0.352941,0.196078}; 

GLfloat t72mat3_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat3_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat3_ambient ,t72mat3_diffuse ,t72mat3_specular 
,t72mat3_shininess); 

} 

break; 

case t72mat4 : 

{ 

GLfloat t72mat4_ambient[] = { 0.039216,0.062745,0.039216}; 

GLfloat t72mat4_diff\ise[] = { 0.196078,0.313725,0.196078}; 

GLfloat t72mat4_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat4_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat4_ambient ,t72mat4_difluse ,t72mat4_specular 
,t72mat4_shininess); 

} 

break; 

case t72mat5 : 

{ 

GLfloat t72mat5_ambient[] = { 0.047059,0.047059,0.047059}; 

GLfloat t72mat5_difluse[] = { 0.235294,0.235294,0.235294}; 

GLfloat t72mat5_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat5_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat5_ambient ,t72mat5_diffuse ,t72mat5_specular 
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,t72mat5_shininess); 

} 

break; 

case t72mat6 : 

{ 

GLfloat t72mat6_ambient[] = { 0.039216,0.094118,0.039216}; 

GLfloat t72mat6_difluse[] = { 0.196078,0.470588,0.196078}; 

GLfloat t72mat6_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat6_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat6_ambient,t72mat6_difluse,t72mat6_specular 
,t72mat6_shininess); 

} 

break; 

case t72mat7: 

{ 

GLfloat t72mat7_ambient[] = { 0.039216,0.086275,0.039216}; 

GLfloat t72mat7_difiuse[] = { 0.196078,0.431373,0.196078}; 

GLfloat t72mat7_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat7_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat7_ambient ,t72mat7_diffuse,t72mat7_specular 
,t72mat7_shininess); 

} 

break; 

case t72mat8 : 

{ 

GLfloat t72mat8_ambient[] = { 0.062745,0.062745,0.062745}; 

GLfloat t72mat8_difluse[] = { 0.313725,0.313725,0.313725}; 

GLfloat t72mat8_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat8_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat8_ambient,t72mat8_diffuse ,t72mat8_specular 
,t72mat8_shininess); 

} 

break; 

case t72mat9 : 

{ 

GLfloat t72mat9_ambient[] = { 0.039216,0.047059,0.039216}; 

GLfloat t72mat9_difluse[] = { 0.196078,0.235294,0.196078}; 

GLfloat t72mat9_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72mat9_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72mat9_ambient,t72mat9_difliise,t72mat9_specular 
,t72mat9_shininess); 

} 

break; 

case t72_3mat0 : 

{ 

GLfloat t72_3mat0_ambient[] = { 0.039216,0.078431,0.039216}; 

GLfloat t72_3mat0_difluse[] = { 0.196078,0.392157,0.196078}; 

GLfloat t72_3mat0_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72_3mat0_shininess[] = { 0.000000 }; 

// Make the brass material calls. 
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makeGLMaterialCalls(whichFace,t72_3matO_ambient ,t72_3matO_difluse ,t72_3mat0_specular 
,t72_3matO_shininess); 

} 

break; 

case t72_3matl: 

{ 

GLfloat t72_3matl_ambient[] = { 0.039216,0.094118,0.039216}; 

GLfloat t72_3matl_diffuse[] = { 0.196078,0.470588,0.196078}; 

GLfloat t72_3matl_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72_3matl_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72_3matl_ambient ,t72_3matl_difiuse ,t72_3matl_specular 
,t72_3matl_shininess); 

} 

break; 

case t72_3mat2 : 

{ 

GLfloat t72_3mat2_ambient[] = { 0.039216,0.062745,0.039216}; 

GLfloat t72_3mat2_difluse[] = { 0.196078,0.313725,0.196078}; 

GLfloat t72_3mat2_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72_3mat2_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72_3mat2_ambient ,t72_3mat2_diffuse ,t72_3mat2_specular 
,t72_3mat2_shininess); 

} 

break; 

case t72_3mat3 : 

{ 

GLfloat t72_3mat3_ambient[] = { 0.039216,0.070588,0.039216}; 

GLfloat t72_3mat3_difluse[] = { 0.196078,0.352941,0.196078}; 

GLfloat t72_3mat3_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72_3mat3_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72_3mat3_ambient ,t72_3mat3 diffuse ,t72_3mat3_specular 
,t72_3mat3_shininess); 

} 

break; 

case t72_3mat4: 

{ 

GLfloat t72_3mat4_ambient[] = { 0.039216,0.039216,0.039216}; 

GLfloat t72_3mat4_difluse[] = { 0.196078,0.196078,0.196078}; 

GLfloat t72_3mat4_specular[] = { 0.000000,0.000000,0.000000}; 

GLfloat t72_3mat4_shininess[] = { 0.000000 }; 

// Make the brass material calls. 

makeGLMaterialCalls(whichFace,t72_3mat4_ambient ,t72_3mat4_difliise ,t72_3mat4_specular 
,t72_3mat4_shininess); 

} 

break; 

case SUN: 

{ 

// Here are the pewter material values ... 

GLfloat SUN_ambient[] = { 1.0, 1.0, 1.0, 1.0 }; 

GLfloat SUN_difluse[] = { 1.0, 1.0, 1.0, 1.0 }; 
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GLfloat SUNspecularQ = { 1.0, 1.0, 1.0,1.0 }; 
GLfloat SUN_shininess[] = { 500.0 }; 

GLfloat SUN_emission[] = { 1.0, 1.0,1.0,1.0 }; 

// Make the pewter material calls. 
makeGLSunMaterialCallsfwhichFace, 
SUNambient, 

SUN_diffiise, 

SUNspecular, 

SUNshininess, 

SUN_emission); 

} 

break; 

default: 

break; 

} // end switch statement. 


} 


// Here is the function called from tumOnMaterial 
void makeGLMaterialCalls(GLenum whichFace, 
GLfloat *ambient, 

GLfloat *diffuse, 

GLfloat *specular, 

GLfloat *shininess) 


// Make the material calls. 

glMaterialfv(whichFace, GL AMBIENT, ambient); 
glMaterialfV(whichFace, GL DEFFUSE, diffuse); 
glMaterialfvfwhichFace, GL SPECULAR, specular); 
glMaterialfv(whichFace, GL SHININESS, shininess); 


} 


// Here is the function that makes the actual gl calls for sun material, 
void makeGLSunMaterialCalls(GLenum whichFace, 

GLfloat *ambient, 

GLfloat *diffuse, 

GLfloat * specular, 

GLfloat *shininess, 

GLfloat *emission) 


{ 


// Make the material calls. 
glMaterialfvfwhichFace, GL_AMBIENT, ambient); 
glMaterialfvfwhichFace, GL_DIFFUSE, diffuse); 
glMaterialfv(whichFace, GL SPECULAR, specular); 
glMaterialfv(whichFace, GL SHININESS, shininess); 
glMaterialfv(whichFace, GL_EMISSION, emission); 

} 
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^*******^****:M********************************************************* 

* This is makeRasterFont.C * 

* * 


* This file contains two functions to support the reading of X fonts * 

* as raster bitmaps for use with OpenGL character display * 

*******j^*i^***>M*******************************************************/ 

#include <Xm/Xm.h> //Get the Motif stuff... 


#include <GL/GLwMDrawA.h> // We are going to use an OpenGL Motif Draw widget 

#include <X1 l/StringDefs.h> 

#include <Xll/keysym.h> 

#include <GL/gl.h> // Get the OpenGL required includes. 

#include <GL/glu.h> 

#include <GL/glx.h> 

#include <iostream.h> // Get the C++ I/O functions. 


#include <stdlib.h> // Get the exit() function. 


/* * * ******** ** ****** ***** ****** * ******************************** ** ******** 


* makeRasterFont * 

* -- make a set of display lists for a particular X font. * 

* - returns base of lists for use in drawing fonts. * 

************************************************************************/ 


GLuint makeRasterFont/Widget w, char* fontname) 

{ 


XFontStruct *fontInfo; // ptr to an X font. 

Font id; // Font id. 

unsigned int first, last; // first and last char of the font. 

GLuint base; // The base of the display lists. 

// See if we can load the X font... 

fontlnfo = XLoadQueryFont(XtDisplay(w), fontname,); 

// Did we get a font? 
if(fontInfo — NULL) 

{ 

cerr «"makeRasterFont: font not found = " « fontname « endl; 
exit(0); 

} 

// Get the font id. 
id = fontInfo->fid; 


// Get some ptr info for the font, 
first = fontInfo->min_char_or_byte2; 
last = fontInfo->max_char_or_byte2; 

// Compute the display lists for the font, 
base = glGenLists(last+l); 
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if(base — 0) 

{ 

cerr «"makeRasterFont: out of display list."« endl; 
exit(0); 

} 

// Create bitmap display lists from an X font. 
glXUseXFont(id, first, last-first+1, base+first); 

// Return base 
return base; 


} // end of makeRasterFont 


/****************************************************************************** 

* * 


* drawCharstring(GLuint base, char *s) 

* — draw the char string with the specified set of display lists. 

* — base is the value returned from makeRasterFont, essentially the base of a set 

* of display lists. 


* * ******** ** * * **** * **** **** * * ** * * * ****** ** ***** * *** ** *** * 4c * * * * ** * * ***** * ** **** 


void drawCharstring(GLuint base, char *s) 

{ 


// Save the List Base Setting. 
glPushAttrib(GL_LIST_BIT); 

II Specify the display list base for CallLists. 
glListBase(base); 

// Draw the chars... 

glCallLists(strlen(s), GL UNSIGNED BYTE, (unsigned char *)s); 

// Restore the List Base setting. 
glPopAttribO; 


} 
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C**************************************** 3^**************************** 

* This is file NPSimage.C * 

* This file contains the methods that implement the NPSimage class. * 

***** + ******^******^********* + ************* ****************************/ 


#include "NPSimage.h" 


#include <string.h> 


#include <iostream.h> 

// Get the I/O stuff. 

#include <GL/gl.h> 
#include <GL/glu.h> 

// Get the SGI routines. 

#include "image.h" 

// SGI image structures. 


// External def for SGI image routines. 

extern "C" IMAGE *iopen(const char *filename, char *access,...); 
extern "C" void getrow(IMAGE *image, unsigned short *buf, int y, int z); 
extern "C" void putrow(IMAGE *image, unsigned short *buf, int y, int z); 
extern "C" void iclose(IMAGE *); 


// temp arrays to hold scratch info for processing an sgi image. 
// RGBA order for the indices, 
unsigned short buf[4][4096]; 


// Here is a constructor with no specs. 
NPSimage: :NPSimage() 

{ 


// If an NPSimage is previously defined, delete it. 

// NPSimage: :~NPSimage(); 

// We haven't been given an image size. 

// Create a default image. 
size[0] = 8; 
size[l] = 8; 
size[2] = 4; 

// Pt the image data at a simple 8x8x4 image, 
imagedata = new unsigned int [size[0] * sizefl]]; 

// No name for the image yet. 
name = new char [8] ; 
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strcpyfname,"default"); 


} 


// Here is the destructor for an NPSimage. 
NPSimage: :~NPSimage() 

{ 


// Call the del() function. 

delO; 


} 


// Delete the image, 
void NPSimage: :del() 
{ 


// Delete the image data if it exists. 
if(imagedata != 0) 

{ 

delete imagedata; 

} 

// Delete the image's name if it exists. 
if(name != 0) 

{ 

delete name; 

} 


} 


// Create an empty image of a particular size. 

// We pass in the name and size of the image to create. 
NPSimage: :NPSimage(char *name_img, int *size_img) 
{ 


// If an NPSimage is previously defined, delete it. 
NPSimage: :~NPSimage(); 

// Create an image. 

// Copy the size. 
size[0] = size_img[0]; 
size[l] = size_img[l]; 
size[2] = size_img[2]; 
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// Allocate an image of the appropriate size, 
imagedata = new unsigned int [size[0] * size[l]]; 

// No name for the image yet. 

name = new char [strlen(namejmg) + 1]; 

strcpy(name, name_img); 


} 


// Create an empty image of a particular size. 

// We pass in the name and size of the image to create. 
NPSimage::NPSimage(char *name_img, int xsize, int ysize, int zsize) 
{ 


// If an NPSimage is previously defined, delete it. 
NPSimage: :~NPSimage(); 

// Create an image. 

// Copy the size. 
size[0] = xsize; 
size[l] = ysize; 
size[2] = zsize; 

// Allocate an image of the appropriate size, 
imagedata = new unsigned int [size[0] * size[l]]; 

// No name for the image yet. 

name = new char [strlen(name img) +1]; 

strcpy(name, name img); 


} 


// We need to read in an image from a particular file, 
void NPSimage::read_from(const char *filename) 

{ 

IMAGE *image; // a ptr to an SGI image structure. 

int x,y,z; // Loop temps. 

unsigned int *ptr; // Ptr to the image longs. 


// open an sgi image 
image = iopen(filename,"r">; 
if(image == NULL) 

{ 
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cerr «"NPSimage::read_from: can't open input file" « filename « endl; 
return; 

} 

// If an NPSimage is previously defined, delete it. 

NPSimage: :~NPSimage(); 

// We have the image open and the old image deleted... 

// Create an empty image of the right size. 

// Copy the size. 
size[0] = image->xsize; 
size[l] = image->ysize; 
size[2] = image->zsize; 

// Allocate an image of the appropriate size, 
imagedata = new unsigned int [size[0] * size[l]]; 

// Copy the filename into the image, 
name = new char [strlen(filename) +1]; 
strcpy(name, filename); 


// get a pointer to the NPSimage longs, 
ptr = imagedata; 

// for each row of the image. 
for(y=0; y < size[l); y=y+l) 

{ 


for(z=0; z < size[2]; z=z+l) 

{ 


// Read each row in reds, greens, blues, alpha order. 
getrow(image, buflz], y, z); 


} 

// we must now step across the row and set each long integer 
// of the NPSimage format by combining the info from the sgi 
// rows. 

for(x=0; x < size[0]; x=x+l) 

{ 


// Someday we should generalize the below code to make it useable 
// for images with only 1 and 2 channels. 

// compute the RGBa long to plug in and plug it in. 
if(size[2] = 4) 

{ 

// Data reversed for OpenGL implementation. 

*ptr = (bufI0][x] « 24) + (buf[l][x] « 16) + (bufl2][x]« 8) + (bufl3][x]); 

} 
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else 

{ 

// RGBA image with alpha forced to Oxff 

// Originally: *ptr = bufiOJ[x] + (bufll][x] « 8) + (buf[2][x] « 16) + (Oxff « 24) 
// Images are stored in memory differently in OpenGL (RGBA) instead of 
// ABGR as in IRISGL. 

*ptr = (buf[0][x] « 24) + (buf[l][x] « 16) + (buf[2][x] « 8) + (Oxfi); 

} 

// step the ptr to the next long. 
ptr++; 


} 


} 

// Close the image file. 
iclose(image); 


} 


// We need to write out the image in SGI format, 
void NPSimage::write_to(const char ^filename) 

{ 

register IMAGE *image; // a ptr to an SGI image structure. 

int x,y,z; // Loop temps. 

unsigned int *ptr; // Ptr to the image longs. 


// open an sgi rgb image for writing. 

image=iopen(filename,"w",RLE(l),size[2],size[0],size[l],size[2]); 

// get a pointer to the NPSimage longs, 
ptr = imagedata; 

// for each row of the image... 
for(y=0; y < sizefl]; y=y+l) 

{ 


// we must now step across the row and decode each integer 
// of the NPSimage format into the 16 bit shorts sgi requires. 
for(x=0; x < size[0]; x=x+l) 

{ 


// Get the colors from the longs. Reversed for OpenGL. 
buf[0][x] = (unsigned short)((*ptr & OxfiTOOOOOO)» 24); 
buf[l][x] = (unsigned short)((*ptr & OxOOffOOOO)» 16); 
buf[2][x] = (unsigned short)((*ptr & OxOOOOffOO)» 8); 
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buf[3][x] = (unsigned short)(*ptr & OxOOOOOOfl); 


// step the ptr to the next long. 
ptr++; 


} 

// Write out the rows of the image. 
for(z=0; z < size[2]; z=z+l) 

{ 


// Write each row in reds, greens, blues, alpha order. 
putrow(image, buflz], y, z); 


} 


} 

// we must close the output sgi image file. 
iclose(image); 


} 


// We need to display the image, 
void NPSimage::displayO 
{ 

// Here is the OpenGL call to send the image to the open window. 
glDrawPixels(size[0], sizefl], GL_RGBA, GLUNSIGNEDBYTE, imagedata); 


} 


// The purpose of the = operator is for copying an NPSimage. 
NPSimage& 

NPSimage: :operator=(NPSimage &img) 

{ 


inti; //Looptemp 


// Delete the old image, if any. 

NPSimage: :~NPSimage(); 

// Copy the size. 
size[0] = img.size[0]; 
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size[l] = img.sizefl]; 
size[2] = img.size[2]; 


// Allocate an image of the appropriate size, 
imagedata = new unsigned int [sizefO] * size[l]]; 

// Copy the image data from img to "this". 
for(i=0; i < (size[0] * size[l]); i=i+l) 

{ 

imagedatafi] =img.imagedata[i]; 

} 

// Copy the filename into the image, 
name = new char [strlen(img.name) +1]; 
strcpy(name, img. name); 

// Return a ptr to the new image, 
return *this; 


} 


// The purpose of the is_equal function is for comparing 2 NPSimages. 
// We return 1 if the images are equal, otherwise 0. 
int NPSimage::is_equal(NPSimage &img) 

{ 


inti; //Looptemp 


// Compare the sizes. 

if(size[0] != img.size[0] j| size[l] != img.sizefl] 
|| size[2] != img.size[2]) 

{ 

// The sizes are not equal, return 0. 
return 0; 

} 

// Compare the image data. 
for(i=0; i < (size[0] * sizefl]); i=i+l) 

{ 

if(imagedata[i] != img.imagedata[i]) 

{ 

// There is a difference between the images, 
return 0; 

} 


} 


// If we get here, the images are equal. 

// Note: we don't check to see if the names are equal, 
return 1; 
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// The purpose of the colorreplace function is for 
// changing all occurrences of the old color to the new color, 
void NPSimage::color_replace(unsigned int oldclr, unsigned int newclr) 
{ 


inti; //Looptemp 


// Compare the image data. 
for(i=0; i < (size[0] * size[l]>; i=i+l) 

{ 

if(imagedata[i] == oldclr) 

{ 

// We have the old color, replace it. 
imagedata[i] = newclr; 

} 

} 
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I******************************************************************************* 


* This is TextureDisplayList.C 

* 

* This program is written to support the texturing need of the software Here we 

* used the display list phenomena to make the switching of texturing faster. Without 

* display list it is almost imposible to switch in a short time period though we developed 

* this program under Reality Engines. This Program is also very flexable when more 

* textures are needed simply call H makeGLTextureFromImagefile("???.rgb") M function 

* with necessaiy "rgb" file. 

* 

* 

* 


* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


* Authors: Ltjg. Mustafa YILMAZ 

* Ltjg. Mehmet GORGULU 

* 

****:M*****#******************************************************************/ 


* 

* 

* 


include <GL/gl.h> // Get the OpenGL required includes, 

^include <GL/glu.h> 

#include <GL/glx.h> 

#include "TextureDisplayListh" 

#include "texturesupport_funcs.h" // Get the texture support functions. 
#include "drawsupport_funcs.h" // Get the drawing functions. 
#include "eotda_globals_for_all_programs.h" 


void CallDisplayListForCube() 

{ 

CUBE = glGenLists(l); 
glNewList(CUBE,GL_COMPELE); 
drawCube(0.0,0.0, 0.0, 17998.0); 
glEndList(); 

} 


void CallDisplayListForTexturingO 

{ 

Index = glGenLists(ll); 

glNewList(Index,GL_COMPILE); 

makeGLTextureFromImagefile(”sea.rgb");glEndList(); 

glNewList(Index + l,GL_COMPILE); 

makeGLTextureFromImagefile("carpetl.rgb");glEndList(); 


glNewList(Index + 2,GL_COMPILE); 

makeGLTextureFromImagefile("cement4.rgb");glEndList(); 


glNewList(Index + 3,GL_COMPILE); 

makeGLTextureFromImagefile("sand3.rgb");glEndList(); 
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glNewList(Index + 4,GL_C0MPILE); 

makeGLTextureFromImagefile("brick3.rgb");glEndList(); 

glNewList(Index + 5,GL_C0MPILE); 

makeGLTextureFromImagefile("clouds2.rgb");glEndList(); 

glNewList(Index + 6,GL_C0MPILE); 

makeGLTextureFromImagefile("rocks.rgb");glEndList(); 

glNewList(Index + 7,GL_C0MPILE); 

makeGLTextureFromImagefile("window.rgb");glEndList(); 

glNewList(Index + 8,GL_C0MPILE); 

makeGLTextureFromImagefile("roof.rgb");glEndList(); 

glNewList(Index + 9,GL_C0MPILE); 

makeGLTextureFromImagefile("wood.rgb");glEndList(); 

glNewList(Index + 10,GL_COMPILE); 

makeGLTextureFromImagefile("city.rgb");glEndList(); 
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* This is "drawsupportjuncs.h" header file.This header file is included only * 

* in "drawsupport.C" file.. * 

***************************♦**„,*****„, ♦****%**********,|,**** J( , % **** 1(1+ * +1|t * J(1 ^ 

void drawTarget(int Target_no, float x,float y,float z); 

void drawGround(float x, float y, float z ,int Terrain,int Mode); 

void DrawTheSun(int TIME); 

void DrawTheColorScale(int MODE); 

void DrawTheSyclicColorScale(); 

void drawSphere(float x, float y,float z, 

float radius, int nslices, int nstacks); 

void draw_a_rectangle(int x,int y); 

void DrawTheMoutains(int Mode); 

void drawCube(float x, float y, float z, float sidelength); 

void drawTargetShadow(int Target_no, float x ,float y ,float z); 

void update_green_blue(double * RGB ARRAY, int MODE); 
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**************** **** ** ******^*************************************** * * 

* This is file eotdafuncs.h It is the header file for functions defined * 

* in eotdamain.C file. * 

* * 

* Authors: Ltjg. Mustafa YILMAZ * 

* Ltjg. Mehmet GORGULU * 

************************************************************************/ 


l^^t******************************************************************** 

* — CALBACKS— * 

*^*tti^t*****************************************************************/ 

static void rotatexCB(Widget, XtPointer user_data, XtPointer); 
static void rotateyCB(Widget, XtPointer userdata, XtPointer); 
static void rotatezCB(Widget, XtPointer user data, XtPointer); 
static void translatezCB( Widget. XtPointer user_data, XtPointer); 
static void targetCB(Widget, XtPointer user data, XtPointer); 
static void BackgroundlCB(Widget, XtPointer user_data, XtPointer); 
static void Background2CB(Widget, XtPointer user_data, XtPointer); 
static void Background3CB(Widget, XtPointer user data, XtPointer); 
static void Background4CB(Widget, XtPointer user data, XtPointer); 
static void SensorCB(Widget, XtPointer user data, XtPointer); 
static void initCB (Widget w, XtPointer, XtPointer); 
static void exposeCB(Widget w, XtPointer, XtPointer); 
static void resizeCB (Widget w, XtPointer, XtPointer call data); 
static void inputCB(Widget, XtPointer, XtPointer call_data); 
static void quitCB(Widget w, XtPointer, XtPointer); 
static void newviewCB(Widget, XtPointer, XtPointer); 
static void DisplayCB(Widget, XtPointer, XtPointer); 
static void BackToInputScreenCB(Widget, XtPointer, XtPointer); 
static void materialCB(Widget, XtPointer user data, XtPointer); 
static void TargetAreaCB(Widget, XtPointer user data, XtPointer); 



static void TargetHeadingCB(Widget, XtPointer user_data, XtPointer); 
static void TargetSpeedCB(Widget, XtPointer userdata, XtPointer); 
static void FireStatusCB(Widget, XtPointer user data, XtPointer); 
static void EngineStatusCB(Widget, XtPointer user data, XtPointer); 
static void IsolationStatusCB(Widget, XtPointer user data, XtPointer); 
static void RelativeApertureCB(Widget, XtPointer user data, XtPointer); 
static void ApertureDramCB(Widget, XtPointer user data, XtPointer); 
static void MagnificationCB(Widget, XtPointer user_data, XtPointer); 
static void OpticalTransmittanceCB(Widget, XtPointer user data, XtPointer); 
static void ElectronicBandwidthCB(Widget, XtPointer user_data, XtPointer); 
static void SensorHeightCB(Widget, XtPointer user data, XtPointer); 
static void DetectivityCB(Widget, XtPointer user data, XtPointer); 
static void DetectorElementsCB(Widget, XtPointer user data, XtPointer); 
static void FiringAngleCB(Widget, XtPointer user_data, XtPointer); 
static void FiringDistanceCB(Widget, XtPointer user data, XtPointer); 
static void ChangeWindSpeedCB(Widget, XtPointer user data, XtPointer); 
static void WindDirectionCB(Widget, XtPointer user_data, XtPointer); 
static void AeresolCB(Widget, XtPointer user data, XtPointer); 
static void ChangeTimeHourCB(Widget, XtPointer user data, XtPointer); 
static void ChangeTimeMonthCB(Widget, XtPointer user_data, XtPointer); 
static void ChangeTimeYearCB(Widget, XtPointer user_data, XtPointer); 
static void LongitudeCB(Widget, XtPointer user data, XtPointer); 
static void LatitudeCB(Widget, XtPointer user data, XtPointer); 
static void CloudinessCB(Widget, XtPointer user data, XtPointer); 
static void FogCB(Widget, XtPointer, XtPointer); 
static void FogDensityCB(Widget, XtPointer, XtPointer); 
static void IR_ScaleCB(Widget, XtPointer, XtPointer); 
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y********>M************************************************************** 

* —FUNCTIONS— * 

**#**#******************************************************************/ 


GLboolean drawWPO; 

void draw_the_scene(); 

void Create_The_Scales(Widget Parent); 

void BuildlnputScreenDialogO; 

void create_target_screen_widgets (Widget Parent); 

void createsensorscreenwidgets (Widget Parent); 

void createbackgroundscreenwidgets (Widget Parent); 

void create_others_screen_widgets (Widget Parent); 

void create_bulletin_boards(Widget); 

void concatenate_accumulative_matrices(float refx, float refy, float refz); 

void unit(GLfloat *m); 

void set_initial_viewing_valuesO; 

void compute_new_viewing_values(); 

void make_the_window_for_the_controls(); 

void create_toggle_widget(Widget parent); 

void load_the_viewpoint(); 

int FillInInitialTargetArrays(char * InputFileName); 
void make_pulldown_entry( 

Widget menupane, // A menupane of the menubar, 

char *entiytext, // Display text for the pulldown entry. 

XtCallbackProc callback, // Name of procedure to call 

// when this menu entry is selected. 
XtPointer user_data); // Data to be sent the callback. 


void Read_Screen(); 
void Read_TV_Screen(); 
float Calculate_Extinction(); 

int Fill_Up_Color_Table(int Vertice_no,int ColorMode); 
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void Modily_Color_Array_With_A_Given_Mode(int Vertice_no,int SELECTEDMODE, double 
redvalue); 

float Find_Total_Radiance(float kt); 

float Temp_At_Sensor(float totalradiance , float mu); 

int LoadModifiedGroundColor(float Range, float mu); 

int Find_Index_For_Max_Value(); 

void Update_ROW_COLOMN(int, int); 

void Print_Lock_On(int MAXX, int MAXY); 

void Call_IR_Functions(); 

void Call_TV_Functions(); 

voidPrint_Numbers_On_Color_Scale(); 
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#ifndef _EOTDA_GLOBALS. 
//define_EOTDAGLOBALS 


/♦♦♦♦fr******************************************************************* 

* Global widgets definitions * 

**************j^*****jM*************************************************/ 

Widget 

toplevel, // The top level (shell) widget. 

Main Window, // The topmost container widget (excluding the shell). 

// This type of widget is used as a container for a 
// collection of widgets that dynamically resize 
// together. 

InputWindow = 0, // This is the main widget that we use for input 

// window. 

Display Window = 0, // This is the main widget that we use for display 

// window. 

ControlWindow, 

DisplayButton, 
bulletinboardl, 
bulletin_board2, 
bulletin_board3, 
bulletin_board4, 

BboardForExit, 

BboardForReset, // Bulleten boards for reset and exit buttons. 

BboardForlnput, // B.B. for InputScreen button.. 

MainDisplayFrame, // Freame that holds the Display Window 

frame, // Frame widget attached to Form widget. 

framel, 

frame2, 

frame3, 

frame4, // Frame widget attached to bulletin board 

//widgets. 

glw, // The GL widget that we use to draw into. 

info, 

scalel, 

scale2, 

scale3, 

scale4, // Scale widget. 

controls, // Form widget for the toggle controls, 

rc, // The bulleten board widget for control window. 

rc2 , // The bulleten board for pushbuttons in display 

// screen 

sep, // The seperator widget. 
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SensorRCll, 

TargetRowColumn3, 

TargetRowColiuiui4, 

TargetBB5, 

TargetBB6, 

TargetBB7, 

FogBB, 

WindSpeedLabel, // The label Gadget holds the wind speed. 

TimeLabelHour, // The label Gadget holds the Hour 

TimeLabelMonth, // The label Gadget holds the Month. 

TimeLabelYear, // The label Gadget holds the Year. 

FiringAngleTextWidget, 

FiringDistanceTextWidget; 


/*********************************** ************************************* 
* WINDOW VARIABLES * 

********************** 4 ^************************************************/ 


static XtAppContext app_context; // Applications context, 
static XtWorkProcId workprocid = NULL; // Id of the work procedure 
static GLXContext glx_context; // A GL context (see initCB). 

Display *global_display; // Global display. 

Window global_window; // Global window. 

static XmStringCharSet charset = (XmStringCharSet) XmSTRINGDEFAULTCH ARSET; 

/************************************************************************ 

* CONSTANTS * 

#define sigma 5.6697*pow(1.0,-12.0) 

#define PI 3.1415926545 

#define Speed_Of_Light 3*pow(1.0,8.0) 

#define BoltzmanConstant 1.38*pow(l. 0,-23.0) 

#define he 1.986*pow(1.0,-25.0) 

#define IntegralConstant 1.1916*pow(1.0,-16.0) 

#define PT 100 
#define STROKE 101 
#define END 102 


I************************************************************************ 

* CLASS DEFINITIONS * 

************************************************************************/ 


class winddata { 
public: 


value 

direction ; 


int 

int 

}; 
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class timedata { 


public: 


int 

hour ; 

int 

month ; 

int 

year ; 


}; 

typedef struct charpoint { 
GLfloat x, y; 
int type; 

}CP; 


/************************************************************************ 
* VARIABLES * 

**************♦*********************************************************/ 


int 

n, i ; 

static float 

tz ; 

static float 

pitch ; 

static float 

heading; 

static float 

roll ; 

static float 

density; 

static float 

R_values[27][4] 

float MAXRED = 0.0, 

Max Contrast_Of_Scene 

int MAX_X = 0, 


// Used to specify number of args. 

// translation on in the z direction. 

// rotation on in the x direction. 

// rotation on in the y direction. 

// rotation on in the z direction. 

// used for fog density. 

= {0.0}; // This array is keeping the red color values 
// for each box which is 100 by 100 pixel. 


MAX Y = 0, 


Max_Contrast_X, 

Max_Contrast_Y, 

ROW = 0, 

COLOMN = 900, 

OLDMAXX = 1, 

OLDMAXY = 1, 

INDEX_FOR_MAX_RED, 

READONES = 0, 

DISPLAY_ENABLED = 0, // We use it to control displaye screen functions. 

SEARCHEDONES = 0, 

BackGround = 35, // Globally set background foR default value which is SOIL 

Target = 29, // Set default target value which is TANK 

HOUR = 0, // Used in positioning the sun function. 

MODE = 202 ; // Default scale is temperature scale. 


GLint FOG_FLAG; 
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/************************************************************************ 

* VARIABLES FOR INPUT SCREEN * 

**********>m#***********************************************************/ 


static float 


Firing Distance = 400.0,// Default distance 

FiringAngle = 10.0 , // Apply a default value (in degrees) 

// measured from +y axis towards x plane. 


Target_Area, 

TargetJHeading = 0, //Heading from North in degrees 
Target_Speed, 

RelativeAperture, 

Aperture_Dram, 

Magnification, 

Optical_Transmittance, 

Eln_Band_ Width, 

SensorHeight, 

Detectivity, 

Wind_Speed, 

WindDirection, 

Latitude, //Lat and Long, values would better be 

Longitude; //defined in a class so degrees and 

// minutes values could be stored seperately. 


/************************************************************************ 
* VARIABLES FOR TARGETS * 

************************************** 4 ^********************************/ 


// Globals for targets status. 

int ENGINE ON = 1, // For Tank Targets. 

FIRED =1, //diddo. 

HEATJSOLATED = 1, // For buildings. 

PARKED = 1 ; // For Aircraft. 


// These variables are used to control loop control variables upper limit, which 
// is different for different targets. 

int Tanklndex = 1384 , 

Airplane_Index = 208 , 

Shiplndex = 629, 

Buildinglndex = 38 , 

VERTICE_NUM = Tank_Index; //Default value for num_of_vertices. 


#endif 



/************************************************************************* 
* * 

* This is eotdaglobals.h file.Here all necessary global variables and constants are * 

* defined which could be refered in all .C files * 

* * 

* Authors: Ltjg. Mustafa YILMAZ * 

* Ltjg. Mehmet GORGULU * 

* * 

********+***************************************************************/ 
#ifndef _COMMON_ 

#define_COMMON_extern 

#endif 

_COMMON_ 

float Vertex[2500] [3], // Array to hold the coordinates. 

Color [2500] [4], // Array to hold RGB Alpha values. 

Normal[2500][3], // Array to hold Normal array values. 

Temperature [2500][1], // Array to hold Temperature for each veretex. 
Original_Temp[2500][l], // Array to hold initial temp. 

Viewpoint[3], //-x-y-z coordinates for viewer. 

Refpoint [3], //x-y-z coordinates for the point looked towards. 

avarage_red_datas[110][3], 

GroundColor[4] ; 

__COMMON_ 

GLuint Index, 

BACKGROUND, 

MOUNTAIN , 

CLOUDS , 

CUBE ; 

_COMMON_ int SENSOR = 39; //Default sensor is IR. 

#define IR39 

#define TV 40 

#define LASER 41 

#define GRAY SCALE 200 

#define SYCLIC_SCALE 201 

#define TEMPERATURESCALE 202 

static int _IR = IR; 

static int _TV = TV; 

static int LASER = LASER; 

static int _GRAY_SCALE = GRAY SCALE ; 

static int _SYCLIC_SCALE = SYCLIC_SCALE ; 

static int TEMPERATURE SCALE = TEMPERATURE SCALE ; 


163 



/************************************************************************ 

* This is "global_targets.h " header file. Here are the global definitions for targets and * 

* environment textures. * 

************************************************************************/ 


#ifndef GLOBALTARGETSH 

# define GLOB ALT ARGETSH 

#define TANK 29 

^define TRUCK 30 

#define SHIP 100 

#define HELO 32 

#define HOUSE 33 

#define PLANE 34 

#define SOIL 35 

#define GRASS 36 

#define SEA 37 

#define ASPHALT 38 


static int TANK = TANK; 
static int _TRUCK = TRUCK; 
static int SHIP = SHIP; 
static int HELO = HELO; 
static int HOUSE = HOUSE; 
static int PLANE = PLANE; 
static int _SOEL = SOIL ; 
static int GRASS = GRASS ; 
static int SEA = SEA; 
static int ASPHALT = ASPHALT; 

#endif 
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#ifndef _GL_MAGE_H 
#define _GL_MAGE_H 

#ifdef_cplusplus 

extern "C" { 

#endif 


/* 

* Defines for image files.... 

* 

* Paul Haeberli -1984 

* Look in /usr/people/4Dgifts/iristools/imgtools for example code! 

* 

*/ 


include <stdio.h> 

#define MAGIC 0732 

/* colormap of images */ 

#define CMNORMAL 0 


#define CMDITHERED 1 

#define CM SCREEN 2 


#define CM_COLORMAP 

#define TYPEMASK 
#define BPPMASK 
#define ITYPEVERBATM 
#define ITYPERLE 
Adeline ISRLE(type) 

#define ISVERBATM(type) 
Adefine BPP(type) 

#define RLE(bpp) 

#define VERBATM(bpp) 
#define IBUFSIZE(pixels) 
#define RLE NOP 


/* file contains rows of values which 

* are either RGB values (zsize — 3) 

* or greyramp values (zsize ==!)*/ 

/* file contains data which is a screen 

* image; getrow returns buffer which 

* can be displayed directly with 

* writepixels */ 

3 I* a colormap file */ 


OxffOO 

OxOOff 

0x0000 

0x0100 

(((type) & OxffOO) == ITYPE RLE) 

(((type) & OxffOO) = ITYPE_VERBATM) 
((type) & BPPMASK) 

(ITYPE_RLE | (bpp)) 

(ITYPE VERBATM | (bpp)) 
((pixels+(pixels»6))«2) 

0x00 


#define ierror(p) (((p)->flags&_IOERR)!=0) 

#define ifileno(p) ((p)->file) 

#define getpix(p) (~(p)->cnt>=0 ? *(p)->ptr++ : ifilbuf(p)) 

#define putpix(p,x) (~(p)->cnt>=0 \ 

? ((int)(*(p)->ptr++=(unsigned)(x))) \ 

: iflsbuf(p,(unsigned)(x))) 

typedef struct { 

unsigned short imagic; /* stuff saved on disk.. */ 


165 


unsigned short type; 

unsigned short dim; 

unsigned short xsize; 

unsigned short ysize; 

unsigned short zsize; 


unsigned long min; 
unsigned long max; 
unsigned long wastebytes; 
char name[80]; 

unsigned long colormap; 

long file; /* stuff used in core only */ 

unsigned short flags; 

short dorev; 

short x; 

short y; 

short z; 

short cnt; 

unsigned short *ptr; 

unsigned short *base; 

unsigned short ""tmpbuf; 

unsigned long offset; 

unsigned long rleend; /* for rle images */ 

unsigned long *rowstart; /* for rle images */ 
long ""rowsize; /* for rle images */ 

} IMAGE; 

IMAGE *icreate(); 

/* 

* IMAGE *iopen(char ""file, char "mode, unsigned int type, unsigned int dim, 

* unsigned int xsize, unsigned int ysize, unsigned int zsize); 

* IMAGE *fiopen(int f, char ""mode, unsigned int type, unsigned int dim, 

* unsigned int xsize, unsigned int ysize, unsigned int zsize); 

* 

* ...while iopen and fiopen can take an extended set of parameters, the 

* last five are optional, so a more correct prototype would be: 

* 

* IMAGE *iopen(char ""file, char ""mode,...); 

* IMAGE *fiopen(int f, char "mode,...); 

* 

* unsigned short *ibufalloc(IMAGE ""image); 

* int ifilbuf(IMAGE * image); 

* int iflush(IMAGE "image); 

* unsigned int iflsbuf(IMAGE "image, unsigned int c); 

* void isetname(IMAGE ""image, char ♦name); 

* void isetcolormap(IMAGE "image, int colormap); 

* int iclose(IMAGE *image); 

* 

* int putrow(IMAGE ""image, unsigned short *buffer, unsigned int y, unsigned int z); 

* int getrow(IMAGE ""image, unsigned short ""buffer, unsigned int y, unsigned int z); 

* 

*/ 


IMAGE *iopen(const char ""file, char ""mode,...); 
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IMAGE *icreate(); 
unsigned short *ibufalloc(); 

#defme IMAGEDEF /* for backwards compatibility */ 

#ifdef_cplusplus 

} 

#endif 

#endif /* !_GL_IMAGE_H_ */ 
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